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Die Spannungsverteilung in der Haut eines kugelférmigen Behilters bei 
tangentialer Krafteinleitung durch dquidistante Stiitzen in einem beliebigen 
Breitenkreis nach der Membrantheorie 


Von F.-J. Mainzer 


1. Einleitung. Im Schrifttum wird haufig darauf hingewiesen, daB die Membrantheorie nicht 
imstande sei, den Spannungszustand geschlossener Rotationsschalen unter Wahrung ihres geome- 
trischen Zusammenhanges bzw. unter Vermeidung von Singularitaten in der Spannungsverteilung 
zu beschreiben. Vielmehr miisse hierzu die Biegetheorie herangezogen werden. Im Falle einer 
idealisierten linienférmigen Einleitung von auBeren Kraften ist dies allerdings unumginglich, selbst 
dann, wenn die auBeren Krafte tangential zur Kugelhaut gerichtet sind. Es la8t sich jedoch nach- 
weisen, daf} die Membranspannungen in der Haut eines kugelférmigen Behalters unter Beanspru- 
chungen, wie sie bei der Speicherung von Gasen oder Fliissigkeiten unter Druck und bei geeigneter 
flachenhafter Einleitung der Stiitzkrafte entstehen, iiberall stetig sind und beschrankt bleiben. 

Die Anwendung der im folgenden zusammengestellten Ergebnisse auf ein rechnerisches Beispiel 
und die Gegeniiberstellung der Membranspannungsverteilung mit einem entsprechenden, aus der 
Biegetheorie hergeleiteten Spannungszustand soll in einer spateren Untersuchung erfolgen. 


2. Die Rotationsschale bei beliebiger Belastung. Die bekannten Differentialgleichungen der 
Membrantheorie fiir das Gleichgewicht der am Schalenelement einer beliebigen Rotationsflaiche 
angreifenden auBeren und inneren Krafte lauten!: 


ly . — ”) ae = r, Ngcosy +r,r, Ysing =0, 
Pp 
ON. 4 : ; 
Ty san “+ Te oa (Noo sing) +1 Nog csp +i, X sng =0, (1) 
rzN,—1 NO=—147,Z. 


Wie in Abb.1 veranschaulicht, bedeuten hierbei r, und r, den Meridian- und Ringkriimmungs- 
halbmesser, 9 den Winkel zwischen einem Meridianschnitt und einer festen, die Drehachse ent- 
haltenden Ebene, g den Winkel zwischen einer Schalennormalen und der Drehachse, X, Y, Z 
die Komponenten einer auBeren Kraft je Flacheneinheit an der Stelle (0, py), N,, No die inneren 
Normalkrafte je Langeneinheit des Breitenkreises bzw. des Meridians an der Stelle (0, y) und 
N,»» No, die inneren Schubkrafte je Langeneinheit des Breitenkreises bzw. des Meridians an der 


Stelle (0, ¢). 

a) Die Kugel bei rotationssymmetrischer Belastung. Im Falle einer Kugel werden 
die Kriimmungshalbmesser r, und r, einander gleich, und die weitere Voraussetzung der Rotations- 
symmetrie aller 4uBeren Kriafte bewirkt eine ebensolche der inneren Krafte, so daB alle Ableitungen 
nach %@ verschwinden. Damit vereinfachen sich die Differentialgleichungen (1) mit r;=rg=r au 


$i (N, sin ¢) a N,cosp+rYsing=0,  } 

3 (N,9 sing) + Ny, cosy +rXsing=0, (2) 
ay po P ao” 

ay -|- Ng =—=—r My . 


Bei der (im folgenden zumeist nur im Endergebnis aufgetiihrten) Integration der vorstehenden 
Gleichungen sind die auftretenden Integrationskonstanten jeweils so den gegebenen déuBeren 
Kraften anzupassen, da die inneren Spannungen iiberall in der Kugelhaut beschrankt bleiben. 


1s. 2. B. W. Fliigge, Statik und Dynamik der Schalen. 5. 21 ff. Berlin 1934. 
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Auf diese Weise kénnen im einzelnen die rotationssymmetrischen Beanspruchungen der Kugel- 
schale infolge Eigengewicht, Fiillung, Innendruck und Schneelast ermittelt werden!. Der Voll- 
stindigkeit halber seien sie hier, bezogen auf einen beliebigen, festen Unterstiitzungsbreitenkreis O; 
0<@O<ax mitgeteilt. 

1) Die Kugel unter Higengewicht. Die Komponenten der iiberall parallel zur Rotations- 
achse wirkend gedachten Flichenlast der Kugelschale tangential zum Meridian baw. senkrecht 
auf die Haut sind 


Cea, Z=y,t cos@, 


Y=y,tsing, 


ged dap 


avap 


—  ~y 


7 =73 Sinp 
af=-ndprap 


MN, pao. e Dlapa -r,r,snpdodd 


Abb. 1. Krifte am Schalenelement einer Rotationsfliche. 


wobei y, das spezifische Gewicht und t die Wandstarke der Kugelschale bedeuten. Sie erzeugen 
oberhalb des Unterstiitzungskreises ® die inneren Krafte je Langeneinheit 
a, ryit 


N. = 
? 1+ cos ’ 


1 fir 0<~<O@ (3) 


Mio, = =P OS 0D) = = ——— 
. YI nl a 1+ cos 
und unterhalb des Unterstiitzungskreises die Schnittkrafte 
ean t 
l1—cosp ’ 


i 
Ng =—ry, (cos gy + ora 


Ny, = 


fir, Oo = 7. (4) 


; 2) Die Kugel unter Fallungslast. Oberhalb des Fiillwinkels gy verschwinden mit den 
iuBeren Kraften auch die inneren. Unterhalb des Fillwinkels wird hydrostatischer Druck an- 
genommen, der nur senkrecht auf die Schale wirkt: 


X=Y=0, Z=—ryy (cos gy — cosq). 


1 Vel. W. Fliigge, a. a. O. S. 29f. 
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Er ruft oberhalb des Unterstiitzungskreises die Schnittkrafte 


ry 

No = Gain? g [2 008° Y — 008 Wo (3 cos’ p — cos? y)],) F 
fir gSps@ (5) 
Ng — Te Vu (cos Yo — cos 7) ans No 
und unterhalb des Unterstiitzungskreises die Schnittkrifte 

Tay 
No = % sint gp [3 008 Go sin” p + 2 (cos* p + ‘1 fir OD<g<n (6) 
No =P yx (cos Po — cos y) — 


hervor; y,, kennzeichnet hierbei das Caan Gewicht der Fiillung. 


3) Die Kugel unter Innendruck. Der Innendruck steht iiberall senkrecht auf der Kugel- 
haut und ist hoc 
A=Y=0, Z=—p. 


Demzufolge sind auch die inneren Normalkrafte je Langeneinheit in jeder Richtung einander 
gleich: 


N,=N,=yrp fir 0SpSa. (7) 
4) Die Kugel unter Schneelast. Wird fiir die Schneelast p, eine konstante Dicke der 
Schneedecke iiber der oberen Halbkugel zu Grunde gelegt, so haben die Lastkomponenten 
X=0, Y=p,sng, Z=—p, cos 
oberhalb des Aquators g =2/2 die Schnittkrafte 


(SDS Se 
x 1+ cos ’ © = 
l fiir ==. (8) 
No => — Pp; (cos Pp —— cacrd 


zur Folge. Zwischen dem Aquator und dem nicht dariiberliegend angenommenen Unterstiitzungs- 
kreis sind die Schnittkrafte 


r Ds 
N, =-= 
P sin? p ” 7 ; 
Ae igget ae cei (9) 
Nae : 
0 sin? p 


Unterhalb des Unterstiitzungskreises sind dann die aus der Schneelast stammenden inneren Krafte 
gleich Null. 

5) Die Schnittkraftdifferenzen im Unterstiitzungsbreitenkreis, Der Schnittkraft- 
verlauf infolge Eigengewicht, Fiillung und Schneelast weist in Héhe des Unterstiitzungskreises O 
fiir alle 9} jeweils einen konstanten Sprung auf. Wird eine jede Schnittkraftdifferenz wie folgt 


festgelegt: 

(AN, jo-o = (No sos — No@ses ®))p = (10a) 
und entsprechend 

(ANs)p—0 = (No @sos — Noeso<o))p=0> (10b) 
ergibt sich als Differenz der Meridianschnittkrafte infolge Kigengewicht, Fiillung, Innendruck und 
Schneelast der Reihe nach 


2ryyt 
(AN Jeno = 222, 
rey 
(AN,)o=0 = ; <a e (2 + 3 cos py — cos? Wo) (11) 
(AN,)¢ =@ 0 ? 
I Ps 
(ANo)a =o ee > ) 


Diese Schnittkraftdifferenzen JN, waren tangential an den zuniichst linienformig zu denkenden 
Stiitzring abzugeben. Ihre Vertikalkomponenten sind die je Laingeneinheit des Stiitzringes tiber- 
6* 


84 F.-J. Mainzer: Die Spannungsverteilung in der Haut eines kugelformigen Behalters Ingenieur-Archiv 


tragenen Lasten, die, mit dem Stiitzringumfang multipliziert, in der Summe das Gesamtgewicht 
der Kugel einschlieBlich Fiillung und Schneelast darstellen. 

Die Differenzen der Breitenkreisschnittkrafte stimmen auf Grund der dritten Gleichung (2) 
bis aufs Vorzeichen mit denen der Meridianschnittkrafte tiberein: 


(AN5)p = 0 = —(AN,)y=0- (12) 


b) Die Kugel bei aquidistanter Belastung. Im Hinblick auf die Ersetzung des umlau- 
fenden Stiitzringes durch eine Anzahl k= 2 von Einzelstiitzen, die tiber den Unterstiitzungsbreiten- 
kreis Aaquidistant verteilt sind, wird als Belastung in Meridianrichtung gewahlt 


Y= yy, cosnke. (13 a) 


Diese Gleichung gibt Veranlassung, die restlichen 4uBeren und die inneren Krafte in Abhangigkeit 
von # wie folgt anzusetzen!: 


Xa een | 
L =a COSmicys 
Ny = No, cosn kd, 
Np = Ny, cosnk®, 
y = Noo , Sin nk Vie 


Nog = Noo,, Sia oe 


(13b) 


An die Stelle von (1) treten somit drei neue Differentialgleichungen, die die Gleichgewichts- 
bedingungen fiir rotationsasymmetrische, periodisch verteilte auRere Lasten enthalten: 


ANo nt sin @) Sune Ny ; Y. sj 
r3— 7 ee aes 6, — 1 No, 8 P +1 Te np ng = 0, 
7) nit n 
(Ngo, Sin ~) 3 
k ; 
r,—— ap ——— =F, TURING tt Neg, COS @ eto Xy sino 0 (14) 
ne nk 


Die Elimination von Ny, . fiihrt im Falle der Kugel auf zwei voneinander unabhingige lineare 
T 
Differentialgleichungen erster Ordnung mit veranderlichen Koeffizienten und Stérfunktion fiir 
die Summe 


und die Differenz 


und zwar lJauten diese 


dU, ih \ 
cE, MANY rerd ee ao nk +- cos p 
dp 4 (2 tes + sin —p, Uink = Ni (Xu (er | sin p Lik ? 
U. 
~~ 2nk 2 et _ nk VAC nk—cosm, 
dp ( C £4 sin Pp Ua nk == if Xk Yo, | sin p as Zey: . 


Ihve Lésungen gewinnen nach Ausfiihrung einiger Integrationen die Gestalt 


\ nuk 
P 
(eww = 7 ; ; 
2 é - 7 ; nk 
Uh i tee St as | (Xn sie Blech Sis Ee a Soe Zu) sin (P (is 2) ay > 


sin? p i sin — 
(te a (15) 
- 2 


iz = : ; : nk —cos@p me /~ ee 
U, ne i. T I Xiak + Yur —- ance van sin? @ (cts xi a : 


1s. W. Fliigge, a. a. O. S. 37ff. 
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Um allgemeineren Belastungsfallen gerecht zu werden, sind die obigen Ansatze durch Summa- 
tion endlich oder unendlich vieler Glieder der Lasten zu erweitern: 


X= J) X,, sin nk?d, } 
H 


[o.9) 


Yo) YS Cosi. 
I 


(16a) 
Ss 

Z=)/2Z,,cosnk?d. 

Dementsprechend lauten die Ansitze der Schnittkrifte 


N. = J) Ny, cosnkd, } 
1 nk | 
No, cosnk?, | (16b) 


Noo = nk sin n k ) :J 


Hierbei stellen dann die zu einem bestimmten Index nk gehérenden Glieder der drei Fourier- 
Reihen jeweils die Lésung fiir die n k-te Harmonische der Lasten X,, sinnk#), Y,, cos nk 0, Z,, 
cosnk@® dar. 

Die in den Symmetriebedingungen der Ansatze (13) enthaltene Einschrinkung entfallt, wenn 
die gerade cos-Funktion durch die ungerade sin-Funktion ersetzt wird und umgekehrt. Dann 
ergeben sich neue Lésungsformeln fiir N,, Ny und Ny5,, die den alten iiberlagert werden kénnen, 
so daf sich die inneren Spannungen fiir jeden Fall einer periodischen Belastung berechnen lassen. 

1) Die Spannungsverteilung bei linienférmiger Stiitzung. Bei einer Ersetzung des 
umlaufenden linienférmigen Stiitzringes im Breitenkreis ® durch k Aquidistante und zunichst 
wiederum linienférmig anzunehmende Stiitzen von der Winkelbreite 2 y < 2 a/k sind den Lésungen 
der Differentialgleichungen (2) fiir die rotationssymmetrische Belastung zusatzlich Lésungen der 
Differentialgleichungen (14) fiir einen rotationsasymmetrischen Spannungszustand mit den Flachen- 
lasten X,,, sinnk} = Y,, cosnk? = Z,, cosnkd =0 fir jedes nk hinzuzufiigen; n = 1, 2, 3,...; 
k = Stiitzenzahl. Die Lésungen der Differentialgleichungen (14) fiir dieses zu ttberlagernde Span- 
nungssystem lauten dann unter Beriicksichtigung von (15) 


U, oe U. l nk nk 
= nk 7s Fas a Digi - vy : P 
Non 24 = = 2 sin” Y lee [eve 2 Sie Cok tg a Y 


U, = 1 ( i nk nk 
oe nk 2nk et UF aN tens 7a as 
No0nK eo Sa sin’ > [ote <2 4 Con '8 Z| : 


Diese Gleichungen gelten sowohl oberhalb als auch unterhalb des Unterstiitzungsbreitenkreises P 
jedoch mit verschiedenen Konstanten C, |, C, 
nk 


(17) 


Damit die Spannungen an den Polen der Kugel endlich bleiben, muf fiir die Kugelschale ober- 
halb des Unterstiitzungsbreitenkreises @ 


Come. (ian 23,8.) 
und fiir die Kugelschale unterhalb des Unterstiitzungsbreitenkreises 
Cpe eee Do a ae) 


71 Ic: 


nk 


in. = 
‘s Die restlichen Integrationskonstanten Gee und C,,, bestimmen sich aus der Forderung, daB 
das zu iiberlagernde Spannungssystem die Summe 3’ AN, der Meridianschnittkraftdifferenzen 
gema® den Gleichungen (11) zwischen den Stiitzen gerade aufhebt und in den Stiitzen selbst — 
entsprechend ihrer geringeren Tragbreite — erhéht. : 
Wie in Abb.2 (S. 88) veranschaulicht, mu dann aus Gleichgewichtsgriinden 
a—y 
Bae ager) te Tine Bi 
sein, wenn p, die senkrechte Last je Langeneinheit im umlaufenden Stiitzring, d. i die betas 
komponente der Summe )'AN,, der Meridianschnittkraftdifferenzen gemaB den Gleichungen (11), 
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Pp» die senkrechte Last je Lingeneinheit in den einzelnen Stiitzen und p;=p,—P, die der rotations- 
symmetrischen Lastverteilung zu tiberlagernde rotationsasymmetrische senkrechte Last je Langen- 


einheit in den Stiitzen bedeuten. 
Die Spriinge der zu iiberlagernden Spannungen betragen demnach, als Funktion p(0) betrachtet, 


fir v2a—y<0<v2a+y ; pO) = Ps =" Pi 


fir v2a+ty<0<(v+1)2a—y : pi) =—p,- 


Unter Verwendung der Fourier-Darstellung 


Pp) = Ds, aE sinnkycosnk (18) 


werden die zu iiberlagernden Spannungsspriinge im Unterstiitzungskreis 


= 1 ea Cit — @\nk 
> (A Non) _,cosnki} = Ven? p Od Cin (ete z) — @., (ts a cosnk} 


na 


CO 
1 Zp aes 
or Sek ; kd. 19a 
eae a as oe (19a) 
Da die Schubspannungen Ny, | im rotationssymmetrischen Spannungszustand tiberall gleich 
Null sind und im rotationsasymmetrischen Zustand keine Spannungsspriinge von Ny, , im Unter- 
n 


stiitzungsbreitenkreis zwischen den Stiitzen auftreten diirfen, mu auBerdem gelten 


. ; 1 4 @\nk = @\nk] , 
7a (Nee sinnk } = Pain? p pa [Cin (cts z) a Coe (te z) sin nk}=0. (19b) 
So ergibt sich fiir jedes n = 1, 2, 3,... ein inhomogenes Gleichungssystem von zwei Gleichungen 
mit den zwei Unbekannten C,,,, und C2,, und mit nichtverschwindender Determinante. Durch 
Kinsetzen seiner Lésungen in die Gleichungen fiir N,, Ny und Ng, lassen sich diese wie folgt 
schreiben: 


oS —P \ nk 
: ctg : 
; 2 Pie ae SUD EURO ORs Se 
sin? @ nky Bea; 
n=1 ctg a 
Pe te Q nk 
in @ 2 snnkycosnkd 
ONT ee eke aS GB DRS ELBE 
? sin? ae nky WEE OB 
n=l tg —— 
N»5 eae at: IN ? (20) 
nk 
5 ea ctg —— : 
_ py sin @ 2 snnkysnnk } 
te Dae 
=] } 
n= ctg ay 
ut eee ete ue ; : 
N,. —Pusin 2 snnkysinnk } 
oe sin? - D nky (WEST S 2) 
n= tg 2 


Die vorstehenden unendlichen Reihen kénnen mit Hilfe der nachfolgenden Beziehungen auch 
in geschlossene Ausdriicke verwandelt werden!: 


ie) . 
>) un See = aret u sin x 
n=1 ru AM ari, ju) S1 
oe und 
cos n x it 
ee aan 5 In (1 -2ucosx + u?), ucosx +1 


1's. z. B. Hiitte, Band I, 28. Aufl., S. 109 unter m) und k). Berlin 1955, 
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; 1 ee 
snnkycosnk& = sinnk (y + 8) + > sin n k (y — 8), 
sinnkysinnkj ——_ . : 
y =— a cosnk(y + 8) + x cosnk (y—%). 


po promempes in der Kugelhaut lauten dann, auBer an den Randern der linienférmigen 
u == i | 
zen (das sind die Stellen 9 = » 2 x ry; v=0,...,k—1)im Unterstiitzungsbreitenkreis D, 


| ip k k 


te tel ah 
; @ sin k (y + 2) 7 sin k (y — 8) 
N, _— Pisin ® ae) tes 
%~\<p<0 2h ysinty| TCS Ne I are tg : cae 
1—| —g@ | cosk(y + 8) 1—| —g | cosk (vy —9) 
- ae By 
ie . ke = 
tg tg 
=e sin k (y + 9) ars sin k (vy — #) 
N = Pan? |. ae: Be Tor 
°p<p<x = 2k y sin? aS 5 + arc tg : 
| ctg => ctg 2 
| 1—\ —~@ | cosk(y + 9) 1—| ——~g | cosk (vy — 8) 
1 Ce Ss. ctg co 
~ tel) 
ke 2k 
tg 4 tet 
T—2| —@ | cosk(y + 9) + 1g 
=a Pp, sin l So ee: 
9% <_<0 4ky sin? p y , Ok 9 \ 2k? 
Bo tg 5 
1—2 ~@ | ©os k(y— #8) + ae 
eo tg 
k 2k 
tg F etg F 
1—2(—g | cosky +0) +(—@ 
, me Pp, sin @ ks 2 ctg a 
? SQpss Ft Ak in2 ¢ 7 k : = _ ;, pias 
PS=PSx vy sin” ~ etg £ wee 
1—2 ——@ | 69s k (vy — #) +) — s 


Wahrend die beiden Gleichungen fiir N, zeigen, da an den fraglichen Stellen die zu iiber- 
lagernden inneren Normalspannungen beschrankt bleiben, lassen die Gleichungen fiir Ny, erkennen, 
daB die zu iiberlagernden inneren Schubspannungen an den Randern der linienférmigen Stiitzen 
tiber alle Grenzen wachsen. 


2) Die Spannungsverteilung bei flachenhafter Stiitzung. Die vorgenannten Singula- 
ritaten heben sich auf, wenn die linienférmige Stiitzung der Kugel im Unterstiitzungsbreitenkreis ®, 
wie in Abb.3 dargestellt, von einer flachenhaften Stiitzung mittels k Pratzenblechen der Winkelhéhe 
fp’ +f" und der Winkelbreite 2 y abgelést wird. In den Bereichen ® — f’ <= y < @ + Bf" mit 
y2a—ySdbsv20+y v=|l, ..., k) bewirken die Pratzenbleche dann einen tangentialen 
Flachenlastangriff Y= Y(0, gy) = 0 auf die Kugelhaut, deren Membranspannungszustand auf die 
Weise entwickelt werden kann, daB den bereits ermittelten Spannungen gemaB den Gleichungen (3) 
bis (9) und (21) zusatzlich Lésungen der Differentialgleichungen (14) tiberlagert werden. Diese 
miissen so beschaffen sein, da in den Stiitzen je Langeneinheit ihres Breitenkreisumfanges die 
Kraft p, auf die Pratzenbleche iibertragen und zugleich die von der linienformigen Stiitzung her- 
rihrenden Spannungsspriinge im Unterstiitzungsbreitenkreis ® aufgehoben werden. 


F.-J. Mainzer: Die Spannungsverteilung in der Haut ein 


es kugelformigen Behalters —_ Ingenieur-Archiv 
i) 


(9) betrachtet, in | 
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den Intervallen 
y2xa—yVd<v24a+y 
g(t) =O; 


Der Verlauf der in Abb.3 veranschaulichten Stiitzkrafte ist, als Funktion q 


Hierzu gehért die Fourier-Entwicklung 
[o.°) 
sinnkycosnkd = + dy) + > a,, cosnkd. 
et 


(o4 
QO) Pa eee 


y2Zaty<b<(y41)2a—y 
(22) 


f i 


qh Po=-G(H) 


Abb. 3. Flichenhafte Stiitzung mittels 
Pratzenblechen. 


Abb. 2. Die im Unterstiitzungsbreitenkreis ® zu iiberlagernden 
Spannungsspriinge. 

Die Verteilung von q(?) tiber der Winkelhéhe gy des Pratzenbleches richtet sich nach der Vor- 

schrift Y(0,~), die in der Praxis allerdings meist ungewif sein diirfte, und tiber die daher zunachst 


nicht weiter verfiigt werden soll, wenn nur beachtet wird, da®B gemaf Abb.4 


O48” 
{ 0X,¢8in od F = a8 sm: Or dyna 0, m2: 


(23) 


o—p' 
ist. Wegen dF = (rdp)(r sing d?) gilt ferner die spater bendtigte Beziehung 


D+p” 
Deeply Yo, Sin? @idgi ec sin (naa Ooaloe ce 
o—p’ 
Im Falle des konstanten Anteils a) der Fourier-Entwicklung (22) ergeben sich die Schnitt- 


» und Ny, mit Hilfe der Gleichungen (15) allgemein auf der ganzen Kugel zu 


krafte N,, N, 


1 
N, =——, 
sin* 


Noo = 0 . 
Oberhalb —/" und unterhalb @ + 6” ist Y)=0. Zusammen mit den Endlichkeitsbedingungen 
(24a) 


an den Polen der Kugel folgt also dort 
aN os 
Zur Festlegung der Integrationskonstanten in den Bereichen @— fi’ < pS und OS VSO@+f" 


ist zu beachten, daB an den Stellen ©—f’ und @+f"’ keine Spannungsspriinge vorhanden sein 
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diirfen. Demnach werden die Spannungen im Bereich @—f’ <p @ 


Ps 
r eet 
iN, =— Ng = Seca? @ | ve. sin ~p dp (24:b) 
o—p’ 
und im Bereich OS p< @ + fp” 
Dp tp” 
No =—No = = A | Y, sin? yp dp. (24c) 
pee 


Die itberall auf der Kugel giiltigen Lisungen der Differentialgleichungen (14) lauten im Falle 
der veranderlichen Summenglieder von (22) 


1 &2 nk 2 
DN eee Noy = Tae di (8 a [Cin ae A [> nk Sin? p (te ae ay) cos nk #) 
eo QP nk nk ; 
at Pm (te A |Com — rfy nk Sin? —p [erg H dp| cosnk O, 
N 1 ES QP nk (2 nk 
dp — (a) = Sa ee 59 Ink nk sin” W? a sinn kt 
= IN aint p Oy ("8 5) ic ie (ts 4 lp| sin nk 0 
oe nk 
-+- », [ts £) Cs, oe [¥ nk Sin? (et oe °| sinnk? 
n=1 E © 
Die Integrationskonstanten 
G. (t= 122) 


sind fiir die vier Bereiche 
0297S, 
o—pf so, 


O<pSO+p" sla 
und O+p"Sypsun 
im einzelnen wie folgt festzulegen: 
Oberhalb und unterhalb der Pratzen- 
bleche (das sind die Bereiche p-F 
ees tng 7-sin D-dB 
i +24 -sinpdF 
@+ fp" SpsSn) 
verschwinden mit den Flachenlasten Y,,, O+p” 


Abb. 4. Zusammenhang zwischen f : Var sin g dF und aK: 


aus Endlichkeitsgriinden auch die Inte- 
grationskonstanten Cj,, bzw. C, 


“nic * 
Zur Berechnung der noch unbekannten sechs restlichen Integrationskonstanten lassen sich vier 
Bestimmungsgleichungen aus der Forderung gewinnen, daB weder der Spannungsverlauf von N,, 


noch der von Ny, an den Stellen ® — f’ und ® + fi" unstetig werden darf: 


eel. 
? Pe are 
0<p<o—p o—p'<p<so : 
fir g=D—Pf', 
Nog = Noy 
0<~so-pf b—f'<pso 


== NN 
°B< pS 0+ 8" 2p +p <psan 
=/=> =lif B =P >! fiir p —@ As Bas 


9p <gso4p 4 B"S pS 
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Eine weitere Gleichung ergibt sich daraus, da in dem endgiiltigen iiberlagerten Spannungs- 
system die Meridianspannungen N, auch am Unterstiitzungsbreitenkreis @ stetig sein miissen. 
Das bedeutet, daB die neue, nun zu ermittelnde Spannungsverteilung die von der linienférmigen 
Stiitzung herrithrenden Spannungsspriinge im Unterstiitzungsbreitenkreis gemaB Gleichung (22) 
dort gerade aufzuheben hat: 


l aihor : 
AN ae [p ae 2 an sin n ky cos no) fur 7 =O. 


SchlieBlich miissen die Schubspannungen in der neuen Spannungsverteilung an der Stelle 

stetig sein: 
Mop <p <0) Noo <9< 0+") a ee 

Unter sehr allgemeinen Voraussetzungen ist die Darstellung einer Funktion durch eine Fourier- 
Reihe eindeutig, so da8 dann fiir jedes nk mit n = 1, 2, 3, ... ein inhomogenes Gleichungssystem 
von sechs linearen und unabhangigen Gleichungen mit nichtverschwindender Determinante in den | 
unbekannten Integrationskonstanten existiert. Mittels der eindeutigen Lésung kénnen also die 
Spannungen aus der flachenhaften Stiitzung auch im Falle der veranderlichen Summenglieder (22) 
angegeben werden. Die Uberlagerung mit den von der linienférmigen Einzelstiitzung [siehe die 
Gleichungen (20)] und den vom konstanten Anteil der flachenhaften Stiitzung [siehe die Glei- | 
chungen (24)] herriihrenden Spannungen ergibt im Bereich 0 < gy < O—f’: 


(oe) P+ Bp" 
r 7) nk meee nk 
Nv Tent p pa (ts 5) | Y, Sin? y (ote a dp cosnk?, 
- op c 
op" (25a) 
r nm og . 
Nog =— dan? i tg 4 | Y,, sin? y (cts =) dp sinnk?, 
i op 
im Bereich 6—f’ <p @ 
oo O + 
r nk F nk 
Ne Van? p Pp (ts 2 i Y 4 Sin? y (cts A dy cosnk? 
QP 
fore) P P 
nk nk 
—»>' (ote 4 | Y 4 Sin? y (ts =) dy cosnk ) —- 2 | Yo sin? y | : 
n=1 
of on’ 
n opr (25b) 
r n ales nk : 
Non = — ORY: oT (ts ZI | Y 4 Sin? p (cre a dp sinnk } 
~ 
+ D>) (ete z | Yup sin? @ (te 5) dy sinnk?, 
n=1 \ 
op’ 
im Bereich ODS pl@+H": 
a) +p" 
_ pe) Pp ae fe ae Y nk 
va a: > (ts z) | Y,, sim? (ots 5) dy cosnk } 
P 
COM Pp @+-p” 
nk 
es | | Y, Sin? (ts 2 dp cos nk + 2 | Y, sin? » in ' 
n=1 
op ? 
5 op" (25c) 
r o/ nk Sue nk * . 
No, = — Dain® \'s $| Y 4 Sin? p (ete 2] dp sinnk 
es @ 


@ 
nk i nk 
a y (ote 5) | Y 4 Sin? (ts 5) dp sinnk #}, 


o—Bp’ 
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im Bereich © + 8’ Spaz: 


lo) +e" 
<a r 2) nk kK gq \nk 
No ae 2 sin? Pp a (ete 9 | es sin? /P (ts £ dp cos n k 4). 9 
op’ 
oo Spa ’ (25d) 
fa pies ee Fm On aN : a 
Noo 2 sin? p aa (ete 9 ) [ Y yy Sin? p (ts z dy sinnk?. | 
o—p’ 


Bei naherem Zusehen erscheint es schwierig, notwendige und hinreichende Bedingungen fiir 
die Konvergrenz der auftretenden unendlichen Reihen in Kigenschaften der Funktionen Y,,, = Y,,,(@) 
auszudriicken. Jedoch ]aBt sich fiir den Fall beziiglich y konstanten Flachenlastangriffes Y= 3’ Y,, 
cosnk#? mit beziiglich y konstanten Y,, die gleichmaBige Konvergenz der unendlichen Reihen 
und die Beschranktheit aller Spannungen leicht nachweisen. 


Auf Grund der Gleichungen (22) und (23) mu8 dann 


Nt dy sin PB > _ ip, sin B 
9 D+” =3 +p ? 
r fj sin? » dp r f sin? y dp 
Ll =f (26) 
yank sin @ a 2pjsn®D snnky 
ot ae ss Pee. onky 
r f sin? » dp r Jf sin? ydy 
o—p' o—p' J 


o» 
nk 
( q | sin? (ote 4 dy , 
Oy 
und 
eo. k 
(ete at | sin? p (ts a dap 
é, 
mit 
0<6,56,56,.560,52 und O—f'’< 0, = 0, DO =p 
auszuwerten. ; 


Die exakte Berechnung dieser Integrale fiihrt auf Rekurionsformeln die, mit Y,,, und cos nk @ 
baw. sinnk # multipliziert, schwierig zu summieren sind. Wohl lassen sich die Integrale leicht ab- 
schatzen, wobei zwar brauchbare Abschatzungen der Héchstspannungen fiir jeden Breitenkreis 
zu gewinnen sind, der Verlauf der Spannungen in einem jeden Breitenkreis, als Funktion von ? 
betrachtet, jedoch verschleiert wird. 


Mit Hilfe der Substitutionen z = tg + baw. z = ctg folgt 
6 ace 2 tee 
2 - s Cc g es 
6.\nk ae nk 6.\nk gnk 2, 2; 
(ts Z) | se p [ete 5) dp < 2 (ts 2 t ear dz< ay ar 
’ 1 == (ets — 
2 ctg = 2, 
und entsprechend 
tg a 0) 
§.\nk i k 5 ices Bo 2 
(<te 2] | sin? p [tg z) dy S 2 (cts ) | ae da = ACS 
\ : eee ( Z) 
4 0; 2 
te = 
Zu den Normalspannungen NV, und N, = — N, gema® den Gleichungen (3) bis (9) und den 


Schubspannungen N,, = 0, die von den verschiedenen Belastungen bei linienformiger umlaufender 
Stiitzung herriihren, treten also noch folgende Héchstspannungen aus der flachenhaften Einzel- 


stiitzung hinzu, und zwar 
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im Bereich 0X gp S$ O— fi’: 


rps 
N 2 p, sin De Eig S qe 
Cia 2(b — B’) +8" / D+ Bp"? 6k? y? 
sin ( Oy i aie P ris 1 Lt Teles, = [ 
a 2 (27 a) 
pees 
2 Pp, sin D- ny ae 
Noy = = sin? (W— B’) FB" en mW 
; ie ra ae) aie sin? p oe, 1: (as? = p | 
im Bereich ® + BP’ Sp Sa: 
@P ae Be 
t J 
N = 2 pisinD , 2 
MO sin? 2(D —- Bb’) cOsse OFS! 5’ ) oS 
J sin?ydp 1 a a5 
o—B \ (27b) 
ee 
N 2 p, sin ® tg 2 
6g = = @+B) ote" BBN oR ) 
as a J sm? gdp I+ [ts ad J 
p—p' 


im Bereich 0 — ff’ <p <@: 


N. = — max ||N, , IN, a Y, sin? p do , 
: max {| a eee | ee LY, YP ap (276) 
See | lo WN, 
IN pea, Nee ee a > | PP e 4 pS psa) | : 


im Bereich OX p< @+ fp": 


i er 1 IN Bee hae e dy , 

N, Tec Me perry ee | Pp +B” es oe ee ey (27d) 
— | | 

Nop = + max { Nove ale 7 5 Noo er ; 


Wie die Gleichungen (25) zeigen, fallen die Spannungen aus der flachenhaften Einzelstiitzung 
von den Breitenkreisen ® — f’ und ® + f”’ gegen die Pole zu bei einer Stiitzenzahl k = 3 rasch 
nach Null und fiir & = 2 auf einen von Null verschiedenen kleinen Wert monoton ab. Im tbrigen 
hat die gleichmaBige Konvergenz der mit konstanten Y,,, gebildeten unendlichen Reihen zur Folge, 
daB in dem vollstandig tiberlagerten Spannungssystem alle Spannungen der Kugelhaut stetige 
Funktionen der beiden Veranderlichen # und sind. SchlieBlich sei darauf hingewiesen, daB sich 
die Meridianspannungen des tiberlagerten Spannungssystems bei der Krafteinleitung durch aquidi- 
stante Stiitzen keineswegs im Verhaltnis des vollen Stiitzringumfanges 2 rsin ® zur vermin- 
derten Tragbreite 2 ky r sin © erhéhen, sondern bei anzahlmafig hinreichend bemessenen Pratzen- 
blechen bzw. Stiitzen weit darunter liegen; das ist aus den Gleichungen (27) ohne Schwierigkeit zu 
ersehen. 


(Eingegangen am 7. Mai 1957.) 
Anschrift des Verfassers: Dip|.-Math. F.-J. Mainzer, Oberhausen-Sterkrade, Borbeckerstr. 35 
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Experimentelle Bestimmung der Wélbflache einer Vollkreismembran 
unter hydrostatischem Druck 


Von K.-E. Meier-Doérnberg 


1. Einleitung. Auf Anregung von Herrn Professor Karas habe ich zur Erginzung seiner kiirzlich 
veroffentlichten Arbeit! Die re und Kreisringmembranen unter hydrostatischem Druck, 
1. Mitteilung: Homogener Spannungszustand“ bei einer Vollkreismembran die durch hydrostati- 
schen Druck erzeugte Wélbflache ausgemessen. Die MeBergebnisse sollen hier den Berechnungsergeb- 
nissen der zitierten Arbeit gegeniibergestellt werden. Au®erdem soll durch den Versuch gezeigt 
werden, daB die der Berechnung zugrunde liegenden spezialisierenden und vereinfachenden Vor- 
aussetzungen (unendlich kleine austouliiny: icseschidiriicit und Undehnbarkeit der Membran), 
die das gesamte Problem einer Bes rilesronte Ge achemaretien Lésung zuginglich machen, sowie 
homogener Spannungszustand, durchaus mit guter Annaherung in der praktischen Ausfiihrung 
einzuhalten sind. 

Die Versuchsdurchfiihrung geschah in enger Anlehnung an die zitierte Arbeit. Angabe von 
Gleichungs- bzw. Tabellennummern bezichen sich auf eben diese. 


2. Versuchsanordnung. Eine Versuchsanordnung mit waagerecht angeordneter Membran mit 
Gewichtsspannung iiber Rollen und Sandschiittung hat sich aus verschiedenen Griinden als un- 
zweckmaBbig erwiesen. Fiir die Versuchsdurchfiihrung wurde daher die in Abb. | skizzierte An- 
ordnung gewahlt. Hierbei wurde durch einen 
aufpumpbaren Fahrradschlauch, der von einem 
Ringkasten und der Membran umschlossen ist, die 
Membran gespannt. Die auftretenden Reibungs- 
krafte lassen sich durch wiederholtes Uber- und - 
Entlasten weitgehend homogenisieren, so da ihr 
StéreinfluB auf den Spannungszustand als gering 
zu betrachten ist. Diese Anordnung hat neben der 
einfachen Spannungsregulierung den Vorteil, daQB 
der Membranrand gegen das Gehause abgedich- 
tet ist, so daB die Auswélbung der Membran leicht 
durch Fliissigkeitsdruck erzeugt werden kann. 

Das MeBsystem besteht aus einer AnreiBplatte 
als Grundplatte, auf der ein ParallelmefSistander 
waagerecht an einer mit MaBeinteilung versehenen 
Anschlagleiste verschoben werden kann. Der 
MeBstander tragt an seiner lotrecht verschieb- 
lichen Muffe eine MeBuhr mit einer MeBgenauig- 
keit bis zu 0,01 mm, deren federbelasteter Mef- 
stift zur Einregelung eines beliebig kleinen MeB- 
druckes durch eine einstellbare Gegenfederung 
gehalten wird. Abb. 1. Schemaskizze der Versuchsanordnung. 

MeBfebler bis zu + 0,1 mm sind wegen der 
gegebenen Bearbeitungsungenauigkeit, der unvollkommenen Formbestandigkeit und Empfindlich- 
keit der Anordnung zu erwarten. 

Die Auswahl des geeigneten Membranwerkstoffes wurde durch Vorversuche hinsichtlich Festig- 
keit, FlicBverhalten und Elastizitat aus einer Reihe von Gummi- und Kunststoff-Folien getroffen, 
Gewahlt wurde die von Kalle u. Co. hergestellte ,, Hostaphan“-Folie von 0,05 mm Dicke mit einer 
ZerreiBfestigkeit von mehr als 5 kp/cm (1000 kp/em?), die bis 3,5 kp/em (700 kp/em?) Spannung ein 
nahezu ee elastisches Verhalten bei einem mittleren [- Modul von 160 kp/em (32 000 kp/em?) 


ee 
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3. Versuchsdurchfiihrung. a) Membranspannung S und Kenngré8e a. Zur Bestimmung 
der Membranspannung wurde vor und nach den eigentlichen Messungen die maximale Auswélbung 


wax der Vollmembran unter achsensymmetrischer Gleichlast, erzeugt durch Luftdruck, ausge- 


messen. 


Tabelle 1. Die gemessenen Werte weg und die daraus bestimmten Werte ag und 
ihre Abweichung 0 in % vom theoretischen Wert « (Sollwert) 


| nf] é[4] te g [mm] ane 51%] 
0557 | 0 6.52000 0,98 =) 
05 +0,289 6,30 0,977 —2,3 
: —0,289 6,42 0,995 —0,5 
ae +0,464 6,34 0,983 —1,7 
: —0,404 6,45 1,000 0 
Le +0,568 6,30 0,977 —2,3 
; —0,568 6,40 0,993 —0,7 
sie 0,636 6,33 0,982 —1,8 
, —0,636 6,34 0,983 —1,7 
a +0,681 | 6,33 0,982 —1,8 
; —0,681 6,30 0,977 —2,3 
a 0 +0,707 6,30 0,977 —2,3 
3 —0,707 6,32 0,980 —2,0 
= Soa +-0,718 6,30 0,977 —2,3 
ef : —0,718 6,28 0,974 —2,6 
= ae +0,711 6,30 0,977 —2,3 
; —0,711 6,30 0,977 —2,3 
sone +0,090 6,34 0,983 —1,7 
: —0,690 6,30 0,977 —2,3 
Pee +0,650 6,40 0,993 —0,7 
j —0,050 6,35 0,979 —2,1 
ie +0,592 6,43 0,998 —0,2 
; —0,592 6,35 0,979 —2,1 
40,6 40,505 6,40 0,993 —0,7 
: —0,505 6,30 0,977 —2,3 
Sie +0,373 6,38 0,990 —1,0 
—0,373 6,23 0,966 —3,4 
Mie 0,053 6,25 0,969 —3,1 
3 —0,053 6,23 0,966 | —3,4 
+-0,803 0 6,23 0,966 | —3,4 
0 0 12,78 1,98 —1,0 
> Ae +0,201 127 1,97 —1,5 
3 : —0,201 | 12,8 | 1,985) | 0,7 
5 lah +-0,2275 Re Bey —0,15 
eS : —0,2275 | 12,9 | 2,00 0 
ES +0,201 12,73 1,974 —1,3 
J —0,201 12,9 2,00 0 
U cay +0,0816 12,8 1,985 =97 
; —0,0816 | 12,92 | 2,004 +0,02 
+0,414 0 12,9 a a'| dee 2200 0 


Weg max = 13,7mm}; &gmax = 2,12 gegeniiber &%max (berechnet) von 2,1126; 6 = 0,7% 


Mit den Gleichungen der Karasschen Arbeit 


ee ee 
Wmax = a ‘Se Va» (6 a) 
Zo eos 
eed 0 fos (23¢) 


kann daraus die Membranspannung und das Verhiltnis der unter hydrostatischer Belastung ge- 
messenen Durchsenkung w, zu der dimensionslosen GréBe ~ mittelbar angegeben werden. 
Gemessen wurde zu Beginn der MefGreihe w,,,, = 12,8 mm und zu SchluB der MeGreihe 
Wee = 13,0'mm. 
Nach (6a) errechnet sich die Membranspannung, wenn der Mittelwert w/,, = 12,9 mm und 
die Daten der Versuchsanordnung r, = 20,5 em, p’ = 20,5- 10-* kp/em? eingesetzt werden, zu 
"yX 90,53. 10-3 
= qr = aoe 1,66 kp/em 


max 
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und nach (23c¢) im Verein mit (6a) unter Beachtung, daB p’ = y r, ist, 


8S 2 Ue 2 ; 
Ue 8 es 5 ip 0, Loo 10,7, 
fe ro Winax oaks 


wenn w, in mm eingesetzt wird. 


b) Ausmessung der Wélbflache. Die Ausmessung der durch hydrostatischen Druck er- 
zeugten Wolbflache geschieht in enger Anlehnung an Berechnungsbeispiel 1 der Karasschen Arbeit. 
Es werden die Auslenkungen w, der in Tabelle 1 angegebenen Koordinaten fiir « = 1, x = 2 und 
Omax ausgemessen. Die MefSwerte sind durch den Index g gekennzeichnet. 


4, Beurteilung der Messungen. Die in Tabelle 1 zusammengestellten Werte a, sind durchweg 
kleiner als der jeweils erwartete Sollwert. Die mittlere Abweichung 6, liegt 


fiir den Sollwert « = 1 bei 6,, =—1,7 % (2 —0,11 mm bei w,), 


fiir den Sollwert « = 2 bei 6,, = —0,53% (4. —0,07 mm bei w,). 


Wenn auch die Abweichungen noch an bzw. in den Grenzen der eingangs genannten und spater 
durch die Unsymmetrie der MeBfehler bestatigten MeBzuverlassigkeit liegen, bleibt ihre Tendenz 
jedoch offenbar und fordert eine Erklarung, die mit dem Hinweis auf zwei Stichworte gegeben sei: 

1. Zeitabhangigkeit des Spannungs-Dehnungsverhaltens, ersichtlich aus den MeBwerten unter 
Ziff. 3a), wobei der zeitliche Mittelwert fiir w/a, in Ubereinstimmung mit dem Versuchsablauf 
naher an dem ersten MeBwert von way = 12,8 mm lag. 

2. Dehnungsvermégen der Membran unter Beachtung des Reibungsschlusses am Membranrand, 
wodurch bei Belastung durch hydrostatischen Druck der Einflu8 inhomogener Stérspannungen 
kleinere als durch Berechnung gegebene Auslenkungen erwarten laBt. 

Eine Ermittlung von « nach (6a) mit Wyay = 12,8 mm statt 12,9 mm, was nach Ziff. 4.1 und 
4.2 naheliegt, hatte fiir « = 1 etwa eine mittlere Abweichung 6, = —0,8% und fir « = 2, 
6, = +0,3% ergeben. 

Unter Beriicksichtigung der durch den Versuchsaufbau und die Materialeigenschaften der Mem- 
bran gegebenen Fehlerquellen darf das Versuchsergebnis als exakte Bestatigung der Karasschen 
Berechnung der Volikreismembran gelten. 


(Eingegangen am 11. Mai 1957) 


Anschrift des Verfassers: Dipl.-Ing. Karl-Ernst Meier-Dérnberg, Darmstadt, Technische Hochschule, Institut 
fiir Angewandte Mechanik und Technische Schwingungslehre 


96 S. Falk: Die Berechnung geschlossener Rahmentragwerke usw. Ingenieur-Archiv 


Die Berechnung geschlossener Rahmentragwerke 
nach dem Reduktionsverfahren 


Von S. Falk 


1. Einleitung, Ziel der Arbeit. Die vorliegende Arbeit behandelt als Fortsetzung meiner beiden 
fritheren Aufsiatze1, deren Formeln und Abbildungen mit I bzw. IT zitiert werden, geschlossene 
ebene und raumliche Rahmentragwerke. Zwei Beispiele werden ausfihrlich vorgerechnet. Den 
Schluf bildet eine kurz gefaBte vollstandige Theorie des Reduktionsverfahrens fiir beliebige Systeme 
von linearen gewohnlichen Differentialgleichungen (Schwingen, Knicken, elastische Bettung usw.). 


2. Zwei- und mehrfache Koppelfedern. Bislang haben wir nur die einseitige Koppelfeder nach 
Art der Abb. 1 kennengelernt. Um auch die Matrix der beidseitigen (zweifachen) Feder nach Abb. 2 
zu gewinnen, denken wir uns einen durchlaufenden Trager oder ein offenes Tragwerk an zwei be- 
liebigen Stellen M und N durchgeschnitten (Abb. 3), rechnen, in M beginnend, mit den dort vorhan- 
denen 20 Kraft- und DeformationsgréBen als FreigréBen bis zum Punkte N und finden abnlich 
wie unter II (37) 


(dy dy. +. diay (hg kya +++ ha) I 
( * x x * d, dy 
ge ey, aS Oa x oe x | de dy 
cecal earl 
* * St ee ee * * oT peek. oO d, 
ae ees Se (1) 
Bey ps OF el ceo. Oe . sae 
knoe races EAT ctr bes keg 
ee ene ok x ok x hy ky 
0 0 0 0 0 Qed l}/wn 
oder kurz 
©, Dy + Gs. by = (0 —d)y (2) 
Ss Dy + Gly = (€— f) 
und daraus 
— ly = Oy WO —Dd)y + ©,* © dy » 
dT, 6,606 by te, Se phest 
Die Gesamtmatrix 
Cun = e 2 Ss a = (4) 
C163" —"© 1G Cae, 


ist nach dem Satz von Maxwell- Betti nebensymmetrisch und heiBt die Federmatrix der Rahmen- 
schnitte M, N. Im wichtigen Sonderfall des einfeldrigen geraden Balkens mit konstanter Biege- 
steifigkeit lautet sie zum Beisicl fiir die beiden Paare Gc M y) und (w, Q,) 


— 6F —3/1 6/2 —ayl 
| 12/8 6 ==12/8 6p = w\  (M 
ents eet era ey 6ien=2 (| Pee -("). t=(0"). (5) 


—12/8 —6/P 12/8 —6/P 


* S. Falk, Ing.-Arch. 24 (1956) S. 216; S. Falk, Ing.-Arch. 26 (1958) S. 61. 
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Sind insbesondere die Endpunkte unverschieblich, Wy = Wy = 0, so folgt hieraus nach Streichen 
der ersten und dritten Spalte, wenn wir noch die Zeiger M und N durch 0 und 1 ersetzen, 


— M, AS i ae a 1-1 
GEA 3) Po ): Y 


Wenn die Belastung fehlt, entnehmen wir daraus die Gleichung der geraden Momentenlinie (Abb. 4) 


ae x EJ as 
M(x) = M, + (M,—M,)-2 =* (49 +291) — 9 + 69)-7]. (7) 
4 n 
K K 
———P <¢-———-- 
| 
2-4 
wo 
K=czr 
Abb. 1. Einseitige Feder. Abb. 4, Linearer Momentenverlauf iiber einem unbelasteten geraden 
Balken. 
My 
K K 
Se NNN | enn 
_——— a 
1 
L, Cc 
BthAADAAAW) 
K=0 (24-2; M, 
(%1-Zz) z aye 
Abb. 2. Beidseitige Feder. Abb. 5. Drittelpunkt bei fester Einspannung und fehlender Belastung. 


&& 
Abb, 3. Ersatz eines Rahmenteiles M N durch cine beidseitige Abb. 6. Riumliches Tragwerk mit Schnittflache; 
(belastete) Koppelfeder. schematisch. 

Das Biegemoment M(x) verschwindet im Punkt 

a 2htH @) 
3 Yy +39 
Fiir py = 0 folgt daraus «* = 1/3 (Drittelpunkt bei fester Einspannung, Abb. 5) und weiter aus (7) 
2a be 

M(x) = ‘ae if 3 1: (9) 


Wir haben in (3) gleich den Vektor —f,, (nicht + fz) auf dic linke Seite gesetzt; das aus —fy,, fy 
und re bestehende Kraftesystem ist dann im Gleichgewicht. Oder anders ausgedriickt: die Glei- 
chung (3) gibt an, wie eine im Bereich M N gegebene Belastung bei vorgeschricbenen Deforma- 
tionen Dy und Dy auf die Umgebung ,.abgesetat“, d. h. in welchem Verhaltnis sie verteilt wird. 
Das zugehiérige Symbol zeigt Abb. 3 unten. 

Ganz allgemein lassen sich nun durch p Schnitte insgesamt 7 = 0; + 02 +°** + Qp KraftgréBen 
des Tragwerkes bloBlegen, die mit ihren konjugierten DeformationsgréBen iiber eine z-reihige qua- 
dratische Federmatrix zusammenhangen (p-fache Koppelfeder); die Reziproke dazu heist eine 


p-fache EinfluBmatrix. 


3. Der Begriff der gemeinsamen Reduktion. Eine Schnittflache moge ein beliebiges Rahmen- 
tragwerk in zwei getrennte Teile zerlegen und dabei in s Schnittpunkten treffen (Abb. 6). Ahnlich 


u 
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wie frither I. Ziff. 2 kénnen wir jetzt die Gesamtheit der zweimal 
Nie oe (10) 
1 


Kraft- und DeformationsgréfBen in der Schnittflache ¢ reduzieren, und auch hier lauft die Aufgabe 
darauf hinaus, die zunachst unbekannten Kraft- und DeformationsgréBen an der ,,Anfangsflache* 
e, zu finden. Sind diese bekannt, so laBt sich das Tragwerk in einem einzigen Rechnungsgang bis 
zur ,,Endflache“ ¢€, durchrechnen. 


4. Die geordnete Strangschar. Wir betrachten jetzt s parallele gleich lange gerade Durchlauf- 
trager — wir wollen sie in diesem Zusammenhang ,,Strange“ nennen — als Kanten eines geraden 
Prismas iiber einer Grundflache ¢) und legen n dazu parallele Schnitte hindurch (Abb. 7). Koppel- 
federn, die in Abb. 7 der Einfachheit halber durch gerade 
Verbindungslinien dargestellt sind, sollen allein in diesen 
Schnittebenen liegen, und zwar sprechen wir von ,,Nach- 
barkopplung“*, wenn héchstens je zwei benachbarte Ecken 
durch Koppelfedern miteinander verbunden sind (wie in 
der Ebene ¢,, Abb. 7), sonst von ,,durchgehender“ Kopp- 
lung (so in den Ebenen ¢, und ¢,,). Nun fassen wir die Kraft- 
und DeformationsgréBen samtlicher s Strange zu einem 
einzigen Vektor ¢ zusammen 


f= (0, per De ee (11) 


a Sa ees oem ers ee 
Vie AN faa ao 


Abb. 7, Geordnete Strangschar; schematisch. Abb. 8. Einige ebene geordnete Strangscharen. 


und bauen die zugehérigen Leitmatrizen aus den Matrizen II (10) bis II (13) auf. Haben wir zum 
Beispiel s = 2 Strange, an deren jedem allein das Paar (p, My) beteiligt sei, so liefert II (18) im 
Falle unverschieblicher Knoten (w; = 0) 


Gee 0 0; S208 Jase | a0 Ya 
| 0 —2 4L2ES, 0 % | 0 Po 

OeOE Sane? 0 seve ceastiag, (12) 
lorry, 0 Oey 22 a ee | M, 

0 0 0 0 foreman | 1 


Von den Komponenten des Vektors y (11) verschwinden an der Anfangsflache r Komponenten, 
die tibrigen r FreigréBen werden durch die r Randbedingungen an der Endflache festgelegt, und 
mit dem nun bekannten Anfangsvektor Y, wird die ganze Rechnung wiederholt. Diese Rechnung 
verlauft somit auch hier nach dem Schema IT (34) und wird im Prinzip nicht gestért, wenn wir das 
Prisma der Abb. 7 noch auf folgende Weise verallgemeinern: 


1, Auf jedem Strang werden noch beliebig viele Feldgrenzen eingefiigt (Anfederungen, Ablenk- 
winkel, Zwischenbedingungen). 

2. Die Punkte M_,,, werden auf ihren Strangen verschoben, doch nur so, daB ihre Reihenfolge 
ungedndert bleibt; allenfalls diirfen einige Punkte M,,, zusammenfallen. 

3. Die Strange brauchen weder geradlinig noch parallel zu sein, diirfen sich nur nicht schneiden. 

Wir gewinnen hiernach eine groBe Mannigfaltigkeit von raumlichen Netzwerken, die wir alle als 


geordnete Strangscharen bezeichnen wollen. Kinige einfache ebene Netze zeigt die Abb. 8, wo die 
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Strange durchgezeichnet sind. Zwar ist die Strangzahl s keine fest definierte GréBe, doch fallt 
bei den iiblichen Bauwerken der Praxis die Auswahl der kleinstméglichen Zahlen s bzw. r nicht 
schwer. 


5. Die Ausnahme-Strangschar. Nicht jedes Tragwerk ist einer geordneten Strangschar topo- 
logisch aquivalent. Es kann durchaus vorkommen, daB sich einige oder alle Strange schneiden. 
An solchen Schnittpunkten — in der Abb. 9 z. B. den Punkten A, B, C, D, E — sind im allgemeinen 
Zwischenbedingungen zu erfiillen, so etwa in A: 9, = Y, = Fy = vy 3) My = 0. Dab sich jedes 
 tiberhaupt denkbare Rahmentragwerk auf eine Strangschar zuriickfiihren laBt, liegt auf der Hand: 
man braucht ja nur jedes einzelne Feld zum 
Strang zu erklaren. 


aN 
\ XN 


. Abb. 10. Langsstarres Tragwerk mit gemeinsamen 
Abb. 9. Zwei Ausnahme-Strangscharen. Knotenverschiebungen w, 


6. Tragwerke mit einheitlichen Zwischenbedingungen. Bei Tragwerken mit einheitlichen Zwi- 
schenbedingungen kann es oft von Vorteil sein, den Vektor x (11) anders einzufiihren. Das prak- 
tisch wichtigste Beispiel zeigt die Abb. 10: infolge der Langsstarrheit der senkrechten Zwischen- 
stiitzen sind die Durchbiegungen w an der Stelle a (ebenso in b, c usw.) einander gleich; zweckmaBig 
ist daher die Einfiithrung eines konjugierten Paares (w, Q,.;) mit der gemeinsamen Durchbiegung 
w= w, =--: =w, und der Querkraftsumme am Knoten Q,., = Q,; +--:+Q,. Fir ein recht- 
eckiges ebenes Tragwerk nach Abb. 10 mit zwei (und ahnlich mit s) Strangen lautet dann die zu-. 
gehérige Leitmatrix 


I = ee hie 3 lm 0 00) 
3/1, ase) 0 0 (2EJ), 0 Po pees TREO SO | 
3/1, 0 SS (ek 1); 0 0 @» = Si (080 
—6(EJ/P), 0 6(EJ/D, ——2 0 0 M, |6(EJ/P), * *| (13) 
|—6(EJ/P), 6(EJ/), 0 0 = 0 M, |6(EJ/P),* * 
0 0 0 0 0 1 On cay ma OL 
0 0 0 0 0 0 1 0 00 
(w Pa Lr M, he Ue 1) = ¢. 
Dabei wurde zur Abkiirzung gesetzt 
ETI m= 2m, REA sia hk 6 en (14) 


Der mittlere Teil des Kerns der Matrix (13) enthalt die KenngréBen der Paare (p, M) und wird um- 
rahmt von den (zum Teil durch (14) definierten) Kenngréfen des allen s Strangen gemeinsamen 
Paares (w, Q,,,)- Der Ausdruck )’ (w 1)/m stellt gewissermaBen das gewogene Mittel der durch die 
Belastungen hervorgerufenen Durchbiegungen dar. 


7. Einige spezielle Tragwerke der Baustatik. Wir stellen jetzt einige in der Praxis immer wieder- 
kehrende Typen von Tragwerken zusammen und geben ihre kennzeichnenden Eigenschaften an. Es 


bedeuten: 


7 6 
‘ 
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Anzahl der Strange ‘Sie 
keine hee 
Art der Kopplung 4 Nachbarkopplung _ ko 
durchgehende Kopplung k, , (15) 
Anzahl der beteiligten Paare je Strang (as 
Art der Leitmatrix. (Hinweis auf IT. . .) (lias 


Ordnung r des zu lésenden Gleichungssystems r= 2} Q, . 
1 
Damit ist Tab. 1 wohl ohne weitere Erlauterungen verstandlich. 


Tabelle 1. Einige spezielle Tragwerke der Baustatik 


1. Durchlauftrager auf festen Stiitzen. 


gss lia, Oral (INS EB) fa AI, TD? 

2. Gerader Durchlauftrager auf unverdreh- = — sonar hee =e Y 
baren Stutzen. = iy 
Gas ore Or) ve — dle LY: TE: 

3. Gerader Durchlauftrager, elastisch dreh- acini Y 
und senkbar gestutzt. : Z Y 
Simp Ol 2a ly 9) ) ane Pr > 


4. Fachwerkbriicke, Knoten eckensteif. s = 2, ky. 
Falls langsstarr (Knoten unyerschieblich): 
OE (LS) Sree 
falls Biegung mit Langsdehnung: @ = 3, (II. 14), r= 6. 


5. Stockwerkrahmen, Knoten unverschieblich (langsstarr). 
S$, pes (Oh ils (iI. 18), | ean NY, 


6. Stockwerkrahmen, Knoten 
verschieblich (langsstarr). s, kp. 
Entweder 0 = 1, (III. 13), r=s. 


7. Tragerroste (Kreuzwerke). 
Biegen allems (GT = 0)is, ka.10 = 2-0 (LO) eres: 
Biegen und Drillen: Sits = 85 (Goagh 1? Sah 


8. Die praktische Durchfiihrung des Verfahrens. Wir stellen jetzt ein allgemeines Programm 
auf, nach welchem sich jedes ebene oder réumliche Rahmentragwerk (auch Fachwerk) berechnen 
1aBt. 

1, Die Charakteristik des Tragwerkes wird nach (15) festgelegt. Man entwirft eine Skizze der 
einzelnen Strange und stellt die Hauptteile der Leitmatrizen nach II (10) usw. zusammen. Falls der 
Rahmen abgewinkelt ist, miissen diese Hauptteile noch mit den Transformationsmatrizen & I (35) 
multipliziert werden. 

2. Kommen im Strangschema Koppelfedern vor, so entnimmt man deren Federmatrizen einer 
Tabelle oder berechnet sie nach dem in II. Ziff. 6 baw. ILI. Ziff. 2 geschilderten Verfahren. Diese 
stets nebensymmetrischen Federmatrizen werden unten rechts in die Leitmatrix CII (15) einge- 
tragen. 

3. Alle Werte liegen jetzt formelmaBig vor, z. B. 6 E J/P, 21 cos 40° und dergleichen. Diese 


Werte werden dimensionslos gemacht mit Hilfe zweier zweckmabig zu wihlenden VergleichsgréBen 
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P* (Kraft) und 1* (Lange), indem man simtliche vorkommenden GroBen als Vielfaches von Ld he 
P* I*, P*/I* usw. darstellt, die Buchstaben P* und I* aber wahrend der Rechnung fortlaBt (oder 
in etwas nachlassiger, aber gebrauchlicher Redeweise: indem man sie gleich Eins setzt). Bei homo- 
hi ee (und wohl auch sonst) empfehlen sich /* = 1 und P* = 6 ES baw. B= 

4. Nun multipliziert man die r + 1 Spalten unabhangig voneinander an siimtlichen Leitmatrizen 
von oben nach unten vorbei. Nach jeder Matrizenmultiplikation ist die Sumimenprobe zu machen, 
bei groBeren Matrizen schaltet man Zwischenproben ein. Unten rechts neben der letzten Leit- 
matrix UE findet man dann die r NullgréBen vor. (Die wbrigen r GréBen g,; brauchen oft nicht be- 
rechnet zu werden.) 


5. Man lost das Gleichungssystem IT (33) 
et = by aa) (16) 


nach dem Gaufschen Algorithmus oder durch Iteration und findet mit q zugleich den Anfangs- 
vektor Yo. Die Lésung ist aus mechanischen Griinden stets vorhanden und eindeutig. 

6. Mit dem nun zahlenmaBig vorliegenden Anfangsvektor y, wird die Rechnung in einer einzigen 
Spalte, der ,,Endspalte“ wiederholt. Die r Bedingungen rechts miissen von selbst erfiillt sein. 

7. Jetzt sind auch alle etwa abgeschnittenen Zweige (ein -oder mehrfache Koppelfedern) nach dem 
gleichen Schmea in einer einzigen Spalte zu berechnen, da deren Anfangswerte nach II. Ziff.6 bzw. 
III. Ziff. 2 ebenfalls festliegen. 

8. Damit sind die Werte der an der Deformation beteiligten Kraft- und DeformationsgréBen fiir 
alle Feldgrenzen, also nach I (2) usw. fiir den ganzen Rahmen gegeben. Auch etwa unterdriickte 
Paare lassen sich leicht wiedergewinnen IT (20). 

9. Zum SchluB priift man nach, ob die gegebene Belastung und die berechneten Reaktionen zu- 
sammen im Gleichgewicht sind. 

Soll dasselbe Tragwerk fiir verschiedene Belastungen durchgerechnet werden, so sind auszu- 
wechseln: die Lastspalten in den Leitmatrizen, die inhomogene Spalte im Schema II (34), der 
Vektor r in (16) und die Endspalte. Es macht auch keine grundsatzlichen Schwierigkeiten, einige 
Kennwerte E J, 1 usw. als frei verfiighare Parameter mitzunehmen, auch lassen sich nachtraglich 
neue Feldgrenzen einfiihren I (13). 

Ein Umstand ist indessen an dieser Stelle erwahnenswert: Aus mechanischen Griinden liegen 
die Werte der Endspalte bei verniinftiger Wahl von P* und /* etwa zwischen 0,1 und 10, auf alle 
Falle aber in derselben GréSenordnung. Nicht so die Werte der homogenen Spalten: diese schwellen 
in der Regel mit wachsender Felderzahl stark an, was sich oft schon beim vierten, fiinften Feld 
unliebsam bemerkbar macht. Mit anderen Worten: die in der GréSenordnung gleichbleibenden 
Endwerte sind die so gefiirchteten .,,kleinen Differenzen groBer Zahlen“ und das um so mehr, je 
hoher die Nummer des Feldes ist. Bei der Wiederholung mit dem Vektor 7) wird daher im allge- 
meinen eine Gleichung 


Noe pase 0 (17) 
mit einem Fehlervektor j entstehen. Jetzt lésen wir das Gleichungssystem 
: Aa +7 =0 (18) 


und wiederholen die Rechnung mit dem neuen Vektor xj”, diesmal aber mit der letzten Ziffer 0 
statt 1, wodurch die gegebene Belastung automatisch ausgeschaltet wird, und bekommen als Summe 
zweier Endspalten das Ergebnis 

Git hy. (19) 
wobei wir eine gewisse Anzahl neuer geltender Stellen gewonnen haben. Falls notwendig, ist dieser 
Vorgang hinreichend oft zu wiederholen. In fast allen praktischen Fallen jedoch geniigt es, den 
Vektor a aus (18) zu ermitteln und mit seiner Hilfe an Hand der homogenen Spalten die Anzahl 
der geltenden Stellen abzuschatzen, was meist mit einem fliichtigen Blick auf die Rechnung erledigt 
ist. 

Das Anwachsen der homogenen Spalten kann man aber auch grundsatzlich verbindern, indem 
man die Trager zerschneidet und mit Koppelfedern rechnet; allerdings miissen dann mehrere Glei- 
chungssysteme gelést werden. Aber auch durch geschickte Wahl des Anfangsvektors laBt sich das 
Anschwellen weitgehend verhindern, indem man nicht mit den r FreigréBen selbst, sondern ge- 
eigneten Linearkombinationen losrechnet, etwa so, daB die ,,rechte Seite y in (16) klein bleibt 
oder die Hauptdiagonalelemente der Systemmatrix 9| stark iberwiegen (zum Zwecke der Iteration) ; 
doch sind das Einzelheiten, die hier zu weit fiihren. 


102 S. Falk: Die Berechnung geschlossener Rahmentragwerke usw. Ingenieur-Archiy 


Sind mehrere oder alle Leitmatrizen einander véllig oder doch annahernd gleich, so schreibt 
man sie ein einziges Mal auf ein bewegliches Blatt, das an den Kolonnen des Schemas II (34) von 
oben nach unten vorbeigezogen wird. Bei groBer Strangzahl s werden die Matrizen ziemlich um- 
fangreich — sie enthalten dann relativ immer mehr Nullen — und daher leicht uniibersichtlich ; 
man wird in solchen Fallen, namentlich bei nur schwacher Kopplung (insbesonder Nachbarkopp- 
lung) die Strange einzeln baw. gruppenweise nebeneinander herfiihren und die unterdriickte Kopp- 
lung nach jeder Matrizenmultiplikation gesondert nachholen. 


9. Die Programmierung fiir digitale Rechenautomaten. Das Reduktionsverfahren verlangt in 
der hier entwickelten Form lediglich zwei Operationen: das Multiplizieren und das Umkehren von 
Matrizen, beides aber gehért zum Grundprogramm auch der kleinsten Ziffernmaschine. Will man 
die in IT (22) definierte verschrankte Multiplikation umgehen — obwohl auch deren Programmie- 
rung keine besonderen Umstinde erfordert —, so rechnet man nach II (23). Um den allgemeinen 
Rechnungsgang auf einem Automaten zu zeigen, entwickeln wir jetzt das Programm fiir eine ge- 
ordnete Strangschar (Durchlauftrager, Stockwerkrahmen, Tragerrost und dergleichen), und zwar 
werde das Tragwerk gleich fiir b verschiedene Lastfalle berechnet. Das zugehérige Matrizenschema 
sieht so aus: 

< 


iff 


wy eS 
2 
| 
| 
Ss) 
| 
Ep TN RR REE a are ———————————————— 


NUE | OF 
+ — — 
‘. 
b yD) & §) y) & 
(| inet lice ee } 
 Leitmatrix 
& Kern 
© Federmatrix 
(© Einheitsmatrix 
§) Nullmatrix v v 
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Und die Anweisung fiir den Rechenautomaten lautet: 
1. Eingabe: n Matrizen Q; (Leitmatrizen), 
Matrix % (homogene Spalten), 
Matrix % (6 Belastungen). 
. Berechne und speichere das Produkt Bae es, a ie 
. Berechne % = J &B. 
. Greife aus ¥8 r bestimmte Zeilen heraus und fasse sie zur Matrix 9{ zusammen. 
Bilde und speichere Q{-1. 
. Berechne das Produkt eee SES 
. Greife aus 8 r bestimmte Zeilen heraus (die gleichen wie unter 4) und fasse sie zur Matrix 3) 
zusammen. 
8. Berechne D = —Y-1 MM. 
9. Berechne ) = 6+ 89. 


NAVE wp 


Niue ruelgce A‘elle .q) Keusel” 4, care 


7 (iss oe gy 9) 
Die unter 10. ausgedruckten Matrizen enthalten die Kraft- und Deformationsgréfen an simt- 
lichen Feldgrenzen fiir die b vorgeschriebenen Lastfalle in den Dimensionen P*, /*, P*/* usw. 


10. Berechne und drucke aus I; al §) ’ 


Z1,Xp 


Abb. 11. Fiinfstéckiger Rahmen mit unverschieblichen Knoten, dazu Strangschema. 


10. Beispiele. Wir fihren nun zwei Beispiele vor und richten uns dabei nach dem Programm 
aus Ziff. 8. 

a) Erstes Beispiel. Fir den fiinfstéckigen Rahmen von Abb.11 ist der Biegemomenten- 
verlauf infolge der angegebenen Windbelastung zu berechnen. Die Knoten sind unverschieblich. 

1) s=2,k,, o=1, II (18), r=2. Strangschema Abb.11 rechts. Samtliche Koppelfedern sind 
unbelastet. 

2) Es ist C7. =4-2E J/L=8E J/l nach Tabelle II, 2; die Feder C, des zweifeldrigen Rahmen- 
winkels unten rechts ergibt eine kurze Rechnung mit der Matrix I (18) zu C,=7EJ/l. Die Matrix 
der Koppelfeder ist nach (6) 


4 eed, : ote rere. Y, 
M, a es ee Y 
3) Wir wihlen P* =2EJ/P, I* =) get2P*l*, oh, EJ a2, lai, q= 12, 


4) Zahlenschema 1. Die Leitmatrix finden wir unter (12) vor; die Federmatrix &,, unten 
rechts in & ist hier 


0 —C, (baw. 0) —1 —2\ /—1 —5,5 ae 
ae 0 \+(=, = er sae 
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5) Die beiden Randbedingungen rechts (d.h. hier oben) lauten 
M,, = — 18,600 My, — 11,7145 My, — 0,155 = i My, = 9,003 83667 , 
i ro, 
M;,; = — 82,394 My, — 24,795 M,, — 0,163 = 0 M,, = — 0,0193232. 
6) Siehe Zahlenschema 1. In den homogenen Spalten ist mehrmals nach abgeschlossener Ma- 
trizenmultiplikation der Faktor 10 herausgezogen, um ein unndétiges Anschwellen der Zahlen zu 


vermeiden; auf die eingerahmten Endgleichungen ist dieser Rechenkniff ohne EinfluB. Die Werte 
der Endspalte sind einer achtstelligen Rechnung entnommen und aufgerundet. 


P*Y*=G1 2 
———— 


ie 
Abb. 12. Momentenverlauf fiir den Stockwerkrahmen von Abb. 11. 

p Pp 5k 
(fed I; (2 Pp 


SPETET CTT TT TT {| 
Vez Gd ah | ie y 
pros cesta agama 


WH, Z “, 
Z L Z Z Z Z 


Abb. 13. Symmetrisches Tragwerk mit gemeinsamen Knotenverschiebungen w, Abb. 14, Linke Hialfte des Tragwerkes von Abb. 13. 


das 


7) Fir alle eingespannten Rahmenteile gilt infolge fehlender Belastung die Gleichung (9): 
Drittelpunkte in Abb. 12. Fiir die waagerechten Mitteltrager gibt (6) den Momentenverlauf. 

8) Momentenverlauf Abb. 12: iiber den drei belasteten Feldern parabolisch, sonst linear, Dimen- 
sionen: M in P*I/* —2E J/l, q=12 P*/I* =24E J/l, also erscheint M in q (7/12. 

9) Ableiten der Momentenlinie gibt den Querkraftverlauf, womit alle 24 Reaktionen bekannt 
sind. Wie man sich leicht iiberzeugt, ist das gesamte Kraftesystem im Gleichgewicht. 

b) Zweites Beispiel. Fiir das symmetrische und symmetrisch belastete Tragwerk von Abb. 13 
sind Querkraft- und Momentenverlauf zu ermitteln. Infolge der starren, gelenkig angeschlossenen 
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Hoes ame Zwischenstiitzen sind die Durchbiegungen beider Tri ager an den Feldgrenzen einander 
gleich 
1) s=2, ky, 0, =o, =1, dazu das Paar (w, Ore (lo) yar 3. 
2) Reppelfedern sind mich vorhanden. 
3) Wir wahlen 1* =1, P* =6E J/I?, ferner P = P* und haben damit P = 1, I=1, 
EJ =1/6, Clg =O. 
Die in die Leitmatrix (13) einzusetzenden Werte sind nach (14) und den Tabellen I. 1 bzw. 


TW Meese 


Feld 1: 

ft, = E J/P = 1/6 1,1 5 fq =EJimP = 0,8; bad 2a leas 
Feelin iss 6 © 24 24) a, =2EJ/m (202 =0,4; 8, =1/2m(2)) =1,2, 
Seer Che es 9 3 1 Ay eo 
3 oF ae aaa} t)[ py =~ 125. 


Mo: M), 1 
YI 0 0 0 \ ee \ 
9, 0 0 0 
Zahlenschema 1 Fe | 0 1 0 | BS 
M, 1 0 0 0,0038 
1a 0 0 1 
a Ee ea 0 0 1 0 0 0,038) 
ue es ae | O — 0 0 1 0 — 0,0193 
Se |) Qicipc ia eed ig ey eed Us 0 | 0,1411 
Aoi =e ear eee ee —8 eal 0 — 0,0113 
OMIT ak Ose > One ch 0 0 0 0 1 a 
Be I) SY IO 6 = 0 2510 = 0 — 0,0190 
ee aie For G 0 0 =| Pos 0 0,1797 
| oe as ee 14 38 0 — 0,6806 
sie 6” MO 06 1 80 17.5 0 me 
oO” eo v0 0 0 1 
SO MiieretO ie 0s 0. = 6 ( 10 19.5 \ | a 
tg: eekeys eer Qeowiegy 60 1,6 5.1 a — 0,0400 
; ee ree i |-2 393.8 §=503 | — 1,0264 
5 (Tipe Pees a a er eT —80,6 21,2 me J\ 9081 
Th er ee 0 0 0 0 
(— (peacoat Tl"G i F106 eal) 001/006 
eee of 0 0 ea ee 0 0G 0,0537 
‘| hi ee (ee 25) 17,70 16,02 0,23 | |—1,1650 
Oe 6 81,43 23,405 Oe Bee 
(ewe 0e 0 A 0 0 0 0 
Pepe geo Te 3G 0 70 ( 10,155 28425 ~—«0,016 0,0358 
(ee eee 0 0 2,160 2,306 0,036 | | — 0,2724 
&). 0 3s —2 0-—8 1 —2 EBS es =e 600 —11,7145 == 0,155 | 0,0000 
ee ee On 61} | 289-394 24.705 103 e piad 
(ne at ea 080 0 0 
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Feld 2 

Ue — EJ/B == 1/6 {| 1 7 Xe = E J/mP => Orb: By = 1/2 ml => ODE 
m= == 

fy =3 EJ/B = 1/2 6 ' 2°" 3 | =3 EJ/ml = 0,75; 9p, = eon tO i 

ee 2 6-1 i 5 , 

Pa = 9; hi She RG 0,5, Qazes = 9; Q: ees = 6, 


H =05, (Sw p)/m=0 +5 0,375. 


4) Siehe Zahlenschema 2. 

5) Das im Zahlenschema eingerahm- 
te Gleichungssystem (Randbedingungen 
Abb. 14!) hat die Lésungen: 


Por = — 1,864154, 
M,.—= 1,801790, 
Min[Pl] Qoges = — 3,325922. 


6) Zahlenschema 2, Endspalte. 

7) Entfallt. 

8) Abb. 15 und 16. 

9) Das Gleichgewichtssystem der 


Bice Krafte und Momente zeigt Abb. 17. 
Abb, 15. Momentenverlauf fiir das Tragwerk von Abb. 13. 


11. Verallgemeinerungen und Ergin- 
zungen. 1) Die Steifigkeiten 
EJy, EJ,,GI und EF in II (3) 


sind nicht konstant, sondern beliebige 
(oft unstetige) Funktionen der Koordi- 
nate x. Davon werden offensichtlich 
allein die Steifigkeitsmatrizen % in IT (15) 


4,876 


4,320 UIE4 


Gilt) betroffen. Da die Integration der Diffe- 
2646 rentialgleichungen II (3) trivial ist, brau- 

nies chen wir nicht weiter darauf einzugehen; 

; beispielsweise ware in IT (10) die Ma- 

trix %, die dort wegen 90=1 aus dem 


einzigenElement —1/EF besteht, zu er- 


setzen durch 
! 


f a 12P P P 2 dé 
E F¢é) s 


Abb. 16. Querkraftverlauf fiir das Tragwerk von Abb. 13. 


(20) 
f 

6136P 2366P  —_-6,136P 2) Die vier Differentialgleichungen 

II (3) kénnen teilweise miteinander ge- 

1802 PL C ae koppelt sein (wenn z. B. der Schubmittel- 

punkt nicht mit dem Profilschwerpunkt 

2546P zusammenfallt). Dann lautet lediglich 

Abb. 17. Gleichgewichtssystem am Tragwerk von Abb. 13. das Fundamentalsystem anders; die 

Ordnung der Matrizen bleibt die gleiche. 

3) Die gewahlten Richtungen y und z in Abb. II (3) fallen nicht mit den Haupttragheitsrich- 

tungen des Balkenquerschnittes zusammen (schiefe Biegung). Dann sind die Belastungen g, und 
qy mit Hilfe des Ausdruckes 


0,680P 0,680P 


zeugp  2959P 3830P — 2759P 


Jy J,— Jy (21) 
in bekannter Weise zu kombinieren. 
4) Federn mit Vorspannung. Wenn die Anschlu®(koppel)federn in der unverformten (spannungs- 
losen) Lage des Tragers nicht selbst entspannt sind wie in Abb. I (5), so gilt an Stelle von II (4) 
allgemeiner 


k; = — ¢;(d; — d?) baw. f= —C(d—d*). * (22) 


2) Verteilen der Belastung auf die Feldgrenzen. Oft kann es zweckmabig sein, die in einem Feld 
angreifende Belastung auf beide Feldgrenzen zu verteilen. Zu diesem Zweck bringt man in i—1 
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24992 
— 2,8723 
—1,5355 
0,6405 
— 0,8444 
—4,0851 


CT ETE TE ET CORE ig ER ERs 
wd LO 
ee ele Giel 
Bese eee ee ia 
HSe oa oa qe Ge aol SY SS as 
3 CO HO HO ai is 
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ase | 
wm 1A 
s Se So sph thse Ll sek 
SS i a a Lc 
Ww [LA 
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Oe CS SSO Sn Ss Sto 2 lt. 
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ee NS ee a Te en ee 
Il 5p Se eo ee S&S 2 S 
EE EE AE 
Seo a qe es SG Ss S 
yo 7s 
Sf Sas 19 19 QD 
Qo a SS) SS SS er SC eo) NN 
3 ae ae 
te) In 10 We) ie) 
Ve) ANN Wel al Sse r= 
ork St on oD ww ae) 
So AS ooa a So Se Se SSS 
2, 0 A a 
So So otir Se oe eae 
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NAMOoOONSOCN OMOOSONOH 
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F x Ate a 
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os 
e | | | | 
10 wD 
60. ee “. 
a) SS Sh Se SN SS ee Se! 
Wey Ye 
eso ra) 60 A) 8) rt SD 


fe dee as ee 


Se ee SS eS eee 


und i konstante Krafte und Momente an und bestimmt sie nach Tabelle I, 1 baw. I, 1 so, daB 
sie der gegebenen Belastung im Felde i aquivalent sind. 


12. Der Vergleich mit der Eliminationsmethode. Wir stellen zunachst die Grundziige der beiden 
Verfahren in kurzen Worten einander gegeniiber: 
Die Eliminationsmethode | Das Reduktionsverfahren 


eliminiert eine oder mehrere Kraft- bzw. Defor- | bela&t samtliche an der Deformation beteiligten 
mationsgréBen aus | Kraft und DeformationsgréBen in 


ihrem natiirlichen Zusammenhang 
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Die Elimination gelingt nur durch Ubergriffe | Die Grundgleichungen werden fiir jedes einzelne 
von einem Feld auf alle benachbarten (oft auch | Feld getrennt und unabhangig voneinander auf- 
ferner liegende) Felder. Diese sich iiber das | gestellt. Diese Trennung 
ganze Tragwerk hin erstreckende Kopplung 

hat ein lineares Gleichungssystem von verhaltnismaBig 

hoher niedriger 

Ordnung zur Folge. Die Systemmatrix ist 

stets symmetrisch und oft von Streifenform. nicht symmetrisch und im allgemeinen voll aus- 
gefiillt. 
Nachdem dieses Gleichungssystem gelést ist, 


liegen lediglich die eliminierten GréBen vor. | sind nach einiger weiterer Rechnung Kraft- 
und DeformationsgréBen bekannt. 


Die Vorziige des Reduktionsverfahrens sind somit folgende: 

a) Es entfallt der lastige EKliminationsvorgang. 

b) Die Ordnung r des zu lésenden Gleichungssystem ist niedrig. 

Sie ist insbesondere unabhangig vom Grade der statischen Unbestimmtheit und der Felderzahl. 

c) Die gegebenen technischen Daten werden direkt oder doch in einfachsten Kombinationen in die 
Rechnung eingefiihrt. 

d) Die praktische Rechnung besteht allein aus dem Bilden skalarer Produkte, ist deshalb ein- 
férmig und somit leicht programmierbar. 

e) Man gewinnt Kraft- und DeformationsgréBen gleichzeitig, was besondere Vorteile bei der 
Ermittlung von EinfluBlinien bietet. 

Obgleich die vorgerechneten Beispiele genug Uberzeugungskraft besitzen, ist selbstverstandlich 
ein moéglichst objektiver Vergleich des Rechenaufwandes erforderlich. Ein solcher Vergleich lat 
sich exakt jedoch nur fiir Verfahren durchfiihren, deren Operationen von vornherein abzahlbar 
sind, also nicht fiir Naherungsverfahren wie etwa das von Cross und andere, da deren Rechenauf- 
wand in groBem Mafe von der Erfahrung und der Geschicklichkeit des Rechners abhangt. Mehrere 
zum Teil recht umfangreiche und deshalb hier nicht wiedergegebene Vergleiche zeigen in der Tat, 
daB das Reduktionsverfahren in fast allen praktisch vorkommenden Fallen weitaus weniger Ope- 
rationen (Aufschreiben von Zahlen, Additionen und Multiplikationen) erfordert als die Elimi- 
nationsmethode, selbst wenn man annimmt, da der Eliminationsvorgang ein fiir alle mal formel- 
maBig durchgefiihrt wurde und daher keine Mehrarbeit verursacht, wie etwa bei den bekannten 
Drei- bzw. Fiinfmomentengleichungen und anderen [siehe auch I Ziff. 11]. 


13. Die vollstandige Theorie des Reduktionsverfahrens. Zum SchluB wollen wir uns ganz von 
der speziellen Bedeutung des grundlegenden Differentialgleichungssystems II (3) lésen und nur 
noch voraussetzen, daB es linear sei (Biegen, Knicken, Schwingen)*3, elastische Bettung, ge- 
krimmte Rahmen* und dergleichen). Der Rechnungsgang 1a8t sich dann in groBen Ziigen etwa 
folgendermaBen umreiBen: 

1) Das dem Problem zugrunde liegende lineare Differentialgleichungssystem der Ordnung 2 0 
wird gelést; man erhalt ein Fundamentalsystem mit 2 @ willkiirlichen Integrationskonstanten und 
eine Partikularlésung. 

2) Diese im allgemeinen mechanisch sinnlosen Integrationskonstanten werden ersetzt durch die 
Kraft- und DeformationsgréBen am linken Feldrand, was durch geeignete Linearkombination 
geschieht?. 

3) Damit ist auch der Zusammenhang zwischen den Kraft- und DeformationsgréBen an beiden 
Feldrandern hergestellt: 


b= ih. r=(0,61), (23) 


wo D baw. f die zueinander konjugierten Deformations- bzw. Kraftgré®en sind. Falls am rechten 
Ende Federreaktionen angreifen, wird die zunichst quadratische Matrix % zu einer rechteckigen 


1 K. Marguerre, Math. Phys. 35 (1956) S. 28. 
® E. Pestel und G. Schumpich, Schiffstechnik 4 (1957) S. 55. 
8S. Falk, Abh, Braunschweig. Wiss. Ges. 9 (1957) S. 1. 


4 G. Schumpich, Beitrag zur Kinetik und Statik ebener Stabwerke mit gekriimmten Staben. Diss. Hannover 


1957. 
5 S. Falk, Abh. Braunschweig. Wiss. Ges. 7 (1955) S. 74. 
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Leitmatrix mit 2@ +1 Zeilen und 30 +1 Spalten erweitert; die Vorschrift (23) ist dann als ver- 
schrankte Multiplikation aufzufassen. 
4) Jetzt haben wir drei verschiedene Typen von Tragwerken zu unterscheiden: 
a) Durchlauftrager. Fir n aufeinander folgende Felder gibt (23) 


tn = nn es aia ay f= Ete Yo: (24) 

b) Verzweigte offene Tragwerke. Diese werden durch einen einfachen Kunstgriff auf einen oder 
mehrere Durchlauftrager zuriickgefiihrt (Koppelfedern). 

c) Geschlossene Tragwerke. Aus diesen lat sich stets eine Anzahl s von durchlaufenden 
Strangen herauslésen, deren Gesamtheit man nach geeigneter Erweiterung des Vektors y und der 
zugehérigen Leitmatrix ebenfalls wie einen einzigen Durchlauftrager behandeln kann. 

5. Auf jeden Fall sind damit saimtliche tiberhaupt in der Rechnung auftretenden Kraft- und 
DeformationsgréBen nach (24) lineare Funktionen der im Anfangsvektor Y) vereinigten o bzw. 
r=0, +°::-+ 0, Paare von Integrationskonstanten. Aus jedem dieser Paare yverschwindet 


genau eine Gré8e infolge der r Randbedingungen links; die iibrigen gehen als FreigréBen (A,, 


A,, ... A,) = a durch die ganze Rechnung hindurch. Die r Randbedingungen rechts liefern zu 
ihrer Bestimmung das lineare Gleichungssystem 
Wa+tr=o0. (25) 
6) Hier unterscheiden wir drei Falle: 
a) Ist || 0, so folgt a=— Wr. 


b) Falls ry =0 ist (homogene Probleme), stellt |9{| =0 als notwendige Bedingung fiir nichttriviale 
Lésungen eine Bestimmungsgleichung fiir gewisse in der Rechnung mitgefiihrte Parameter dar. 
Die Wurzeln dieser Gleichung hei®en Eigenwerte (kritische Knicklasten, Eigenfrequenzen usw.) 
und die zugehérigen EKigenvektoren q sind bis auf einen willkiirlichen Faktor ebenfalls bestimmbar. 

c) Wenn || =0, aber r40, so sind Lésungen nur méglich, wenn r gewissen Bedingungen 
genugt. 

7) Mit a zugleich liegt nun der Anfangsvektor yy und damit die Gesamtheit der Kraft- und 
DeformationsgréBen an allen Feldgrenzen fest. Das grundlegende Fundamentalsystem liefert daraus 
auch den funktionalen Verlauf der an der Deformation beteiligten Kraft- und DeformationsgréBen 
tiber dem ganzen Rahmentragwerk, womit das Problem erschépfend gelést ist. 


14, Zusammenfassung. Wir haben gesehen, wie sich auf Grund der Gleichungen II (3) oder 
irgendeines anderen linearen gewohnlichen Differentialgleichungssystems und einiger zweck- 
maBig eingefiihrter Begriffe wie Reduktion, Feldgrenze, Koppelfeder, Strangschar usw. die groBe 
Mannigfaltigkeit ebener und raumlicher Rahmentragwerke theoretisch erfassen, aber auch zahlen- 
maBig berechnen 1aBt. Insgesamt acht vorgefiihrte Beispiele und ein allgemeines Programm fiir 
digitale Rechenautomaten zeigen den praktischen Wert des Verfahrens auf; seine Uberlegenheit 
iiber die bis heute gebrauchlichen Eliminationsmethoden beweist ein objektiver Vergleich. 


(Eingegangen am 13. Mai 1957) 
Anschrift des Verfassers: Dr.-Ing. Sigurd Falk, Braunschweig, Celler Str. 93a 
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Nichtlineare Schwingungen bei zwei F reiheitsgraden 
IT. Mitteilung 


Von C. Schmieden 


Heinrich Graf zum sechzigsten Geburtstage gewidmet 


1. Einleitung. In der vorliegenden Arbeit werden die erzwungenen Schwingungen fiir das 
eleiche System von zwei Freiheitsgraden und nichtlinearen Riickstellkraften, aber unter Vernach- 
lassigung der Dampfung behandelt, dessen freie Schwingungen in einer friiheren Arbeit gleichen 
Titels — im folgenden mit I zitiert' — untersucht wurden. 

Um alle Komplikationen der hier méglichen Schwingungsformen zu erfassen, ist es notwendig, 
die Approximation der Liésung ebensoweit zu treiben, wie das in I der Fall war, da es nur so méglich 
ist, die Frage nach der Konvergenz der Fourier-Entwicklung der Lésung wirklich prazise zu beant- 
worten, Denn wenn man auch auf Grund allgemeiner Existenzsatze fiir hinreichend kleine Werte 
des Koeffizienten § der nichtlinearen Glieder und der Amplitude der auBeren Kraft der Existenz 
der Lésung sicher sein kann, so darf man doch aus diesem Tatbestand nicht ohne weiteres die 
Berechtigung herleiten, sich bei einer Naherungslisung nicht einmal um die Gré®enordnung der 
nicht beriicksichtigten Glieder zu kiimmern, insbesondere dann nicht, wenn die Amplituden der 
auBeren Kraft oder des Ausschlages selbst Werte annehmen, die groB gegen § sind. Um aber die 
erforderlichen Rechnungen nicht unnétig zu komplizieren, behandeln wir im ersten Abschnitt zu- 
nichst die sehr viel einfachere Differentialgleichung von Duffing fiir einen Freiheitsgrad, und zwar 
in einer wesentlich weitergehenden Approximation, als das bisher geschehen ist, mit dem Ziele, 
das Konvergenzverhalten der Liésung dieser Differentialgleichung voll aufzuklaren. Bei der Lésung 
des Systems von zwei Freiheitsgraden werden wir uns dann auf diese Rechnung beziehen kénnen, 
soweit sich die Analogie zwischen einem und zwei Freiheitsgraden durchfiihren laBt. Das geschieht 
im zweiten Abschnitt. Der dritte Abschnitt bringt dann die Untersuchung derjenigen Phanomene 
bei zwei Freiheitsgraden, die keine Analogie fiir einen Freiheitsgrad besitzen, also das eigentlich 
Interessante der Untersuchung darstellen. Sie treten dann auf, wenn sich die Kigenfrequenzen fiir 
kleine Schwingungen annahernd wie 1:3 verhalten — allgemein, wenn dies Verhaltnis 1:2 -+- 1 ist. 
Denn dann steht die erste Ultraharmonische — allgemein die n-te — der Oberschwingung in 
Resonanz zur Grundschwingung, ebenso die erste Subharmonische — allgemein die n-te — der 
Grundschwingung in Resonanz zur Oberschwingung. Fiir solche Systeme treten bei den erzwun- 
genen Schwingungen in der Nahe der Resonanzstellen nicht nur, wie tiblich, bis zu drei mégliche 
Lésungen, sondern in der Nachbarschaft der Grundfrequenz ebenso wie in der der Oberschwingung 
bis zu neun mégliche Lésungen auf. Auf die Frage der Stabilitat dieser Lésungen wird in dieser 
Arbeit nicht eingegangen, ebenso auch nicht auf die Modifikation dieser Lésungen durch Hinzu- 
nahme eines Dampfungsgliedes. 


I. Integration der inhomogenen Differentialgleichung von Dujfing. 2. Die Grundlésung. Wir 
gehen aus von der inhomogenen Differentialgleichung von Duffing mit verschwindendem Damp- 
fungsglied 

gtoiqtfpoarg = FoisnQt (1) 
und setzen voraus, da die Amplitude der erregenden auBeren Kraft endlich, also groB gegen den 
Parameter / des nichtlinearen Gliedes ist. Wir stellen zunachst die Grundlésung auf, untersuchen 
deren Konvergenzgrenzen, diskutieren die Ultraharmonischen, also die Lésung in der Umgebung 
der Stellen, an denen sich Q : wy wie 1 : 2 n + 1 verhalt, und zeigen schlieBlich, daB sich eine Ver- 
zweigung der Grundlésung immer dann berechnen Jat, wenn das Verhaltnis von Q : @, eine 
rationale Zahl darstellt, mit Ausnahme der den Ultraharmonischen zugeordneten reziproken un- 
geraden ganzen Zahlen. Verzweigungen treten also nicht nur fiir die Subharmonischen ungerader 
Ordnung, wie schon seit langem bekannt, auf, sondern prinzipiell immer, wenn Q und (Wo in ratio- 
nalem Verhaltnis stehen. Da aber die rationalen Zahlen auf der Zahlengeraden iiberall dicht liegen, 
gilt das auch fiir die Verzweigungspunkte: jeder Punkt auf der Amplituden-Frequenz-Kurve — 


 C. Schmieden, Ing.-Arch. 25 (1957) S. 292. 
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kurz mit »,Response-Kurve bezeichnet — mit einer rationalen Abszisse ist ein Verzweigungspunkt! 
Machen wir nun in (1) die Substitution 


EO tpl ws ot g! s 
so geht sie iiber in 
is . Qt 
te ree), (3) 
Wir setzen nunmehr an 
Q 2 
a =1+x, q = 9, sint + O,sin37-+--- (4) 


und denken uns x und die 0; in Potenzreihen nach dem als klein vorausgesetzten Parameter 
entwickelt : 


Y+e#e=14+%4,1%,8 +p? +°°-, OG. = A, 
O; = Oso + Oy B +°*°, O; = Boo + 85, B + 0° *- 


A iibernimmt dabei die Rolle des Parameters. 


(4a) 


Dabei wird sich zeigen, daB von Ausnahmebereichen abgesehen, die @,,,4, proportional zu f” 
werden. 


Wir haben nun zunichst die dritte Potenz von q zu bilden, wobei wir nur die Terme angeben 
wollen, die im folgenden an irgendeiner Stelle wirklich gebraucht werden. Also [vgl. I, (9)] 


4q? = Q, sint + Q,sin3tT +:°:-+ Qy)sin9Tt (5) 
mit 
 Q =3(G— O26, + 2.0, 6:26, O, O; + O20, — 62 ,), 
Q, =— 61 + 3 (2 10, + O1— 61 6, +20, 0,0, —20; 0, 0, —016,), 
0; = 3 (— @? O, + O, G2 + 2 6? O, — 6? O, + 2 O, O; O, —2 0, 0; 9, +2020;), (6) 
0, = 3 (— 0, 61020, + 260, — 62, +20, O, 0, +20, 05, +2030), 
0, =— 62+ 3(— 6? 6, + 26? 6, —2 0, 0; O; + 2 O, O; O, + 2 2 Oy). 


Ist nun @, klein gegen O, von der Ordnung 3, so sieht man aus (6) sofort, daB dann @,; von der 
Ordnung f/f? usw. wird. Also setzen wir an 


0,=A, 0, =P 95, O;, = POs, O, = §P 0,, Og = * 0; (7) 
dann reduzieren sich die Q,, wenn man Terme der Ordnung /? in ihnen streicht, auf 
0, =3 (44 —p 479, + 2 fp? AP), 
0, =(— 4? + 6f AP, — 3? APY), 
Q; =3(—f A? 0, + B? A OF + 2 f? A? D5) (8) 
0, =3 (— A? 0,— fA 09), Ova 0i 


Damit erhalten wir durch Koeffizientenvergleich nach dem Sinus: 


[sin T] : (l+x)A+A4 PEO = Bs 


B 
[sin 3 7]: —9(1 + x) BO, + BO; 4 oD. 


| 
B c 
[sin 5 T]: — 25 (1 +x) p20, + B29; pie | 


Peas 9490 bape, eo, + £0, =0. 


112 C. Schmieden: Nichtlineare Schwingungen bei zwei Freiheitsgraden Ingenieur-Archiv 


In diesen Gleichungen machen wir nun Koeffizientenvergleich nach f und erhalten damit die end- 
viiltigen Bestimmungsgleichungen fiir die unbekannten Koeffizienten: 


F 
[6°]: eae also De a 
ae aur , A8 
[p"]: =" Aq AO), (8 + 9%) Ys = — > 
ay 
Hg 3 A? bso = 0, — (8 + 9 x) Js, — 9 % Pso + AP ie =e 0, 
[3"] 
— (24 + 25 x9) P59 — 4 A? Os) = 0, (10) 
( 3 bite 
— 1, A + 4 (2.4 0, — A?) = 9, 
G 3 A? 9 q 
— (8 + 9 x) P32 — 9 #1 Dg, — 9 He V39 4 A (2 Vs: — Vso) = 0, 


[6] : ) B 
eos (24 +- Ds 0) er —— 25 4 Ves + Za A (O%, ee A ee -++ y, A Dee) — 0 ° 
ah 
U 


3 oy 4 
48 -|- 49 Ha) es ey A (oe. ae A Oso) a 0 Y 


Fiir die Berechnung von x, ist, wie ein Blick auf die allgemeinen Formeln (6) fiir die Q; zeigt, nur 
die Kenntnis der aus der letzten Gleichungsgruppe anfallenden Koeffizienten erforderlich: 


3 
4 Al | vy ( Aa Ae —+ 4, A ag Onn ——— » vA Don Men) = 0 . (11) 
Die Auflésung dieser Gleichungen (10) und (11) liefert nun 
i = 3 2 9. —_ A? 
eee ee Ms ery te ies 
ae AS el AS gid ne Sere 
“2 — 16° B49 ms? V31 “16 | (8+ 9%)?’ 50 16 © (8 + 9%) (24 + 25%) ” 
Pena: 4 the ve RR ea ee (12) 
U3 64 (GAS Ox)? ? 32 32 (849% 5)® 64 (8 + 9%)? (24+ 25 x) ’ 
45 AP 23 ; 1 . 
ae 256 |(8-+ 9x)? | (8+ 9 m%)? (24-4 25%) | * 


Die Angabe von ?;, und #9 ist iiberfliissig; man bendtigt sie erst fiir x,. Um die Frage der Konver- 
genzgrenze unserer Entwicklung zu beantworten, geniigt aber bereits die Kenntnis der bisher be- 
rechneten Terme, da sie, ebenso wie der Bau unseres Gleichungssystems, eindeutig erkennen lassen, 
daB diese Entwicklung divergieren mu8, wenn einer der Faktoren (8 + 9 x9), (24 + 25 %), ... in 
die GréBenordnung von 6 kommt. Wir wollen die Gleichung der Response-Kurve anschreiben: 


(2) <1 F Dom © p? AS 57 p AS 
i pes Va 


16 (849%)  64(8 + 9%)? 
Sioa eee ee 1 a 
256 |(8 +9 x)® | (8+ 9 x) (24 + 25 v5) zs 


Wo 


(13) 
mit 


q=Asint +f %,sin3t + 670; sin 57 + 630, sin 77 +---. (14) 


Diese Formeln stellen offenbar gut konvergente Reihen, selbst fiir endliche, nicht zu gro8e f, dar, 
solange A hinreichend weit von den kritischen Werten 


on = 
4,=5F, §4,= oF... 
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entfernt ist. Setzen wir daher in der Umgebung der ersten kritischen Stelle 


AS] Paya hy + x6 , (15) 
so erhalten wir 
ae ot 8 x VB 1 eee: «ye \ 3 Gee eene fo) 
& =o 19 7, ale * 1 PR (1 FE) 138 x ( ae nr aa) 
ae 


0 


und eine ahnliche Entwicklung von q, von der wir nur die beiden ersten Glieder anschreiben: 
=\. | Ve F xVB\ . 21 Vp F3 «/B \ 
q=(Fy+xyp) snr aoe 14 rae sin 3T ieee hap +e (17) 


Beide Entwicklungen versagen offenbar fiir hinreichend kleine Werte von x, da namlich in der 


: : w ae mnie 
Nahe der Resonanzstelle der ersten ungeraden Oberschwingung Q ~ ce der Koeffizient von sin 3 T 


von der gleichen GréRenordnung wie A wird und damit unsere Grundannahme iiber die GréBenord- 
nung der Fourier-Koeffizienten hinfallig wird. 

3. Die Response-Kurve in der Nahe der Resonanz der Oberschwingung der 
Ordnung 3 (erste Ultraharmonische). Wird, wie es hier offenbar der Fall ist, @, endlich, 
so treten in der Entwicklung von q? bis hin zu sin 9 t endliche Terme auf; diese Glieder sind also 
zu beriicksichtigen, wenn man die Lésung auch nur bis auf die Glieder der Ordnung f einschlieB- 
lich bestimmen will. ZweckmaBigerweise laBt man jetzt den Koeffizienten 3, die Rolle des freien 


. 8 
Parameters iibernehmen, setzt also mit F = > F, an 


Ol ie pty es 0, = B, O,=B 99 +°:°, O, =P Yo +--+ (18) 
und 
l4+xux=1l4+u4i+%4B+%,B?+°--. (19) 
In den (9) entsprechenden Gleichungen sind jetzt fiir die Q; die vollstandigen Ausdriicke (6) einzu- 
setzen, also z. B. 


: 8 
—(1+~)0,+6,++6(O—BOi +2 BO,—2 BO, 0, + BO,— BO, => F,. 


Der Koeffizientenvergleich in f liefert jetzt 


* 8 * 8 seein 
s [—“F=5 Fh. also Po aes Me aden 
a |-9+%)B+B=0, ako —B+B=0, 
3 2 
— #5 Oy —%41 Fy + 7 (hi BR +2 BR) =0, 
Poms (Fis oR Fs +33) = 0, (20) 
1 peer Ae ee ree I) = 0" 
[ps g V50 50 1 4 0 0 
49 B) oe 
6 On + Or A Bio =0, 
1 
ee tai a0 i B= 0 
mit den Lésungen 
1 
1 = 6p (3 B+ 6 BY — Fi). 
9 27 
On = 4 Bom — 35 (Fi BR + 2 B F,) ; (21) 
27 sears, Pee or 
Deo = 64 B Fy (B— Fy) > Dx = — j69 BP Fo> Don = 39 B?. 
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Damit werden se Glieder in f? berechenbar; wir beschranken uns auf 
Bry = [— F204) + 4 BF, I, — Fy O59 + 2 B Fo (O59 — O79) — B Doo] - (22) 


Setzt man alles ein, so erhalt man schlieBlich: 


il ; 
ue = 3-97 pe F2(4.B— F,) 2U BF, 2) B eo Be) 
+1) BF (B—F,) (2 B— F,) +> Bi Fs + Be. (23) 
Unsere jetzige Lésung ist natiirlich nur wieder zu gebrauchen fiir Werte von B, die groB gegen p 
sind. Setzt man also B = —y//f , so gilt in erster Naherung 
F Honea eae 24 
Also fiir 
a= Vy yo = ay 
wird we 
ial oa 6 Wis 
q = (Fy + xB) sint—y/B sin 3+ =(F, + xB) sine outs sin 37, (25) 


was in dieser Naherung mit (17) iibereinstimmt. Die Konvergenz-Bereiche beider Entwicklungen 
iiberdecken sich also, so da® jede innerhalb ihres Konvergenzbereiches als die Fortsetzung der 
anderen benutzt werden kann. 

Ist B endlich und gro gegen Fy, so ist die Amplitude | q |max wesentlich durch B gegeben, so 
daB wir als Response-Kurve erhalten 


(2 SI Bhs 1a pe See (26) 
Wo 
also als Skelett eine Hyperbel mit den Achsen 
Eyer a eee 
a = 1 432 re b= rar eae Fe (27) 


Fir groBe positive B verlauft die Response-Kurve links von dieser Hyperbel, fiir groBe negative B 
rechts von ihr. 


Kommt B in die GréBenordnung von Fo, so mu man bei der Amplitude auch das Glied mit 
sin T mitberiicksichtigen. Fiir negative B 
B=——2 Fy (x20) (28) 


wird | q |max fiir 7 = 27/2 erreicht, da sich hier wirklich alle Terme der Fourier-Reihe addieren. Es 
kommt 


we (Ge Weer Coes ieee ae oe 

q= Fi (l+*) +559 e eee a x) + [B?] . (29) 
Bestimmen wir das Minimum dieser Funktion, so erhalten wir fiir x 

dq a oe i S21 3 Pad 5 BVA os Pol 

a OF also Bes 5 9% 14 B r* == (29a) 


'B 
mit dem annahernden Wert der Wurzel x ~ — Fol Bb 
AND 
Dies ist also die Stelle, wo die Response-Kurve zwischen cy und w,/3 ihr Minimum erreicht. Zu 
diesem Wert von x gehéren 


Gnin = Fy ee [A]. 2= 9+ 518 +e. (30) 


Zu diesem Minimum yon q gehort der Wert 


6, =¥ + EP +A. (0a) 
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Da nun offenbar die Konvergenzgrenze der Entwicklungen (16) und (17) an der Stelle des Mini- 
mums der Response-Kurve zwischen 9/3 und , liegen mu, folgt, daB mit 


81y2 _, 
es, ey 
diese Entwicklungen konvergieren fiir | x essa cas 

Der Verlauf der Response-Kurve in der Umgebung der Stelle Wo /3 ist also qualitativ der gleiche 
wie an der Hauptresonanzstelle «,, nur daB die Schwingung wesentlich mit der dreifachen Frequenz 
der auBeren Kraft vor sich geht (erste Ultraharmonische). Jedenfalls tritt aber keine Verzweigung 
der Lésung auf. 

Ganz Entsprechendes gilt fiir die hdheren Resonanzen ungerader Ordnung. Ein ganz anderes 
Verhalten liegt bei den Resonanzen der Oberschwingungen gerader Ordnung und bei den Subhar- 
monischen aller Ordnungen vor: hier ergibt sich jeweils eine echte Verzweigung der Lésung. Wir 
wollen das hier nur fiir die Subharmonische der Ordnung 1/3 zeigen, weil hier die Rechnung am ein- 
fachsten verlauft, und anschlieBend fiir die Oberschwingung der Ordnung 2. 

4. Die Subharmonische der Ordnung 1/3. Wir setzen in (3) jetat 

O\2 
- . ie Sonera (32) 


\o 


Xmin & 


Wy Vl +x 
die auBere Kraft ist also durch F'sin3t gegeben. Um die Subharmonische zu erhalten, setzen 
wir an 


O, =A, Os =U +B On °°? O;,=B Os +°*:s Oy =B Io +°-*- (33) 
Damit erhalten wir als Gleichungen nullter Ordnung in [3 
—(l+%,)4+A=0, also %yA=090, | 
— 9 (1 + x9) Pq + Ogq = F. J 
Soll A = 0 sein, so muB x) = 0 und damit 4, = — F/8 sein. 


Die andere mégliche Lésung A = 0, x) + 0 liefert eine Amplitude, die mit einem Gliede pro- 
portional sin 3t beginnt. Man erhalt hier 


: (34) 


FEES DO, reel ae 
— (8 + 9x) Ion = F also % =—— na 


also 


Das ist aber gerade der Anfang der Grundlésung, nur in anderer Bezeichnung: A ist durch Uso, 
t durch 3 T ersetzt. : 
Fiir die Subharmonische, also A +: 0, folgen nun als Gleichungen der Ordnung f 


Bd Ba ( A Diy 2A 02) = 05 


1 : $06 , ‘ 
Ua 9 ey Ug 4 (— 4? +6 Ab + 3 Bj.) = 9. ¢ (35) 
q 3 9, 1 . — 
= esis, (O59 — A) = (Ny, — 48 On — | AP, == (I) — 80 og, — A C0 | 
und in der Ordnung /* 
Ko, AS £ ( Ae Igy aA D9 Dz, — 2 A Vso Iso a lass (P50 coc Da)) =9 (36) 
mit den Lésungen a 
AF aN AF? pete 
Iso Gril (ae + 8 A) ? Dao ag ae Pee Joo Ss? | 
J 299 43 
On, = gy (51 F3 4 21-28 AF + OA Dede). | = 
ty = oy (P+ 4A F 4 324%), 
AF? FS 
m2 [by 44+ A 44e (F+8 4) + 5764 +3 P| 
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Fiir A > 0 wird die maximale Amplitude fiir tT = 7/2 erreicht, also 


g=At=4+i61, A=q-F+lA- 


ies , Oh 
Setzt man diesen Wert von A in die Entwicklung von ical ein, so kommt 
0 


(By =9 +P —4gF +320) +1, 9 > 9, (38) 
0 
was sich umformen lat in 
Q \? 3p Fr WM Bp 
Cie nea ” 


Diese Gleichung stellt eine Hyperbel dar, deren Scheitel bei gq = F/16 liegt und die Q-Achse bei 
i= 3 0, nnedee Die Verzweigung tritt fir A = 0 ein. Hier gibt unsere jetzige Lésung 


Gielen ey 


Wy 
(40) 
51 rales Fr 51 Dm 3 ei 
q=—(F— SP r)singx = — (232 ro) in (3 elas Li i), 
waihrend die Grundlésung mit 
F 51 fp F Bley = 
ety ate apie (41) 


fiir die Frequenz liefert 


ae g Re 3p 
Gq) lage tgeatal tee Se Pla gaol ta) 


Wo 


und fiir die Amplitude 
q— At sine = A* sin(3 le +38 1). (43) 


Die Gleichungen (42) und (43) stimmen also jeweils mit (40) tiberein. Ubrigens gibt es, worauf hier 
nur kurz hingewiesen werden soll, noch zwei weitere Subharmonische der Ordnung 1/3, die aber bis 
auf eine Phasenverschiebung von + 22/3 mit der hier behandelten Fikes Sheet a. 

5. Die Verzweigung Ths Q & @,/2. Die Rechnung verlauft prinzipiell genau so wie im 
letzten Abschnitt. Man mu nur in der Fourier-Entwicklung der Amplitude jetzt valle Glieder, also 
auch die mit geradem n mitnehmen, also ansetzen: 


q = 9, sint + O,sn27+°°-. (44) 


Endliche Terme enthalten @, und @,, Terme erster Ordnung die folgenden Glieder bis n = 6. Als 
freien Parameter wahlen wir @, und setzen es gleich B. Dann liefert die iibliche Rechnung fiir 


Ot n= BB +2 ie ale [B2] » 


3 
Ones 1 ti ; AF mn 
be ~ 4 i ace C3 Ae Teel 
Im Verzweigungspunkt ist B = 0, also erhalten wir die Amplitude 


A* = é F(1 = a cine (46) 


(45) 


und die Frequenz 


Q? 1 | 3B AF? | 
Se TAR} & es: - 
Setzen wir A* als Wert von A in (13) ein, so erhalten wir 


SIP eEAD eek ty, 


was mit (47) iibereinstimmt. 
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6. Die Verzweigung bei 2 = 3/2 @). Um zum SchluB noch zu zeigen, daf man fiir jeden 
rationalen Wert von Q :q, eine Verzweigung berechnen kann, da8 also prinzipiell die Verzweigungs- 
punkte auf der Grundlésung fiir die Response-Kurve iiberall dicht liegen, wollen wir hier den fiir 
solche Verzweigungen einfachsten, aber typischen Fall soweit durchrechnen, da® der in diesen 
Fallen einzuschlagende Weg klar zu iibersehen ist. Wir wahlen dazu Q — 3/2 wy) und setzen an 


C= Osin y + Onsint +O, sin 51+ +++ 4+ Oysin a7 fos, 


Itealtmtmbt, (Sy =(z) ata). 


Wo 


(49) 


Dabei werden, wie man aus dem Glied f q® in der iiblichen Weise erkennt, fiir die Verzweigungs- 
lésung O, und @, endlich, wahrend @, und O, bis O, auBer O, von der Ordnung f ausfallen; O, und 
alle tibrigen O; werden mindestens von der Ordnung f?. Um aber den Vergleich mit der Rechnung 

fiir die Subharmonische der Ordnung 1/3 leicht durchfiihren zu kénnen, geben wir in g? auch noch 
_ die Glieder an, die fiir ein endliches @, ebenfalls endlich wiirden. Setzen wir also (endliche Terme 
sind unterstrichen) 


eee = 0, 0, — 6 0. O20, 20, G2 4 +r -):, 
0,=3(2 6:0, + 62 + 20, 6, 0, +20, 02+ ---), 
Q, = — G} + 3 (O, 8} + 2 626, +2020, + OF +--+), 


3 (50) 
Q, = 3 (— OF O, + 2 O, O, Os + 0, ee) 
Cee © 07-6; — 0; 0, 4 O, OF + + +)s Qc =— §;— 60, 0; 0; ela 
Q; =—3 (020, +O, +---), Q=—30,08+--, @=—O4---,| 
so erhalten wir die Gleichungen 
} 
6; | 0, FQ=0, 
—(1 +x) 0, + 0,+40,=0, =n 
—F(1+4)0,+0,4+10=F, 
SAS) iO, Ox FQ, =0 
usw. 
Als Glieder nullter Ordnung in f folgen nunmehr 
== (1 = 42,) Pan Von = 0 also Hoag 0s 
(52) 


Sh + x9) Ogg + Ogg = F. 


Wiirden wir hier 0, = 0 setzen, so wiirden in unserem Ansatz alle @,; identisch verschwinden und 
wir damit auf die Gleichungen fiir die Subharmonische der Ordnung 1/3 zuriickkommen; denn man 
braucht ja nur statt t/2 wieder 7 zu schreiben, und fiir (Q/a,)* statt (3/2)? (1 + ) den Wert 9 (1 +%) 
zu setzen, so sieht man, daB beide Ansitze identisch sind. Also muB man die andere Initiative 
wahlen und setzen 


Ho — 0 > Don = B 9 (53) 


wobei wir B wiederum als den freien Parameter benutzen. Die weitere Rechnung verlauft nun 
vollig analog zu den Rechnungen der vorigen Abschnitte. Um die Verzweigung sauber auszu- 
rechnen, braucht man die Koeffizienten #,, und #,. Man erhalt 


4 


O30 in TS 


I ile Oi4 = — B O35 2 


3 7p: 2 
Bg, =— 35 Oy (B+ 1405), 4 = (B+ 283,)- 
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Um x, zu berechnen, wiirde man noch 5, und #7, bendtigen; aber um den Charakter der Verzwei- 
eung zu erkennen, geniigt bereits die Kenntnis von x. Es gilt also 


(ey ap tae 2Cs)). 


8 
q =P Oy, sin + + B sin t + (P39 + f Ys;) sin 9 + [6], 


(55) 


also die Amplitude 


q~B fir groBe B, (56) 


also die Amplituden—Frequenz-Relation 


Ga) =1+48[e4 a(t, (57) 


3 Wo 


was sich schreiben laBt 


Gale la (58) 


3 W 
Wir erhalten somit die iibliche Hyperbel, rechts von der kritischen Stelle Q = 3/2 w. Die Verzwei- 
gung tritt ein fiir B = 0, wo unsere jetzige Lésung liefert 


Seer | © 


Fir die Amplitude erhalten wir in diesem Falle einschlieBlich der Terme in 6 


21 
1 = P2(1 —7p 6 9) (60) 
Setzen wir diesen Wert fiir A in die Formel (13) unserer Grundlésung ein, so erhalten wir 
4 21 ,/4F\? 
+= Bel 


Spare eBay eae to a 
— Sete Gls ta Spe setey en. 


was offensichtlich mit (59) tibereinstimmt. 


II. Die Grundlésung bei zwei Freiheitsgraden. 7. Die Response-Kurven fiir b=0. Wir 
wollen in diesem Abschnitt den allgemeinen Charakter der Grundlésung fiir das gleiche System von 
zwei Freiheitsgraden untersuchen, dessen freie Schwingungen in I behandelt wurden. Wir legen 
also das durch eine auBere Kraft modifizierte System I, (1) zugrunde: 


P+e2G +028 4 eau (62) 


I++ opp + Borg =0. 


Wir wollen aber die Grundlésung fiir dies System nicht soweit durchrechnen, wie wir es im ersten 
Abschnitt fiir die Duffing-Differentialgleichung getan haben, sondern uns jeweils auf die ersten 
Glieder beschranken, aus denen man dann in Analogie zu einem Freiheitsgrad Schliisse iiber die 
Konvergenzgrenzen ziehen kann. Alles Wesentliche kann man schon aus der Responsekurve fiir 
kleine Schwingungen erkennen, die wir jetzt zunachst bestimmen wollen. Fiir unseren allgemeinen 
Ansatz [vgl. I, (6) und (7)] sind das die Glieder der Ordnung f°. Setzen wir wieder 
> Q\2 
t=OptVl +x, (5) —it* (63) 


dann wird aus (62) 


(64) 


(1 + x) (0" +079") +0468 = Fsint, 


(1+) (O" +o) +o +Bq? =0. 
Mit 


eis ae aia nearest (65) 
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erhalten wir dann fiir die Glieder der Ordnung p°® 


Xo ie rate a2 (1 -|- ea) A ++ F = 0 “| 
6 
(1 + 29) Po + % A =0,f (66) 
mit den Lésungen 
Tas ¢ 
soe es 0) 7? I =—] = A= ee. (67) 
We (ea 5) 2 1+ x we — a2 (1 +- x9)? * 
Der Nenner verschwindet, fiir die beiden Eigenfrequenzen 
1 1 
1+ %q = oe 1 + X92 = re (68) 
A hy 
| 
| 
| 

| 

je 2 * ! 82 

= z 2 ty age 2 

Vita Vink Vra ® YH 1 V7-a D, 

| 
| 
| 
Abb. 1. Der Verlauf der Amplitude A als Funktion von 2 : w,. Abb. 2. Der Verlauf der Amplitude 9,, als Funktion von 2 : @». 


Der Verlauf von A und #,, iiber 2/c, ist in Abb. | und 2 aufgetragen, aus denen wir fir das folgende 
entnehmen, daB zwischen den beiden Eigenfrequenzen A negativ, auBerhalb dieses Bereiches posi- 
tiv ausfallt, und da das Minimum von |A| fiir 


. 1 
oe ieee, (69) 


erreicht wird. Da wir aber fiir die Responsekurve im nichtlinearen Fall 1 + x als Funktion von A 
benétigen, eliminieren wir aus (66) #,) und lésen nach 1 + x, auf. Wir erhalten 


tema alreeyerey—i+ | fm 


Das Plus-Zeichen vor der Wurzel bestimmt den Bereich, in dem x) > x; ist, also wesentlich die 
Umgebung der Resonanz bei der Oberschwingung, das Minuszeichen, fiir das also x) < x, gilt, 
liefert den Bereich in der Umgebung. der Resonanz der Grundschwindung. Der Wert x) korre- 
spondiert dem Minimum von | A |, also 


1 F F — 
LS An = aos wegen 1 + — == el an, 71 
cee =a 2(1—1—2o?) i a ae y fy 

8. Die Grundlésung im nichtlinearen Fall. Wir bestimmen jetzt mit der gleichen Me- 
thode wie im ersten Abschnitt und in I die Anfangsglieder der Grundlésung soweit, wie wir sie 
brauchen, um die Konvergenzgrenzen dieser Entwicklung festzustellen. Wegen der weitgehenden 
Analogie zu dem Falle eines Freiheitsgrades geniigen dafiir schon die Anfangsglieder. Die Bestim- 
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mungsgleichungen fiir die Koeffizienten proportional / werden 


5 x 3 x 3 
oy Oar -}- mA (08 — ee -+- 4 ee A8 == (), 


ty A. ay rk 
(1 + %) Hy + 1 i ee ee As = 0, 
(72) 
1 (_%_\" 49 = 0 
(8 + 9 29) O31 + 9 0? (1 + %0) Ps1 — A (i ea A =U, 
MD yp 
9 (1 + x9) Dg, + (8 + 9%) Par + A 4h = 0. 
Die Determinante der beiden ersten Gleichungen wird 
x 5 l 
2 at (Lt me) +0 = ppg, Ma 2 + Hs) — a? (1+ me) (73) 
Sie verschwindet, wenn 
1 * 
if a ae pms: —]1 ira (74) 


ist, also fiir die Stelle, wo die A-Responsekurve im linearisierten Falle eine horizontale Tangente 
besitzt. In der Umgebung dieser Stelle ist also die Entwicklung der Grundlésung nicht zu ge- 


brauchen (vgl. Ziff. 9). 

Die Determinante der beiden anderen Gleichungen (72) wird 

(8 + 9 x5)? — 81 a? (1 + 29)? = [8 + 9 % —9a (1 + %)] [8 + 9% +90 (1+ %)]- (7) 

Sie verschwindet, wenn 

1 
I. SP @9) 
wird. Diesen Nullstellen sind offenbar jeweils die ersten Ultraharmonischen der Grund- und Ober- 
schwingung zugeordnet. Auch hier versagt also, genau wie bei einem Freiheitsgrad, die Grund- 


lésung, und man muf} dann eine Entwicklung ansetzen, bei der endliche Terme in ®, und 0, auf- 


treten (vgl. Ziff. 10). Wir geben hier gleich noch die Werte von A und #,, fiir diese Falle an. Es gilt 


1 + % = oder 1 + 4% =—.— (76) 


1 9 F 9 84+9a 

Dee oe ay eee PORT err re 

fea ree Saleen 7 aia: ag NO See 7?) 
ETOYS ie , 164s" 10.0 tsk lied Sahin 


Aus den letzten Werten sieht man, da die Berechnung der Ultraharmonischen zur Oberschwingung 
fiir « = 4/5 in der iiblichen Weise nicht méglich ist. Mit diesem Ausnahmefall beschaftigen wir uns 
ausfiihrlich im dritten Abschnitt. 

Solange « < 40/41 ist, liegt 1/9 (1 —«) und erst recht 1/9 (1 + «) links von 1 + x%, so daB fiir 


solche Werte von « unsere Grundlésung sicher konvergent — fiir hinreichend kleine 6 — ist, wenn 
(1 + %) — (1 + x6) = % —% > Bs (78) 


hingegen hort, von der Stelle x* abgesehen, die Konvergenz jeweils in der Umgebung der Ultra- 
harmonischen, also fiir 


elie i a 1 = : 
P12 eigenen oh Mo = ay Cees co 
sicher auf. 
9. Die Entwicklung in der Umgebung von x. Wegen des Verschwindens der Deter- 
minante (73) muf an dieser Stelle eine Sonderentwicklung in der folgenden Form angesetzt werden: 
ONG l 1 
‘ (80) 


F F 
OS ene Q, = 2 (ev ee Pir pris 
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wobei jetzt Jj, die Rolle des freien Parameters iibernehmen m alt hier fiir die Gli 
uB. Ma halt | 
proportional f bzw. fh? die folgenden Relationen: PONCUDR LD eet Eagan odes 


| 44 (1 —yl—o? ) + op, — ae — 0?) — 3 GI VI =o 05 ) 

fa |: es (1 = —o? ) 3 he neil — 0) (1 ae yl oe ) 
tale) eee 

[6°] : is (9s, + 0? G1) + Og, — : eh = 0, 

pice Oa todo +4 (8) (pa) =e « 
PL 1 On + 0% on) +E Tat + 2 8, — 0g) =O, 

Se 4 Ou + om) — BE (1-4 =" 4, 

3 Fr? i 2 
=" (gaye) Cen tow =o. J 


In den beiden ersten Gleichungen verschwindet die Koeffizientendeterminante von x, und #,,, hin- 
gegen nicht die entsprechende Determinante zu g,, und x,, was erneut zeigt, daB 0, der freie 
Parameter sein muB. Die Auflésung dieser beiden Gleichungen liefert 


3 F? 1—(1— yl—@ )? 


LE i Tier eg peek 
* hae" 82 
pues ee ee ee i ee syle - 
F 1— 2 32 yl geet (= ji—@)? g 


GY, hangt also bemerkenswerter Weise nicht von #,, ab. Ein Blick auf die restlichen vier Glei- 
chungen (81) zeigt, daB ~,. von x, ,, und damit nach (82) quadratisch von 7, abhangt. Bestimmt 
man #7 aus der Relation dg,,/di,, = 0 und tragt diesen Wert in (82) ein, so bekommt man bis auf 
Terme in /? die Koordinaten des Punktes in der (A, {2)-Ebene, wo die A-Responsekurve im nicht- 
linearen Fall eine horizontale Tangente besitzt. Da aber g,, unabhangig von 7,, ist, kénnen wir 
den Wert von A, fiir den die Konvergenz der Grundlésung fiir x) © xj aufhért, auch ohne diese 
Rechnung bis auf Terme in /? angeben: 


1 AE ee Sp Bi ee) 
(a Se GG ara alee en ea or (83) 


Unsere jetzige Entwicklung ist aber auch noch fiir endliche, sogar groBe 0, brauchbar, so daB sich 
die Konvergenzbereiche dieser Entwicklung und der Grundlésung iiberdecken. In dem Uber- 
deckungsbereich stellen beide Entwicklungen nur Umordnungen voneinander dar, was man in der 
gleichen Weise zeigen kann, wie das bei der Diskussion der Duffing-Differentialgleichung in Ziff. 3 
gezeigt wurde. Wir kénnen daher hier auf einen solchen expliziten Nachweis verzichten, zumal das 
ja aus bekannten Satzen der Reihenlehre unmittelbar folgt. 


10. Die Lésung fiir die Ultraharmonische zur Oberschwingung. Man kénnte nun, 
ebenso wie bei einem Freiheitsgrad, die Verzweigungen der Grundlésung fiir die beiden Klassen 
von Subharmonischen zur Grundschwingung einerseits, zur Oberschwingung andererseits, be- 
rechnen, ebenso zeigen, da solche Verzweigungen auch schon dann auftreten, wenn die Fre- 
quenz der auBeren Kraft zu einer der beiden Eigenfrequenzen in rationalem Verhaltnis steht. 
Aber da bei solchen Untersuchungen prinzipiell nichts Neues gegeniiber einem Freiheitsgrad 
herauskommt, wollen wir darauf verzichten. Nur auf die erste Ultraharmonische zur Ober- 
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schwingung wollen wir noch kurz eingehen, um das in solchen Fallen einzuschlagende Verfahren 
zu skizzieren. Wir setzen dazu an (a =— 4/5) 


1 
Lee xn = 9(1—2) a aa eg aa 
D, =P + 9uh +°°'> 0, = 9 + up +°°* (Pro Pio aus (77)) , a 
Pz = P30 > Os, 055 Se Osis = 


und verwenden go als freien Parameter. Die Glieder der Ordnung (3° liefern die Gleichungen 


! ¢ 1 
~~ 9(1—a) (P19 + 0? Pro) + Fo = F, es) (io + Pio) + Pro = 9; Pe 
cay | 1 ay 1 0. ‘ 
taf —pig)—ft=e, mali) | 


von denen die beiden ersten natiirlich identisch erfiillt sind, wahrend die beiden letzten tiberein- 
stimmend ergeben 


Ds9 = — % Ps0 « (86) 
Von den Gleichungen der Ordnung f geben wir nur noch die beiden an, aus denen x, und 03, be- 
rechnet werden kénnen: 


: La 
— & Bz, + 9 & P39 (1 — a)? >t Fo + 9X Pao Vin + 3 a? Hy) = 0, 


a (1—a) (87) 


& Pa + 9 & Pzq (1 — a)? ty A oi Sots +398) =0-| 


Man ersieht unmittelbar, daB auch bei den iibrigen Koeffizienten keine Komplikationen eintreten 
und man die so entstehende Entwicklung, die fiir endliche @5, sicher konvergiert, wiederum wie 


im ersten Abschnitt bis zu Werten von Qy) © /p benutzen kann, und sie damit in den Konvergenz- 
bereich der Grundlésung hineinreicht. Wiederum sind dann beide Entwicklungen nur Umordnungen 
voneinander, und es ist prinzipiell wieder méglich, aus der Bedingung, daB die Gesamtamplitude 
stationar wird, auf die Konvergenzgrenze der Grundlésung in der Nahe von 1/9 (1 — x) zu schlieBen 


(vgl. Ziff. 3). 


II. Der Fall « = 4/5. 11. Das Verhalten der Liésung in der Nahe der Resonanz der 
Grundschwingung. Wie schon in der Einleitung bemerkt, verhalten sich die Eigenfrequenzen 
des Systems (62), wenn « = 4/5 ist, fiir kleine Schwingungen wie 1 : 3, so daB gegeniiber dem Nor- 
malfall, den wir in IT untersucht haben, neue Phanomene zu erwarten sind, die fiir einen Freiheits- 
grad keine Analogie aufweisen. Wir untersuchen in diesem Abschnitt das Verhalten der Lésung 
in der Umgebung der Resonanz der Grundschwingung, die zugleich der ersten Ultraharmonischen 
zur Oberschwingung entspricht. Wollen wir auch fernerhin mit endlichen Amplituden der auBeren 
Kraft arbeiten, so miissen wir unseren Grundansatz I (6) modifizieren; nur, wenn wir die Ampili- 
tude der auBeren Kraft proportional / setzen wiirden, kénnten wir diesen alten Ansatz aufrecht 
erhalten. Da aber in beiden Fallen der Rechenaufwand nicht wesentlich verschieden ist, behalten 
wir ein endliches F bei und setzen an 


16 5 : 
=o + Epis, Lee ig Pe eee 


@, = Apis, O, = BB2849,4+0,.f8+°°-, . (88) 
Ds = spp + Pai 4 Psp FG See UMN ES Micra ae ise te 


Die folgenden @;, @; bis hin zu i = 9 beginnen dann jeweils mit einem Term endlicher Ordnung, 
wahrend erst fir 1 = 11 die ®;, O; mindestens von der Ordnung f8 ausfallen. Dies entnimmt man 


den Gleichungen (6). 


Die Glieder der Ordnung /)1’8 liefern je zwei iibereinstimmende homogene Gleichungen fiir A, 
B baw. fitr 9, gq [vgl. I (17)] mit den Lésungen 


4 4 
B==~A, Ost ae P39 + (89) 
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Der Koeffizientenvergleich der Glieder der Ordnung f° liefert dann 


4 36 5 3 /4\8 > 

a a (5 (A® + A 9 + 2 A G5) = F, 
5 9 3 ‘ 

ny Vee silts ic *y Al =" ae (— AB -|- A? P30 — 2 A 7) = 0 5 


9 29 4 ; 
53, + 4 ps) = Lee rg 5 30 — 95 (A® + 6 A? M9 + 3 v3), 


30 


9 1 
3 Bs, + 4 va — 5 1 P30 : A (AO 4" P30 — 3 Yo) » 
125 16 3 /4\8 5 
Pe G4 (so ea 35760] + Bs + a (3) (A" p39 — A yi.) = 0, 90 
123 og | A 2 2 oo 
oe co" (P50 4 P50) + Ps0 4 a A? xq + A Po) = 9, 
245 16 3 /4\8 
a 9 (Pro + 95 v0) + Oy — Tie ( 4 A Pin = 9, 
245 3 ‘ 
9 (oo + Pr0) + Pro 4 Ay, =9, 


E 16.2? L(A 
= 45 (0 35 P20] + Ooo 4 r(3) Poo —9, 


1 
— 45 (Bg + Poo) + Po0 4 Peo = 0. 


Die ersten vier Gleichungen stimmen mit den Gleichungen I (19) iiberein, nur daB bei der ersten 
dieser Gleichungen jetzt auf der rechten Seite die Amplitude F der duBeren Kraft auftritt. 
Eliminiert man aus den beiden ersten Gleichungen @,, und setzt ~)5 = A x, so folgt 
5 Ae 20 F 


3 
eae a a 


ie == 


(01) 
Setzt man diesen Wert von x, in die dritte und vierte Gleichung ein, so erhalt man durch die gleiche 
SchluBweise wie in I fiir x die Gleichung dritten Grades 


55 3 
g(x) = 41 x8 — 92x? —Al x ae £ x Be oad +3=0, (92) 


die mit der Abkiirzung 
a 7 
i =) a cy C3) 


direkt in die Gleichung I (20) mit €, an Stelle von € iibergeht, so daB wir alle in I fiir diese Gleichung 
gewonnenen numerischen Resultate unmittelbar auf den jetzigen Fall iibertragen kénnen. Die 
numerischen Rechungen wollen wir hier nur fir den Fall € = 0 durchfiihren, méchten aber darauf 
hinweisen, daB® die gleichen Anomalien in ahnlicher Form natiirlich auch fiir solche Werte von a 


. . . . 3 ae . 
auftreten, die in einem Streifen von der GréBenordnung Vp um « = 4/5 herum liegen. 


Unser Ziel ist zuniichst, einen Uberblick iiber die Abhangigkeit des Parameters A von der Fre- 
quenz Q zu gewinnen, also die A-Responsekurve in der Umgebung der Grundfrequenz aufzustellen. 
Diese A-Responsekurve ist natiirlich noch nicht die iibliche Responsekurve, da ja die Amplituden 
der Glieder sin t und sin 3 t hier von gleicher Gré®enordnung sind, und die maximale Amplitude 
aus diesen beiden Amplituden zusammengesetzt werden muB, um die wirkliche Responsekurve zu 
erhalten. 


Dazu berechnen wir zunachst die Skelette der A-Responsekurve, die fiir € = F = 0 aus der 
Gleichung 


f(x) = 41 23 —9 x —41x4+43=0 (94) 
folgen. Die Wurzeln dieser Gleichung werden 


xt = + 0,072 400, «x, = + 1,081551, x, = — 0,934442 . (95) 
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Setzen wir zur Abkiirzung 


1 2 
ANG ner m0 (82 «2 —9x-+ 41), (96) 
so wird : 
A(x}) = 0,33 982 , ZA (a;,) =0l00 009% A(x;) = 1,00842. (97) 
Tragen wir diese Werte in (91) ein, so erhalten wir wegen 
ONE =e i 98 
(3) =ltx=9 + 1B CY) 


fiir diese Skelette drei Hyperbeln, zu denen auch gleich noch die Gleichungen ihrer Asymptoten 


angegeben sind: 


oars 9 /Q\? ae ae ee Q 
fur x; : 3 (°.) — 0,61 167 Vp -A2A=)1, A pte — LiLo ae 
a Q 
faeeat ae 1,90780 fp -42=1, Afi = 0,97134— , (99) 
3 5 la 0 
oe ee VB. A? = A B'8 = 0,99582 
fiir x*: ah 1,81516 yp -A2 =1, p ; a 
F 5 : 
Fiir groBe Betrage von A, also kleine Werte von | Ae | erhalten wir mit 
iySy Lit 
eee Ee 100 
Ni eas (100) 


fiir die Wurzeln der Gleichung (92) 
g(x) =f(x) —Ax = (101) 


in erster Naherung 


Ae ai ree re at f fe Pe ; (102) 
da die Werte von f’(x*) sehr groB werden: 
f'(x;) = — 41,658 , ji (%,) == 034k f ,) = = 83.2). (103) 
Weiter wird 
Ae) =AGH eas aos, (104) 


und hieraus 
F {125 x* A’(x*) 25 
EA (Gt) (1) 


4, = A? A(af) 4 


Schatzt man die eckige Klammer ab, so zeigt sich, daB sie fiir x} negativ, fiir x; und x; hingegen 
positiv ausfallt. Also verlauft fiir groBe positive A die A-Responsekurve bei dem zu aj gehérigen 
Skelett links von der Hyperbel, fiir die beiden anderen Skelette rechts davon; fiir groBe negative 4 
ist es gerade umgekehrt. Damit ist der Verlauf fiir groBe Betrage von A bekannt. Die beiden 
positiven Wurzeln der Gleichung (92) fiir € = 0 fallen in eine Doppelwurzel zusammen, wenn [vel. 
I, (26)] 

wy) fit F 


c= "at ee 0,13 224 A2, et 0,73467 (106) 
ist. Fiir F'/A® > — 0,735 besitzt (92) drei reelle Wurzeln, fiir F/A? < — 0,735 nur noch eine, die 
negativ ausfallt. Demnach besitzt fiir F/A® > — 0,735 die A-Responsekurve drei Aste, unterhalb 
dieser Schranke hingegen nur noch einen. Umgekehrt: fiir positive A gibt es drei Aste, fiir negative 
A und |A| > |A,,| ebenfalls drei, fiir |A| <|A,,| hingegen nur noch einen. In der (A, Q)- 
Ebene tritt daher, wie wir zeigen wollen, fiir ein bestimmtes negatives A eine Verzweigung auf, die 
aber fiir die Lésung selbst keine Verzweigung darstellt, da zu diesem Werte von A verschiedene x 
und somit verschiedene 9 gehéren. Soll ein solcher Verzweigungspunkt existieren, so miissen sich 
die zu x, und x, gehérenden Zweige schneiden, also zu einem Tripel A, F, , fir x, und x, gleichzeitig 
die Gleichung (91) erfiillt sein. Daher miiBte gelten 


82 «2 —9 x, = 82 a5 — 9&4, also ty + 4 == : (107) 
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Da nach (92) 
9 


ae As = ky Al? (108) 
so mul} x, = 9/82 sein. Soll aber x, = 9/82 eine Wurzel von (92) sein, so muB gelten 
ioe R, a! Ald Ol F 
Se ae gg 0009» 4s = — 9,839854. (109) 


Ks ergibt sich also in der Tat ein negatives 4. Wir kénnen aus diesem Resultat zugleich schlieBen, 
daf sich die drei Zweige der A-Responsekurve in der oberen Halbebene nicht schneiden. 
Fiir diesen Wert von F/.A® wird (92) 


Al «3 — 9 x* — 26,839431 %4+3=0 (110) 
mit den Wurzeln 
) 
%1 = 95 = + 0,109756, XK = + 0,873217, Xx; = — 0,763461 . (111) 
Man sieht leicht, daB der Wert von x,, der zu «, gehdrt, kleiner ist als der Wert von ,, der zu x, 
und 2, gehért, also dem Verzweigungspunkt zugeordnet ist. 
Es bleibt noch der Verlauf der Responsekurve fiir kleine Betrage von A, also groBe Betrage von 
A zu diskutieren. Fiir A < 0 hat hier (92) nur eine Wurzel, die in erster Naherung gegeben ist durch 


eee oa) 9 AS 
Son (nari P20 ~~ 195 FF ° BES) 
so daB x, in der gleichen Annaherung gegeben ist durch 
: F 25 
ny~ wa 72 (113) 


mit zunehmender Frequenz klingen sowohl A wie @3. rasch ab. Ist hingegen A > 0, so gibt es drei 
Wurzeln, eine kleine positive Wurzel x,, die bis aufs Vorzeichen mit (112) iibereinstimmt, und zwei 
betragsmaBig sehr groBe Wurzeln, wobei x, positiv, x, negativ ausfallt. Die Wurzel x, liefert das 
gleiche Ergebnis, wie wir es eben im Falle A < 0 festgestellt haben: da sich in (113) ebenfalls das 
Vorzeichen umkehrt, wird x, jetzt zwar negativ, aber mit wachsendem Betrag von x, klingen A 
und @ 39 ebenfalls rasch ab. 


Fir die beiden betragsmaBig groben Wurzeln erhalt man durch Reihenentwicklung mit 


o=1/ya+41 


1 9 81 
XK» = — lol 
2 Go Val | 82° 8-41 4 lo" 
l 9 81 (114) 
ae - | —g + [o*]. 
< o jai | 82 9-41 yAl t Lo! 


Zu der positiven Wurzel x, gehért nach (91) 
AES [74 9 1 fo 25 
4 | 4s BS, gira : 
= 120 @ is eval ze ia) A° 2 t lol 
Fiir die negative Wurzel x, hat man nur das Vorzeichen von 41 umzukehren. In erster Naherung 
erhalt man daraus fiir x, und x, gleichlautend 
pees ea 25 F 
sae (56 Py ee 


Hy (115) 


(116) 


so daB also mit unbeschrankt abnehmenden A schlieBlich x, belicbig zunimmt. Da dabei aber auch 
3) unbeschrankt zunimmt, so ergibt dieser Ast der A-Responsekurve offenbar asymptotisch die 
reine Oberschwingung, stellt also die Ultraharmonische dar. 

Die Aste der A-Responsekurve, die zu den Wurzeln x, und x, in der oberen Halbebene gehéren, 
miissen nach dem letzten Ergebnis also ein Minimum von ~, als Funktion von A aufweisen, das wir 
naherungsweise aus (115) bestimmen kénnen. Eine einfache, aber langere Rechnung zeigt, daB 
dies Minimum eintritt, wenn gilt 


ter i25. 9 OSE aoe wre W7 
i+ 123 $ 
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Diese Gleichung gilt fiir x,; fiir x, ergibt sich die entsprechende Gleichung durch Vorzeichenumkehr 
bei der Quadratwurzel. Die Wurzeln dieser Gleichungen werden naherungsweise 
firm: C=a7a64, fir: C= a6: (118) 
Die Ergebnisse aller dieser Rechnungen sind nun in der Abb. 3 als die A-Responsekurve auf- 
getragen. Auf der horizontalen Achse ist (Q/@ )6-1/* aufgetragen, so daB die Resonanzstelle der 
Oberschwingung sehr weit rechts zu denken ist; die Abbildung ist im iibrigen rein schematisch, 
also die drei Skelette, die punktiert eingetragen sind, 
nicht maBstabgerecht gemaB den Gleichungen (99), 
sondern so eingezeichnet, daB alles Wesentliche der 
Kurven klar zu iibersehen ist. Da im allgemeinen 
P39 = Ax von der gleichen GréBenordnung wie A aus- 
fallt, zeigt auch die Responsekurve fiir die Gesamt- 
amplitude im allgemeinen einen gleichen Verlauf, mit 
einer Ausnahme: auf den Asten der A-Responsekurve, 
die sich der Q-Achse nach rechts hin und von obenher 
asymptotisch nahern, wachst ~3) umgekehrt proportional 


4A 


zu | an, wird also sehr groB. Diese Aste entsprechen 
also der Ultraharmonischen, worauf oben schon hin- 
gewlesen wurde. 

Man erkennt aus der Abbildung, daB es rechts von 
der Eigenfrequenz der Grundschwingung einen Bereich 
gibt, in dem neun verschiedene mégliche Lésungen 
unseres Systems (62) existieren, die der gleichen auBeren 
Kraft zugeordnet sind. 

12. Das Verhalten der Lésung in der Nahe 

Licata ida cee der Resonanz der Oberschwingun g. Da hier die 

Breqicas ier Ccunauniw dae, erste Subharmonische zur Grundschwingung in Reso- 

nanz zur Oberschwingung steht, wird, entsprechend 

zu den Verhaltnissen in der Naihe der Grundfrequenz zu erwarten sein, daB dem normalen Ver- 

lauf der Responsekurve in der Nahe einer Eigenfrequenz sich eine Subharmonische mit drei Asten 

iiberlagern wird. Ganz entsprechend dem Ansatz in Ziff. 4 machen wir jetzt fiir die Frequenz den 
Ansatz 


(C) =90 4H). (119) 


so daB die auBere Kraft gegeben ist durch 
Posi ouve (120) 


Fiir ®; und 0; bleiben die Ansitze (88) erhalten; daher bleibt (89) unverandert und in den Glei- 
chungen (90) ist nur bei der ersten auf der rechten Seite statt F Null zu schreiben, wahrend bei der 
dritten Gleichung rechts das Glied — 5/4 F hinzukommt. Dies neue System laBt jetzt die Lésung 
A = 0 zu, woraufhin alle ®; und O;, deren Index nicht durch drei teilbar ist, verschwinden. Da 
1 (P39 die Rolle des freien Parameters tibernimmt, folgt aus der dritten und vierten Gleichung (90) 
sotort 


Pa a) pene Go, 
abe 1D ae, elke 120 P30 (121) 
und damit fiir die Frequenz in erster Naherung 
Oey Varah Cee ta hide BE ters 
& Sag VB (3 Dee “AO vn} 2 (122) 


/ 


Das ist aber natiirlich nichts anderes als das erste Glied der Grundlésung, die ja ebenfalls in un- 
serem Ansatz enthalten sein muB. Fiir die Subharmonischen muB demnach A + 0 beibehalten 
werden. In diesem Falle ergibt die Elimination von #,, aus den beiden ersten Gleichungen (90) 
an Stelle von (91) jetzt, wenn wir gleich € = 0 setzen, : 

~ 120 


My (82 22 —9x + 41) = A? A(z). (123) 
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Einsetzen dieses Wertes in die dritte und vierte Gleichung ergibt jetzt mit Gay — A & fiir 4 die 
Gleichung 


> 2 ; 
a(x) = 41 2° —9 2 — 41 x + 3—- 5, = f(x) A =0. (124) 


Diese Gleichung unterscheidet sich von der Gleichung (92) nur dadurch, dab A jetzt nicht mehr mit 
x multipliziert ist. Aber hierdurch wird die Lésungsmannigfaltigkeit grundlegend geandert. Es 
gibt jetzt zwei Werte von A, fiir die (124) Doppelwurzeln besitzt: 


A, =—16,194760, xp = + 0655139, x, = — 1,090766,) 


Ao ae -- 16,130515 > X19 = 0,508 798 5 Xs — oF 237108 afi (128) 


In dem Bereich A, <4 < A, hat (124) jetzt drei reelle Wurzeln, fiir 2 < A, nur eine, und zwar nega- 
tive Wurzel, fiir 2 >A, ebenfalls nur eine, jetzt aber positive Wurzel. Fiir groBe Betrage von A, 
also fiir kleine |/ | gibt es also drei Subharmonische mit den dazugehérigen Wurzeln von (124) 

A 
FF)’ 
wobei die Werte von x7 und f’(«}) aus (95) und (103) zu entnehmen sind. Die zugehdrigen Skelette 
gentigen den Gleichungen 


= ay ct (126) 


OE ene ae eg ; 
& =5 +978 A A(st). (127) 
also, zusammen mit den Gleichungen der Asymptoten 
rene cee Lye ON 3 );— ; aye) 
pen & —0,61167 4278 =1, Apr = 0.571825, 
¢ - il 2 — 
fir xt: : (°.) —1,90780 22/6 =1, Api = 0,32378 ~ , (128) 
D \@o Wo 
pnts 1 /2\? oe es Ss ; Q 
fiir xt: : | —1,81516 A*Y¥B =1, A fr = 0,33194— 
Fir kleine Betrage von A wird gemaB (123) 
ne 
ty = A* A(x;) — A® A(x}) AeA A) (129) 


Co ae 
Das Vorzeichen des Koeffizienten von ) wird fiir x* positiv, fiir x* und x* negativ, so daf, falls 

: ‘ 2P 1 3g neg : 
A > 0, die A-Responsekurve fiir x, und x, links, fiir x, rechts von der Skelettkurve verlauft. Fiir 
A < 0 ist es gerade umgekehrt. 

Auch hier gibt es in der (A, Q)-Ebene fiir ein gewisses negatives A eine Verzweigung, die aber 
genau wie im vorigen Abschnitt fiir die Gesamtlésung keine Verzweigung bedeutet, da die zu diesem 
A gehérenden Werte von verschieden voneinander sind. Man erhalt 

§ 30 


Ay = — 1,554209 , Coa 5 = + 0,109756, +, = + 1,062379, x, =—0,952623 (130) 
mit den zugehérigen Werten von x, 
a Al ays an Per! A? ; 
fir x;: t= F594"? fiir x55: ¢, == 123,98781 730 (131) 


Fir kleine Betrage von A, also groBe Betrage von A gibt es jeweils nur eine Wurzel, die in erster 
Naherung gegeben ist durch 
A a / 3 
ARUN Ee Ba ES Py VE . (132) 
1M sag P30 V 


Diesem Zweig der Subharmonischen ist offenbar fiir 4 > 0 die groBe positive Wurzel von (124) zu- 
geordnet, er muB sich also mit zunehmendem A der zu x, gehérenden Skelettlinie nahern, und zwar 
wie wir oben gezeigt haben, von rechts her. 

Ist hingegen A < 0, so gehért dieser Zweig zu der groBen negativen Wurzel, muB sich also in 
der unteren Halbebene asymptotisch der zu x, gehérigen Skelettlinie nahern, und zwar ebenfalls 


von rechts her. 
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Im Grenzfall A — 0 wird nach (132) @g9 © jF und somit nach (121) x, ~ 2/3 F?/3. Den gleichen 


3 T . oT) . 
Wert fiir x, erhalt man auch aus der Grundlésung, wenn man dort Py ~ jF eintragt. Die Ver- 


zweigung der Grundlésung erfolgt also fiir x, ~ 2/3 F?”. 
In Abb. 4 ist nun die A-Responsekurve in derselben Weise wie im vorigen Abschnitt schematisch 
aufgetragen. Die auch im Normalfall auftretende Subharmonische entspricht dem Ast, der die 


A 


Abb. 4. Die A-Response-Kurve in der Umgebung der Frequenz der Oberschwingung. 


(-Achse schneidet; die beiden jeweils in endlichem Abstand von der -Achse beginnenden Aste 
sind die Subharmonischen zu den beiden anderen Eigenschwingungen, die zur Grundfrequenz ge- 
héren. Da in der Responsekurve der Gesamtamplitude zu diesen Subharmonischen auch noch die 
drei Aste der Grundlésung hinzukommen, gibt es also hier rechts von dem Verzweigungspunkt in 
der ‘Tat, wie in der Einleitung schon ausgesprochen wurde, insgesamt neun mégliche Lésungen. 


(Eingegangen am 24. Mai 1957.) 
Anschrift des Verfassers: Prof. Dr. C. Schmieden, Darmstadt, Technische Hochschule, Mathematisches Institut. 
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Uber die Beulsicherheit von Kreisplatten mit kreiszylindrischer Aelotropie 


Von S. Woinowsky-Krieger 


1. Differentialgleichung der Biegung einer aelotropen Platte. Das Material der Kreisplatte be- 
sitze eine kreiszylindrische Aelotropie, deren Symmetrieachse mit der geometrischen Symmetrie- 
achse der Platte zusammenfallen mége. Alle Ebenen, die speziell diese Achse enthalten sind dann 
in bezug auf das elastische Verhalten des Plattenmaterials gleichwertig. Die folgenden Ausfiih- 
rungen seien auf achsensymmetrische Belastungs- und Deformationszustande der Platte beschrankt; 
auch soll von allen Vereinfachungen der Kirchoffschen Plattentheorie, die Isotropie des Materials 
ausgenommen, Gebrauch gemacht werden. 

Zwischen den an einem Plattenelement wirkenden Normalspannungen o,, o, (Abb. 1) und den 
zugehérigen Dehnungen ¢,, ¢, bestehen nun die Beziehungen 

Or oe Che Ur 


ey Sa fe Tle Tegan Fe (1) 


worin EF, FE, und E’ drei elastische Konstanten bedeuten. Die Dehnungen sind ihrerseits durch die 
Beziehungen 


mit der radialen Verschiebung u verkniipft, und wenn man 
auch die Plattendurchbiegung w einfiihrt, so hat man noch 


: : 0 
die Relation 
dw 
Le ear A (3) 
Die Biegemomente der Platte 
h/2 h/2 
mp [ o.2 dz), Me [O28 az, 


—h/2 —h/2 
Bn So 
@> >, 


Abb, 1. Plattenelement im Koordinatensystem. 


wo h die als konstant vorausgesetzte Plattendicke bezeichnen 
soll, lassen sich jetzt sukzessive durch (¢,, ¢,) und alsdann 
durch w ausdriicken. Sie ergeben sich schlieSlich in der Form 


d?w vy, dw 1 dw d?w 
m, >= D, (a i: > Me = D,( ean le a) (4) 
worin ist 
E,h® Ei, h? 
ee en en EG 9 9) ©) 
: E 
p= a und Y= ae ; (6) 


Die Bedingung des Gleichgewichtes der an einem Plattenelement (Abb. 1) angreifenden Krafte- 
paare liefert noch die radiale Scherkraft 


[ dw 1 d’w 1 dw 

ge 1? (ae as) Ds ae (7) 

Projiziert man schlieBlich die am selben Element einschlieBlich der Flachenlast p r dr dd wirkenden 
Krafte auf die Achse z, so erhalt man die Gleichgewichtsgleichung 


d 1 
(+ 7)% +P =o, 


die nach Substitution von (7) auf die folgende Differentialgleichung! der Biegung der achsensym- 
metrisch belasteten aelotropen Platte fihrt: 


d*w 2 d*w 1 dw 1 dw Pp 
ee de pel ae eS = = 
k ( r? dr | 43 i) De (8) 


dr* r dr® | 
1 Vgl. auch S.G. Lechnickij, Anisotrope Platten, S. 207. Moskau 1947 (Russisch). 
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Die Aelotropie des Materials ist hierbei durch den Parameter 
_ he (9) 
gekennzeichnet. 


2. Beuldruck einer Platte aus aelotropem Material. In Abwesenheit einer transversalen Be- 
lastung hangt die Belastungsfunktion p(r) in der Gleichung (8) nur von dem Kraftespiel in der Mit- 
telebene der Platte und von den Kriimmungen dieser Ebene nach der vorausgesetzten Deforma- 
tion ab. 

Die Plattenscheibe sei nun einem gleichmaBigen Radialdruck von der GréBe n je Langeneinheit 
des Umfanges ausgesetzt (Abb. 2). Das Gleichgewicht eines Scheibenelementes driickt sich bei 
achsensymmetrischem Zustand der Scheibe durch die Be- 
ziehung 


se eral (10) 


dr ' r 


zwischen den radialen Kraften n, und den tangentialen 
Kraften n, aus. Dan, = ho, undn,=ho,, so ergibt sich 
aus den Gleichungen (1) und (2) 


du are ny Ue 1 /n Ny 

ae Dah ( — e) pee aber —F). _ 
Die Radialverschiebung u la4Bt sich aus (11) sofort eli- 
minieren, alsdann auch n, mit Hilfe der Gleichung (10). Fiir 


n, besteht noch die Grenzbedingung (n,),_. = —n. Dies 
Abb, 2. Platte unter gleichférmigem Radialdruck. alles fiihrt zu dem Endergebnis 


r \k—1 r \k—1 na/l k\/r\k 
n, ee) 5 n= —nk(7) = u ="(p—=) (7) (12) 


mit dem schon friiher eingefiithrten Parameter 


i, 
Ez (13) 


Fir k = 1 liefern die Ausdriicke (12) den Spannungszustand n, =n, = —vn einer isotropen 
Scheibe; bei k > 1 ergibt sich ein unbegrenztes Anwachsen der Spannungen gegen den Mittelpunkt 
r = 0, womit die Anwendung der Lésung, die ja eine volle Kreisscheibe voraussetzt, sich in jenem 
Bereich von selbst ausschlieBt. 

Unter der Annahme einer achsensymmetrischen Beulform der Platte ergeben nun die Druck- 
krafte (12) eine transversale Belastung von der GréBet 


k= 


pee liad dw\ _ [/r \F—l die k [rv \k—2 dw 
Py dr (r As 7) ie (7) ees (7) dr |* ag) 
Geht man hiermit in die Gleichungen (8) ein und setzt noch zur Abkiirzung 
na : r 
D, Aa = =0; (15) 
so erhalt man die Differentialgleichung fiir die Auslenkung w in der Form 
2 5 hk? i 
TVnsihg care 0 Paya Reine re k—2 = 
wh +2 + (Re a + (ke + 7)ui = 0, (16) 
wobei Striche bzw. rémische Ziffern Ableitungen nach 9 bedeuten. Fiihrt man schlieBlich 
yp =— w'/a = —dw/dr als eine neue abhangige Variable ein, so geht die vorstehende Gleichung 
in die Gleichung dritter Ordnung 
a 1 ne 2 pee 
Ferd deer ele eal] oe 


ee C. B. Biezeno und R. Grammel, Technische Dynamik, Bd. Sn Oo2ee Aut: Berlin/Gottingen/Heidelberg 
Jo2. 
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iiber. Wir werden vor allen Dingen mit der Differentialgleichung 


er’ +ey + (Pokt!— ko =0 (18) 


zu tun haben, deren Lisungen auch Lésungen der Gleichung (17) sind. Mit 


Rey eae 
ere (19) 
geht diese Gleichung in 
ap dy 4k? 
2 a“ 6 
whe ne dx ! [= wr pase, 20) 
d. h. in eine Besselsche Differentialgleichung iiber, so daB sich fiir y eine Lésung 
Ecce 
a ioe ries AG (21) 


ergibt, worin C, eine im Folgenden unbestimmt bleibende Konstante bezeichnet. Setzt man noch 
zur Abkiirzung 2 k/(k + 1) =, so erhalt man eine zweite Lésung der Gleichung (20) entweder in 
der Form C, Y,,(x) oder in der Form C, J_,,(x), je nachdem ob ys ganzzahlig (dies nur wenn k = 0, 
1, oo) oder nicht-ganzzahlig ist. In beiden Fallen fiihrt die zweite Lésung zu unendlichen Deforma- 
tionen im Mittelpunkt der Platte und 
ist daher auszuschlieBen. Eine dritte 
Lésung der Gleichung (17), endlich, 
hat die Form einer Potenzreihe 
P= EQ + ay ott? 
“a, At or +3 +--+), 

deren erstes Glied, dem Ausdruck (7) 
gemaB, eine gegen die Plattenmitte 
unbegrenzt anwachsende Querkraft 


Abb. 3, Beuldruck in Abhangigkeit von D_/a. Abb. 4, Eine durch Rippen versteifte Kreisplatte. 


ergibt. Fir eine Platte ohne Bohrung und ohne einen Stiitzpunkt in der Mitte missen wir 
wiederum C,; = 0 setzen. Die Tangentenneigung gy(g) der achsensymmetrisch verbeulten Platte 
ist also endgiiltig durch die Gleichung (21) bestimmt. Die zugehérige Beulkraft mége unter den 
folgenden beiden Annahmen berechnet werden: 

a) Die Platte ist am Rande eingespannt. Die Randbedingung ist y = 0 fir 9 = 1; 
bezeichnet man also mit j,, die erste Nullstelle der Funktion J,,(x), wo 4 = 2 k/(k + 1), so bestimmt 
sich der zugehdrige Eigenwert A aus der Gleichung (19) zu A = (& +- 1)j,,/2, und die erste der Glei- 
chungen (15) liefert den kritischen Druck 


(kk +1)? jn Dy Y 
my, =o aa (22) 
Praktisch geht man natiirlich yon u-Werten aus, fiir die die Nullstellen der Funktion J (x) tabu- 
liert sind und bestimmt danach den jeweiligen Wert k = u/(2 — yz). Die Relation (22) ergibt somit 
eine Kurve, die in der Abb. 3 zur Darstellung gebracht ist. 
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b) Die Platte ist freigestiitzt. Die Randbedingung folgt aus den Formeln (4) und schreibt || 


sich | 
YP Pp 23 
e: re ae Res ( 


Unter Beriicksichtigung der Differentialformel 


Fix) = J,—12) == Sa) 


x 


und der Bezeichnung (19) 1aBt sich der Gleichung (23) die Form i 
[x J, —1(%) — (u —%) Ju(Je=1 = 9 (an 
geben. Ist i, = 0 die kleinste Wurzel dieser Gleichung, so ergibt sich fiir den kritischen Radial- ] 
druck der Ausdruck \] 
k +1) in Dy | 
= eee a: (25) |] 

Die zugehérige Kurve findet man ebenfalls in der Abb. 3. 


3. Beuldruck einer gleichmaBig gespannten Platte. Wie es vor allem in den Gleichungen (9) und |/ 
(12) zum Ausdruck kommt, wurde bisher eine aelotrope Platte behandelt, die aus homogenem Ma- 
terial besteht und eine tiberall gleiche Dicke hat. Es liegt nun nahe, in ahnlicher, diesmal aber nur 
angeniherten Weise, die Stabilitat einer durch radiale und konzentrische Rippen, etwa nach 
Abb. 4 ausgesteiften Platte zu untersuchen, die an sich aus einem isotropen Material besteht. Sind 
die Rippenabstande einer solchen Platte klein im Verhaltnis zu ihren Gesamtbemessungen so kommt 
ihr elastisches Verhalten demjenigen einer aelotropen Platte etwa gleich. 


Wir fiihren jetzt an Stelle eines einzigen, durch die Gleichung (9) definierten Parameters deren 

zwei, namlich 
aa aa und ks ; é (26) 

Hierin sind I, und I, die auf die Einheitsbreite der Platte bezogenen Tragheitsmomente der Platte 
in tangentialer baw. radialer Richtung, Ff, und F’, die betreffenden Querschnittsinhalte. Wahrend 
nun die im vorigen Abschnitt berechneten Eigenwerte fiir den Sonderfall k’ = k galten, sollen sie 
im Folgenden fiir den Sonderfall k’ = 1 ermittelt werden. Liegt also der tatsachliche Wert von k’ 
zwischen k und 1, so kann der zugehérige Beuldruck naherungsweise durch lineare Interpolation 
zwischen den Beulwerten, die diesen beiden Grenzwerten von k’ entsprechen, gewonnen werden. 

Im Grenzfall von k’ = 1 verhalt sich die Platte hinsichtlich der in ihrer Ebene wirkenden Krafte 
wie eine isotrope Scheibe, und ein gleichformiger Umfangsdruck von der Intensitat n ergibt einfach 
n, =n, = —n. In der Differentialgleichung (8) ist somit 


een ee nie 1 a 


dr? | +r dr 


anzunehmen. Mit den Abkitrzungen (15) erhalt man schlieBlich, an Stelle von (16), die Differen- 


tialgleichung 
2 Kk He Ke? 
apLY. ul 1 
wit mee | ( 2 Z| wit 4 ie ze Or (27) 
Fiihrt man wieder mp = — w!/a als eine neue Variable ein, so lat sich die vorstehende Gleichung 


in die Form 


DOR e ate 2) 
(io | 2) | Es + (# 7) 9° (28) 


umschreiben, wonach sich fiir p abermals die Besselsche Differentialgleichung 


d 
a + % bags (CE *) oi : (29) 


mit x= /¢@ ergibt. Die an der Stelle x = 0 singulare Lésung dieser Gleichung und ebenso eine dritte 
noch vorhandene Lisung der Gleichung (28) bleiben mit derselben Begriindung wie in dem vorigen 
Abschnitt ausgeschlossen. Die endgiiltige Lésung ist also 


yp =CJ,{A @) ; (30) 
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_ sofern C eine Konstante bezeichnet. Die Berechnung der Eigenwerte soll 


. wiederum in zwei Fallen 
durchgefiihrt werden: 


a) Die Platte ist am Rande eingespannt. Die Randbedi so =0 bei o—1 fiihrt mit 
Riicksicht auf (15) auf die Beulwerte Ree ck ee 


ei 
Nkr = Jk op > (31) 


_ wobei j, die erste Nullstelle der Funktion J,(«) bezeichnet, D,~ EI, und E der Elastizitatsmodul 
des Materials ist, wahrend, mit Riicksicht auf die Rippen, v +0 gesetzt werden darf. 


Tun a’ 
0 
60|— 


40 


20 


Abb. 5. Beuldruck in Abhangigkeit von Di Ja, k’ = 1. 


b) Die Platte ist freigestiitzt. Die Randbedingung (23) fiihrt diesmal zu der Beulgleichung 


Pie EO 0: (32) 
Ist 1 =i, die kleinste Wurzel dieser Gleichung, so ergibt sich fiir die Beulkraft der Ausdruck 
m, = ib. (33) 


Numerische Ergebnisse, auf die die Gleichungen (31) und (33) fiihren, sind in Abb.5 durch Kurven 
dargestellt. 


4, SchluBbemerkung. Es mége noch bemerkt werden, daB die beiden Differentialgleichungen (16) 
und (27) auch zur Lésung des Problems einer kreiszylindrisch aelotropen Platte herangezogen 
werden kénnen, die auBer dem Radialdruck (oder auch Zug) in ihrer Ebene noch einer Transversal- 
belastung ausgesetzt ist. Wenn q die Intensitaét dieser Belastung bezeichnet, so mu der rechte 
Teil der homogenen Gleichungen (16) und (27) durch ein Glied q/D, erginzt werden und ihre 
Lésungen durch ein diesem Gliede entsprechendes partikulares Integral. Ware schlieBlich anstatt 
eines Radialdruckes n ein gleich groBer Radialzug vorgeschrieben, so hatte man diesem Umstande 
durch eine Abanderung des Parameters /? in — /? Rechnung zu tragen. 


(Eingegangen am 7. Juni 1957) 


Anschrift des Verfassers: Prof. Dr.-Ing. S. Woinowsky-Krieger, Quebec (Canada), Université Laval, 
Boulevard de l’Entente 
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Warmespannungen in Hohlzylindern mit temperatur- 
abhangigen Stoffwerten 


Von R. Trostel 


1. Allgemeines. Mit dem vorliegenden Aufsatz soll ein Beitrag zum Thema Warmespannungen 
im Bereich héherer Temperaturen geliefert werden, wo man von der Temperaturabhangigkeit der 
Materialwerte nicht mehr absehen kann. Eine allgemeine geschlossene Lésung dieses Problem- 
kreises gelingt nur fiir den Spezialfall elastisch-inkompressiblen Materials (y = 1/2), der nachfolgend 
behandelt wird. Er erlaubt eine verhaltnismaBig einfache und iibersichtliche explizite Darstellung 
der Spannungen durch Integrale, die sich gegebenenfalls, wie die unter Ziff. 4 behandelten Beispiele 
zeigen, geschlossen, im allgemeinen aber sehr bequem numerisch bewaltigen lassen. Auch der fiir 
die Beriicksichtigung der Temperaturabhangigkeit der Materialwerte charakteristische Sach- 
verhalt, da die aus einer rein mechanischen Belastung herriihrenden Spannungen von der Ver- 
teilung eines gleichzeitig vorhandenen ungleichmaBigen Temperaturfeldes abhangig sind, laBt sich 
gut erkennen. 

In Ziff.2 wird unter der Voraussetzung elastisch-inkompressiblen Materials der Spannungs- 
zustand fiir ein beliebiges ebenes achsensymmetrisches Temperaturfeld #(r) und fiir einen beliebig 
mit } veranderlichen Elastizitatsmodul E(#) und Warmeausdehnungskoeffizienten «,(?) angegeben. 
In Ziff.3 folgt die Bestimmung einer allgemeinen Lésung fiir das Temperaturfeld unter Beriick- 
sichtigung der Temperaturabhangigkeit der Warmeleitzahl. Die Berechnung eines Warmespan- 
nungsproblems geht dann, wie in den unter Ziff.4 abgehandelten Beispielen gezeigt wird, folgen- 
dermafBen vor sich. 

Nach Ermittlung des Temperaturfeldes ? = J(r) sind bei Vorgabe der Temperaturabhangigkeiten 
E=E(9) baw. «,=o,(9) die Funktionen E(9(r)) = E(r) sowie «,((r)) =a,(r) gegeben. Sie kénnen 
nunmehr, zusammen mit 7(r) in die fiir die Spannungen unter Ziff.2 gefundenen expliziten Integral- 
lésungen eingesetzt werden, womit nach Integration dann die Spannungen folgen. In den in Ziff. 4 
behandelten Beispielen wurde der Warmespannungszustand in einem Stahlrohr untersucht, wobei 
jedoch auf den Umstand, daB hier die Querdehnung nicht 1/2 ist, keine Riicksicht genommen 
werden konnte. Trotzdem erlauben diese Beispiele im Vergleich ihrer Ergebnisse mit denjenigen, 
die sich nach der itblichen Theorie unveranderlicher Materialwerte unter Verwendung von y = 1/2 
ergeben, eine Aussage iiber den EinfluB8 der Temperaturabhangigkeit der Materialwerte 
generell zu machen. In diesem Sinne kommt man an Hand eines Vergleiches der in Abb.5 dar- 
gestellten Kurven zu dem Resultat, das die Beriicksichtigung der Temperaturabhangigkeit der 
Werkstoffwerte bei Stahlrohren auf gréBere Warmespannungen fihrt, als sie sich unter Ver- 
wendung unveranderlicher Materialwerte ergeben. Die Unterschiede zwischen den Spannungen sind 
jedoch gering, wenn bei dem zugehérigen Temperaturfeld die Randtemperaturen vorgegeben 
sind. In diesem Falle ergibt namlich die genauere Berechnung des Temperaturfeldes gegen- 
iiber einer Rechnung mit konstanter Warmeleitzahl einen nur kaum merkbaren Unterschied, 
so da die Verdnderungen der Spannungsbilder praktisch allein von der Veranderlichkeit von 
& und a, herriihren. Diese sind jedoch nicht erheblich. Die Verhaltnisse kénnen sich jedoch 
erheblich andern, wenn fiir das Temperaturfeld eine andere Randwertaufgabe vorgelegt ist (z.B. 
Warmeiibergang am AuBenrande), so dafs die genauere Temperaturverteilung von der mit kon- 
stanter Warmeleitzahl bestimmten bereits merklich abweicht. Hier zeigt sich an Hand der Abb.6, 
da man mit der iblichen Theorie unverdnderlicher Materialwerte schon merklich kleinere 
Warmespannungen ermittelt, und daB man nur noch dann mit den Ergebnissen der iiblichen 
Theorie unveranderlicher Werte E und «, auskommt, wenn man wenigstens dem tatsadchlichen 
Temperaturverlauf Rechnung tragt. 


2. Ermittlung des Spannungsfeldes. Wir betrachten den mittleren Teil eines sehr langen Hohl- 
zylinders, in dem ein ebenes rotationssymmetrisches Temperaturfeld herrsche und an dessen Mantel- 
flachen r=r, baw. r=r, die Driicke p, baw. p, wirken mégen (Abb.1). Die Rotationssymmetrie 
der Belastung und des Temperaturfeldes bedingt dann das Entstehen eines achsensymmetrischen 
Spannungszustandes, der dariiber hinaus im mittleren Rohrbereich von der Zylinderachsenordinate z 
unabhangig sein mu8. Samtliche SpannungsgréBen, und damit im Sinne des Hookeschen Gesetzes 
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auch samtliche Verzerrungskomponenten, kénnen daher nur von r abhangig sein. Insbesondere 
erhalt man also fiir die Axialdehnungen 


e, = = fo) (1a) 
und nach Integration 
w =w(r,z) =zf(r) + g(r).~ (1b) 


Nun ist jedoch im mittleren Bereich eines sehr langen Rohres kein Rohrquerschnitt z = konst. 
vor einem anderen ausgezeichnet. Daher miissen diese Rohrquerschnitte wihrend der Deformation 


eben bleiben (Abb.2), womit 
d(t) = Konet, =A,,.. g(r) = konat. = A, (Le) 

verbunden ist und aus (1b) 
w = w(z) = A,2z+ A, (1d) 


hervorgeht. Die Konstanten A, und A, kénnen aus den Lagerungsbedingungen des Zylinders an 
seinen Enden bestimmt werden. 

Nach dieser vorbereitenden Betrachtung gehen wir zur Berechnung des Spannungszustandes 
tiber, der im mittleren Rohrbereich und wegen der Rotationssymmetrie des Problems allein durch 


Abb. 2. 


Radialspannungen o, =¢,(r) Ringspannungen o, =o,(r) und Achsialspannungen o, =¢,(r) charak- 

terisierbar ist. Zu ihrer Ermittlung stehen die Gleichgewichtsbedingung fiir das Kérperelement 

in Radialrichtung zur Verfiigung 
d 

a 1 9,)—% —ai) oder? = (6, —-o,) 5 (2) 


sewie die durch die Warmeausdehnungsglieder verallgemeinerten Hookeschen Gesetze, die mit 
y = 1/2 die Gestalt 


dv 1 1 a) 

bee = F@® |°— 2 (o, +¢,) ert a,(T) dT, (3a) 
v 1 1 A(r) 

Ey = - == E() Oy 5) (o, 0,) | se, aiihe (3b) 
dw 1 1 al 

6, = = At = Bq | — x (+ o0)| + oT) AT (3c) 


haben, wenn mit 9(r) die Temperatur gegeniiber dem spannungslosen Ausgangszustand (dessen 
iiberall gleiche Temperatur wir mit Null annehmen wollen) bezeichnet wird. 
Zunachst erhalt man nach Addition von (3a) bis (3¢) 
d 
ete+4=et+—+4= 


r 


ar) 


§ (ee) +A =f a(t at 
=O 


woraus nach Integration mit der Integrationskonstanten C, 


Weare é 
emoey =F] r(f ovat) age (4) 


T= 


136 R. Trostel: Warmespannungen in Hohlzylindern mit temperaturabhangigen Stoffwerten — Ingenieur-Archiv 


hervorgeht. Beachtet man nun noch, da aus (3c) 


4 
0, =F (6, +59) +A, E(B) — E(B) f oT) aT (Ga) 


folgt, so gewinnt man schlieBlich aus (3a) oder (3b) mit (4) und (5a) 


r? r? 


(7) (7) 

4 E E dd 4 E 

gg =a | a{T) dT — raul | os vn) Pas [ro pdr Ce 

T=0 E PSN 

Hiermit und mit der Gleichgewichtsbedingung (2) folgert man somit fiir o, die Differentialgleichung 
do, E/4 2. 1,a\ 
als afr (0) ar). 

deren Lésung mit der Integrationskonstanten C, 


Tea i E E ese: | rou) ir) dr (6) 


lautet. Die Konstanten C, und C, werden aus den Randbedingungen 


o(r=r)=—p, und of(r=r,)=—p, mu 


CG, a 9 mae _- ms (7a) 
E E 
73 dr pdr 


E dp 4 Ea, 
Ce [[os( [ea ar) ee Ge] [ea (7b) 


bestimmt, so daB letztlich fiir die Radialspannungen 


a E "E do 
fe [fe fife 
r; E / db un Us 
CG, == p; +. (P; — P,) moo + 2 fa([raa dr ie, Ta oe ¥ = (8a) 
[ae i [pe [ial | tea &) & 


und damit aus (5b) fiir die Ringspannungen 


E dd Cie Wis ¥; 
v= 6,2 | Pog dr + 3 Cp — P+ (pe) ra sas 
E 
Jae 
ee ie emer 3 do B d 
: at [ae a [mem Ge at | ([ moi a) & 
LEH ie dd " UF y 
+2/a(| ron Ge Oar eT 
r E E dd 
t fae [ al [ ema aa 


hervorgeht!. Sowohl in (8a) als auch in (8b) ist die Konstante A, nicht mehr enthalten. Auch bei 


1 Fiir eine numerische Integration sind die in den Doppelintegralen auftretenden inneren Integrale jeweils als be- 
stimmte Integrale in den Grenzen r;...r aufzufassen, 
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dieser Betrachtung zeigt sich also, daB die Lagerungsverhaltnisse des Zylinders an seinen Enden 
auf die Radial- und Ringspannungen keinen Einflu8 haben, sondern sich nur auf die Achsialspan- 
nungen o, auswirken. Fiir diese folgt namlich aus (5a) mit (8a) und (8b) 


gels oe | Z dr 
rs 7: 
1 Or) (r) ” 
Oa 5 (c, Oy) Ef od tay -- EA, SS EA, ae se i a,(T)dT — p,; + (P:— Pa) = ars i ae 
=e '"—() a 
E 
lie dr 
Bet Ef, 40 "E 
Ta gt [ae 5) | rau + [al | rr? of —— ar) dr 
E a rs 
+ 2 | 2 (| 7 v0 Ge | a ae : ANNES 
Tr. aay 
; r3 dr | a({ dp a dt 


und man kann nun zwei Fille behandeln: 


1) Das Rohr ist zwischen zwei starren Wanden gema Abb.3 eingeklemmt. In diesem Falle 
liegt ein ebener Verzerrungszustand mit 


ow 
Of = =) = 0 (8d) 
vor. 
(Eel aes Sree Se 
er ae ae 2) 
Abb. 3. Abb. 4. 


2) Die Ausdehnung der Zylinderenden werde in Achsialrichtung durch eine, auf die Zylinder- 
enden wirkende resultierende Kraft P gemaB Abb.4 behindert. In diesem Falle miissen die in 
einem Rohrquerschnitt z = konst. wirkenden Achsialspannungen die Resultierende P haben. Dem- 
nach folgt A, hierfiir aus der Gleichung 


r 
a 


1B 
[rad =ae 


Ts 
U 


Nach Einsetzen von (8c) folgt hieraus mit Beachtung von 


Ta Ta "a Ta Ta 
| (6, +09) dr = | ro,dr + | rode = | ro,dr + r 4 (ro,) dr 
"t % % ay z: 
"a "a 
=| ro,dr +180,|"*— | rodr= port — pat 
o — 


schlieBlich 
Ta Hr) 
li ts Sa | ee ie . - 
A, == —(- +$@.2-p.) + rb | x(T) dT \ dr \. (8e) 


T=0 
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3. Das Temperaturfeld. Unter Verwendung des Fourierschen Warmeleitungstheorems 


(q = —AV®#) liefert die Warmemengenbilanz am Kérperelement im stationaren Falle, sofern keine 
Warmequellen vorhanden sind, 
Vq=—VAVI)=0, also VAVI)=VAVI+AAD=0 


bzw. im vorliegenden rotationssymmetrischen Falle 


ee 1 ‘| 1 d (ra) =0 


dr dr ' dr? ' r dr r dr 


woraus 


ri = G, ibaa 19) a9 =“ ar 


und nach Trennung der Veranderlichen und Integration 
{4(d) d? = C,Inr + C, (9) 


folgt. Bei Vorgabe der Temperaturabhangigkeit 2=/(9) der Warmeleitzahl ist Gleichung (9) 
integrierbar, und die Konstanten werden aus den Randbedingungen des Problems bestimmt. 


4.Beispiele. a) Wir betrachten ein Stahlrohr mit freien Enden, dessen Mantelflachen 
r =r, = 12cm baw. r =r, = 18cm die Temperaturen ?, = 350° C bzw. 0, = 280° C gegeniiber 
der Temperatur 3 =0°C des spannungslosen Ausgangszustandes besitzen sollen. Die Temperatur- 
abhangigkeit der Materialkonstanten kann im Bereich bis etwa 400°C nach dem ,,Werkstoff- 
Handbuch fiir Stahl und Eisen“, 3.Aufl., Diisseldorf 1953 mit geniigender Genauigkeit durch 


NG) eae (10a) 
XP) = O% +% 10, (10b) 
E(0) = E, — E,# (10c) 


approximiert werden, wobei z.B. fiir einen unlegierten Stahl an Hand des in diesem Buche ge- 
gebenen Kurvenmaterials 


Ay =ACO = OSC) = 0,12 [cal em sec + °C], A, = 0,000.07 [cal em sec-1 °C FJ - 
Kg OO 0 WO) 2 Oe Oe |. Opp = 0,001 10 Ge]. 
Ey = EO = 0° CG) = 2.1,- 10* [keem-2]; B= 3,2 (keen *~ C4] 
abzulesen sind. 
Zunachst folgt mit (10a) aus (9) nach Integration 


A i Asian 
Ayi—F#=C,nr+C, baw. Ar) =a V(r) i, (C,Inr + C,). 


/ 


Die Integrationskonstanten ergeben sich aus den Randbedingungen OC =T)=o, und r=) 
==, 20 


Caen) 
¥\ Co [saad ewan i 
C, = fact aD Sek Ra ea EN EY gras altel 


4 4 jin £2 4 _ ia ie 
9 ese Oy r 0 v 
ee ee a 0 
hervorgeht. 7 


Selbstverstandlich erhalt man fiir /,>0, d.h. fitr 29/A,> 00 aus (11) die Lésung fiir das Tem- 
peraturfeld bei konstanter Warmeleitzahl, indem man zunachst (11) mit %, =A,/A, baw. x, =x, =A,/A, 
in der Form 


A In — In — 
B(r 1 — Bi, %_— 9.) =P — |/ (4 — 9)? — + (% — 0)? 
4 In In“ 


schreibt, und nunmehr fiir x, 0; baw. x, 0, das Temperaturfeld #(r, x,—0;, x,—@,) nach 
MaBgabe der Taylorschen Formel in der Form 

oo a 

Blrs 1% — 95 — 9) BO ep) — | — I 


Ty %y9) Xp 


approximiert. Man erhalt wegen 


oe an ate ee 
/ lige In — 


| A 
P(r, my) =P | 


sowie 


od | iP od 


sth ; eek —' 
OX4\r, x, X, 1 To OXo|r, x1, Xp in Ta 
nee fe 


Ti % 


schlieBlich fiir das Temperaturfeld bei konstanter Warmeleitzahl die bekannte Relation 


Oy (r) = —— eee (12) 


Ein Vergleich von (11) und (12) fiir die vorliegenden Zahlenwerte zeigte, daB die genauere Be- 
rechnung des Temperaturfeldes bei vorgegebenen Randtemperaturen kaum nennenswerte 
Abweichungen gegeniiber der Naherungslésung @,(r) fiir konstante Warmeleitzahl A) ergibt. Die 
Ergebnisse kénnen jedoch recht verschieden sein, wenn man nicht die Randtemperaturen, sondern 
z. B. eine Aussage iiber die durchflieBende Warmemenge vorgibt (siehe hierzu das nachste Beispiel b). 

Mit Riicksicht auf die geringen Unterschiede zwischen (11) und (12) benutzen wir fiir das 
folgende die Naherungslésung #,(r) fiir das Temperaturfeld, womit die Méglichkeit einer geschlossenen 
Integration der in (8) auftretenden Integrale gegeben ist. Nachdem man (12) sowie (10a) bis (10c) 
in (8a) und (8b) eingesetzt hat, gewinnt man mit den Abkiirzungen 


o Ba 
r [eee 1— 
Se ea he 
aa = [dar Le) = Ey 1; ce wee Ar] 8 ae 2 
re rr ’ ri < 
Oa Ba \ 
Er, 4 = E, 82 /99\2 ott Oi oO 
Pils | ee sala us 2In '@ iL Bp (3) aN Mt N Oa mae 
Ti ri 
4 ; = l 
Dir} = | al 72 a ar) dr =| as ®,(r) dr = (13) 
o Oa Ta 
fee ee —— , m= 4 
Bi Oy Di Bi Ee ne Pa 
= — Ey o% 8: —~ og a oe ab, _ oa : a: ) | 
ie J 0 a 
tj iat DO; 


fiir die Radial- und Ringspannungen mit p; = p, = 9 


{| Pa(r P(r) { 
o,(r) = 2 ®3(r,) Eres ~ B,(ra))’ (14a) 


E 
pt Pl) or) + Or) 


: leas 14b 
D,(ra) P3(ra) ( ) 


) 


og(r) = 2 P3l(r,) 
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wihrend man fiir die Achsialspannungen nach (8c) 


iL 
Sat Br) 5 P(r) + Or) 


D,(ra) P3(ra) 
erhalt. Die Konstante A, ist im Falle des ebenen Verzerrungszustandes Null und folgt fir den 
Fall freier Zylinderenden, wenn die resultierende Achsialkraft P sowie Innen- und AuBendruck 


p; und p, Null sind 


(14) 


O(r) 
o({r) =—E { o(T)dT + 4, E + 2®,(r,) 
T=0 


lg d(r) E EF, % 
E| f a(T) dT) dr eal : = 
we Ee a BS tee? Ey Js 2 Ey Xo e. 
jee ee (144) 
JrEd oo Ee 
( 1 a ; 
wobei fiir die Integrale J, mit Jy = f r+ dr = > (r2—ri) die Rekursionsformel 


gilt. . - . 
Aus (14) erhalt man mit E,=«,,;=0, also fir E=E , %= 0% die aus der iiblichen 
Theorie mittels unveranderlicher Materialwerte fiir » = 1/2 folgenden Relationen 


1 (") In Ui 
i: Tr; 
ot) = — Ey 049 (81 — 8.) —az — ra | (15a) 
t=) 98 
Ta : 
Tne 
Wet (7) 1+ Les ie 
a(t) = — Eye (8: —8,)| ——G —*), (15b) 
ale In” 
1 att 
o,(r) =—2 Eyo%49(0:— 9,) (15¢) 


Die Auswertung der sieben Gleichungen (14) und (15) fiir die Zahlenwerte des gewahlten Beispieles 
ergaben die in Abb.5 dargestellten Kurven. Dabei wurden die nach der iiblichen Theorie gemaB 
(15a) bis (15c) ermittelbaren Spannungen gestrichelt dargestellt. Die ausgezogenen Linien geben 
den genaueren Spannungsverlauf nach den Gln. (14) an, und man erkennt, da die tatsachlichen 
Warmespannungen sogar gréBer sind als die in itiblicher Weise berechneten Werte. Die Unter- 
schiede sind jedoch nur geringfiigig, so daB man sich also im allgemeinen, sofern die Rand- 
temperaturen vorgegeben sind, mit der tiblichen Theorie, auch im Bereich héherer Temperaturen, 
begniigen kénnen wird. 

b) Als weiteres Beispiel betrachten wir dasselbe Rohr wie unter a) mit derselben Innentem- 
peratur J; = 350°C, jedoch soll nunmehr am AuBenrande Warmeibergang in Luft der Temperatur 
J, = 20°C stattfinden, die das Rohr quer mit einer Geschwindigkeit w = 20 m/sec anstrémt. 
Die Randbedingung am AuSenrande lautet dann unter Heranziehung des Newtonschen Abkiihlungs- 
theorems 

q{r =r,) =—A i = —« [0,—Hr,)], also O(r,) + se 


=" a dr |r = 1, 


an Oy ? (16) 


wobei H(r,) =, die Randtemperatur des Rohres und « die sog. Warmeiibergangszahl bedeuten. 
Diese folgt z.B. nach Merkel zu 


«= 0,092 |= ob vue.” 17) 
’ ag dat aaa? ( 
Es ergibt sich, wenn in erster Naherung 3, = 310°C gesetzt wird, 
OAD B alae [cal crise secet uGet |e 


* Hierin sind: Az, die Warmeleitzahl, CpL die spezifische Warmemenge und yz das spezifische Gewicht 
der Luft fiir die mittlere Temperatur Oy = (8q + Oz)/2. 
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A 


| 
| 


od 


; v =UH0 
ae 7; =7¢cm. t= 20°C 
H,=280°C m=18em 7000 
Abb. 5. Warmespannungen zu Beispiel a: Abb. 6. Warmespannungen zu Beispiel b: —.—.— tbliche Theorie 
——— genauere Theorie; unter Verwendung des genauen Temperaturfeldes; ——— genaue 
— — — — ibliche Theorie. Theorie; —-— — — iibliche Theorie. 
Mit A?) =A, —A,0, erhalt man damit aus (16) 
pi esse a) ea 
o a dr\r = "o iy 


woraus mit (11) schlieBlich die Relation 


Ag 9, 2 oy & ra [Ag FOUN Ag q\" 4 g & Ta Ag 
(2 9.) + 2 i, aan cs (2 9, (? i] - 2 i, re ne 9, 


mit der Lésung 


Parnes Vl ala a yee alae é 9) | (; 6) = 310,5° C (18a) 


He SO 8 


—- Ag 
ue i" Uae a, Ai 
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hervorgeht. Bei Zugrundelegung einer konstanten Warmeleitzahl A =A, erhalt man dagegen aus (16) 


Ao raln is 
8, = 9, = 0: ———_ + 9, ———_ 7 = 3178". (18b) 
fei 2 role 
a3 x ri a 


Hier liegt also in den Randtemperaturen eine Differenz von 7°C vor, die bereits erhebliche 
Unterschiede in der GréBe der Warmespannungen nach sich zieht. Da die Warmespannungen nach 
der iiblichen Theorie unveranderlicher Materialwerte bekanntlich proportional der Temperatur- 
differenz A) =0;—%, sind, ein Umstand, der auch bei der genaueren Rechnung noch annahernd 
erhalten bleibt, bedeuten die 7°C Unterschied in der AuBentemperatur bereits eine um 


Ad Di baad Da —— 
I Ue Mise Nhe Saeki 1,22 fache 


Steigerung der Warmespannungen nach der iiblichen Theorie. Diese Steigerung erhéht sich 
noch etwas, da die Beriicksichtigung der Veranderlichkeit von E und o, (siehe Beispiel a) eine 
VergréRerung der Warmespannungen gegeniiber der Theorie unverdnderlicher Materialwerte 
ergibt. Diesen Sachverhalt zeigen die Kurven von Abb.6, worin die gestrichelten Linien die 
Warmespannungen nach der iiblichen Theorie unter Heranziechung des sich fiir konstantes A =A, 
ergebenden Temperaturfeldes 
ln + BayIn — 
dn(r) = —__— (19a) 
igo 
ri 
bedeuten. Die strichpunktierten Linien stellen die Warmespannungen nach der tblichen Theorie 
unveranderlicher Werte E und «, unter Verwendung des genauen Temperaturfeldes dar, das wieder 
im Sinne des Beispieles a) nahezu exakt durch 
8; In 2 + Ba, In = 
,(r) = ——___— (19b) 
c a 


In — 


Y% 


approximierbar ist. Die ausgezogenen Linien sind die Warmespannungen die sich bei Be- 
riicksichtigung der Veranderlichkeit von EF und «, mit dem genaueren Temperaturverlauf (19b) 
ergeben. Man erkennt nun sehr gut, daB die gréBeren Abweichungen durch die Abweichungen im 
Temperaturfeld bedingt sind. Benutzt man also in der iiblichen Theorie konstanter Werte E und 
x, das genauere, die Veranderlichkeit der Warmeleitzahl enthaltende Temperaturfeld, so erhalt man 
bei Stahlrohren auch im Bereich héherer Temperaturen eine noch zufriedenstellende Uber- 
einstimmung der errechneten Warmespannungen mit den tatsachlichen Gegebenheiten. In einer 
weiteren Arbeit wird demnichst iiber eine allgemeine Behandlung stationdrer Warmespannungs- 
probleme bei beliebiger Temperaturabhingigheit von FE, x, und A sowie fiir beliebige, jedoch kon- 
stante Querkontraktionswerte » berichtet. 


(Eingegangen am 14. Juni 1957.) 


Anschrift des Verfassers: Doz. Dr.-Ing. Rudolf Trostel, Technische Universitat Berlin, 
Lehrstuhl f. Mechanik (Prof. Dr. Ing. J. Szabo) Berlin-Charlottenburg 2, Hardenbergstr. 34. 
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The Flexure of Infinite Rectangular Plates of Varying Thickness 
By H. D. Conway 


1. Introduction. The following investigation concerns the bending of a rectangular cantilever 
plate of varying thickness, clamped along a longer edge and bent by any system of loading applied 
to a parallel but otherwise free edge. The plate is considered infinite in length. 

The first investigation of a similar cantilever plate problem but one having constant thickness 
seems to have been first made by MacGregor! who treated the case of a single normal concentrated 
load applied at a point on a longer side remote from the shorter sides of the plate. MacGregor 
considered that the solution would shed light on the analysis of the flexure of helical gear teeth. 
His theoretical deflections were checked experimentally with excellent agreement. Other theoretical 
analyses of similar problems have been made by Bitiner?, Jaramillo? and Koitert. Finite difference 
methods have been applied to the problems by Holl. 


2. Analysis. We consider here the cantilever plate of varying thickness, a sketch of the plate 
being shown in figure 1. Making the usual assumptions of plate bending theory, it is not difficult 
to show® that the bending of a plate having a thickness variation which is a function of y only 
is governed by the equation 


eD oa @D hhh z 
DAAw + 2 = PAG eee = — Aw ly yp = 
dy oy ey” 
2D Pw 
= 1 
(L—») FO oS =0 (1) 
with the usual notation. To obtain a solution to /// 


//) 
///) 


the equation, a very convenient choice for the ///// 
variation with y of the flexural rigidity D is 


D =D, e” (2) 


where Dy and 6 are constants. The factor e’Y then Abb. 1. Infinite Plate of Varying Thickness. 
appears in each term of equation (1) and can be 
cancelled throughout. By selecting suitable values for Dy and 6, itis thought that the variations 
in equation (2) can be frequently made to conform quite well to those encountered in practice, 
for example, in gear teeth. : 

Substituting the function (2) in equation (1), it follows that 


ow 


AAw 4 2b 2 Aw + b?Aw — b (1 —») = 0. (3) 


Ox? 


Consider, for example, the periodic loading along the free edge whose intensity, Q, per unit length 
is given by the equation 


Q = Pcosyx. (4) 
It will be shown that the deflection of the plate can be written in the form 
w = Fly) cosy x. (5) 


where F(y) is a function only of y. Substitution of the value (5) in equation (3) then gives the 
differential equation for F(y) as follows 


FV 4+ 26F™ + (62? ~2y?) F’— 26°F’ + (yt— 7 By’) F=0. (6) 
For solutions of the form F = A e+’, the characteristic equation is 
74 +2623 + (b? — 2 y?) A? — 2 by? + (yt — v By”) = 0 (7) 


1 C, W. MacGregor, Mechanical Engineering 57 (1935) p. 225. : 

2 E. Bittner, Momententafeln und EinfluBflachen fiir kreuzweise bewehrte Eisenbetonplatten, Wien 1938. 
3 T. J. Jaramillo, J. Appl. Mech., Trans. ASME 72 (1950) p. 67. ; . 

4 W.T. Koiter and J. B. Alblas, Proc. Kon. Ned. Akad. v. Wetens., 57 B (1954) pp. 250, 259, 549. 

5 D. L. Holl, J. Appl. Mech. 4 (1937) 8. 

6 S, Timoshenko, Theory of Plates and Shells, p. 195, New York 1940. 
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and the problem is then one of extracting the four roots of this quartic. We write equation (7) 
in the form 

(A + & + Pi) (A + o% — By) (A + O%2 + Ba) (A + & — By) = 9 (8) 
which, on multiplication gives 

AS + 2 13 (0% + Og) +A? (of + 0% — Bi — BE + 4% 9) + 

FA (2 of dy + 2 of of — 2 Oty Bt — 2 0% 3) + (9) 
+ (og — Bi) (0 — Be) = 0. 
Comparing the coefficients of 2° in equations (7) and (9) we note that 


04 +O, = 6b (10) 


and equation (9) becomes 
yale, oe ae os ig 2 4X, — BP? — B2) +A (2 a 1, b — 2 x, BF — 2 x, 63) + | (11) 
+ (of — Bi) (03 — B2) = 9- 
Comparing next the coefficients of A? and / in equation (11) with the corresponding ones in equation 
(7) we see that 
Oy = One DIZ (12) 
and thus equation (11) becomes 


M4 2bP +2 (H+ —H—P) + aly 0 + PD| + (GA) (F—B)= 0-8) 


Finally comparing equations (7) and (13), we see that they are identical provided 
é : b2 : 
PEP = yay 


(14) 
ae 


From equation (14) we obtain 


B= Zt storys, 
aa ita e 
and these, together with equation (12) ee the roots of equation (7) which are then 
ee: oo Veta +4 byyy, 
le =— 5 — VP +4 + 4by/r, 
Ig =—¥ + yO +4 —4b y/o, io 
(yo elit + 4y2—Aby ly. 


bo 


For y=0.25, which is an appropriate value for several materials, including steel, we obtain 


s set 2 ee +42 + 2yb> | 

janie Carter 

1 =— 3 +5) Pp apadyD » ee 
beer 5 t/Pp4paaye- 


The deflection of the plate is then given by 


wo =(Aedy + Beby + Coby + Deus) cos yx (18) 


a 
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where A, B, C and D are found from the boundary conditions on the two longer sides, which are 


(ala) 
a= Oe By 
a a 19) 
ae oe (Pa Ow ( 
=€4 | — ba Seipe he =) ; 
a ? oy? ie Ox2 0, Dye {ays | y| Ox? = P cos VX. 


It will be seen that the solution given in equation (18) is also valid for a plate bent by a corres- 
_ ponding distribution of moments applied to the free edge of the plate. Further, for any periodic 
distribution of normal forces and moments which can be expanded in a Fourier cosine series, the 
deflected form is readily obtained in the form of a Fourier cosine series. Similar solutions in Fourier 
sine series are obtainable and the general problem of periodic loading is formally solved. Of course, 
this solution applies only to the region remote from the shorter sides of the plate, since we are 
ignoring the boundary conditions on these sides. 


3. Numerical Example. As a numerical example, we write y = 7/2 and consequently the loading 
intensity along the free edge is Q@ = P cos (z x/2 a). The thickness constant b is chosen as — 2.07944/a, 
thus making the plate thickness at the loaded edge half of that at the clamped edge. 

Assuming vy = 0.25, we find from equations (17) that 


d,a = — 0.34420, J,a = 2.42364, J,a = — 1.23658, A,a = 3.31602. (20) 


The constants in equation (18) determined by boundary conditions in equation (19) are then 
found to be 


A =— 0.13903 Pa’/D,, B = 0.05488 Pa?/Dy, C = 0.10102 Pa*/D,, D =— 0.01687 Pa?/D,. (21) 
The deflection along the loaded edge is finally 
Pa? TX 
Wy —a = 0.08536 Dy cos 5 (22) 
4. Conclusion. It will be seen that the deflection and moment M, are both zero on edges 
x=-+a, +3a,... Hence the above solution is also valid for a plate of finite length simply 
supported on any two of these edges. 
So far as the variable thickness version of MacGregor’s single concentrated load case is con- 
cerned, the loading is written in the Fourier integral form 


oO 


12 ‘ 
Oren ere [ cosy d(y a) (23) 
0 


and the deflections are given by the Fourier integral 
iS al | (A, ehy + B,ekyY + C, e¥ + Dye) cos yxd(ya). (24) 
a 
0 


The evaluation of this integral must be performed numerically after the value of the thickness 
constant b has been selected to best fit the actual variation. This may be done by computing 
the integrand over a certain range for various values of ya and evaluating the integral in this 
range by Simpson’s rule. , 

Beyond this range, the integrand may be written in simpler form and the integral evaluated 
in closed form. This involves considerable labor and, since the calculations must be repeated for 
each value of b selected, they will not be made here. Numerical calculations of Fourier integrals 
such as the above have been made by Girkmann’ and the author®. 


(Eingegangen am 25. Juni 1957.) 
Anschrift des Verfassers: Professor H. D.-Conway, Thurston Hall, Cornell University, Ithaca (New York) USA. 


K. Girkmann, Ing.-Arch. 13 (1943) p. 273. 
H. D. Conway, J. Appl. Mech. 22 (1955) p. 353. 
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Mathematische Behandlung gesteuerter Abkiihl- und Anwaérmvorgange 


Von H. Jackel 


1. Einleitung. Die vorliegende Arbeit beschaftigt sich mit der Lésung des allgemeinsten 
gesteuerten Warmeleitproblems in festen Kérpern, die beziiglich der Warmeleitfahigkeit isotrop und 
homogen sind. Dieses Problem ist dadurch gekennzeichnet, daB verschiedene Teiloberflachen eines 
Kérpers, dessen Temperaturverteilung zu bestimmen ist, von Medien umgeben werden, deren 
Temperaturen irgendwelche Funktionen des Ortes und der Zeit sind. Hinzu kommt noch, daB im 
Kérper selbst innere Warme entwickelt werden kann. 

Die Untersuchung derartiger Probleme findet man erstmals in einer Abhandlung von J. M. C. 
Duhamel!. Ohne jedoch die dort in drei langeren Abschnitten ausgefiihrten Gedankengange zu 
benutzen, soll im Folgenden ein anderer Weg beschritten werden, der durch eine Arbeit von 
A. Kneschke? nahegelegt wird. In dieser ist unter Verwendung einer EKinfluBfunktion ein Theorem 
zur Ermittlung der Temperaturverteilung eines Kérpers mit konstanter Anfangstemperatur und 
einer lediglich von der Zeit abhangigen Temperaturfunktion des die gesamte Kérperoberflache um- 
gebenden Mediums abgeleitet. Dieses Theorem, insbesondere der Begriff der EinfluBfunktion, ist 
nun einer solchen Erweiterung fahig, da®B die Lésung der allgemeinsten Aufgabe gelingt. 


Abgesehen von der einfachen physikalischen Interpretation, die sowohl die EinfluBfunktion als 
auch der Lésungsansatz gestatten, sind die Ergebnisse Jeicht zu tiberschauen und fiir den Praktiker 
bequem zu handhaben. Dariiber hinaus liefert das Theorem nicht nur, wie etwa die bei H. S. 
Carslaw® angegebene allgemeine analytische Fassung der Duhamelschen Gedankengiange, die Lésung 
fiir das Kérperinnere, sondern sie gilt auch auf der Kérperoberflache. Dies ist fiir die Behandlung 
technischer Probleme, bei denen die Temperatur auf der Begrenzung des Kérpers eine Rolle spielt, 
von Bedeutung. SchlieBlich kénnen aus der allgemeinen Lésung ohne Schwierigkeit beliebige 
spezielle Ergebnisse abgeleitet werden. 


Zur Erlauterung der Theorie werden ein Warmeleitvorgang in einem Kreiszylinder und die 
Temperaturverteilung in einem Betonquader wahrend des chemischen Abbindeprozesses behandelt. 


2. Formulierung der Aufgabe. a) Probleme ohne innere Warmeentwicklung. Wir 
betrachten einen hinsichtlich der Warmeleitfahigkeit isotropen und homogenen Kérper K, dessen 
Gesamtoberflache F' sei. Seine Temperatur zur Zeit t=0 bezeichnen wir mit f(P). Dabei ist unter 
P irgendein Punkt von K zu verstehen. Die Gesamtoberflache F bestehe aus r Teilflachen F, 
(k = 1, 2,...,r). Die Temperatur ®, der Grenzschicht des die Teilflache F,, umgebenden Mediums 
andere sich mit der Zeit t und dem Ort P, also ®, = @,(P,t). Die ®,(P, t) werden Steuerfunk- 
tionen genannt, weil sie die gesuchte Temperatur T(P,t) des Kérpers wesentlich beeinflussen. 


Setzen wir Fouriersche Warmeleitung voraus, so muB T(P,t) in K der Differentialgleichung 


oT 

4 @ATH=DiT(P,»)} =0 (1) 
gentigen. Dabei bezeichnet A den auf die Ortskoordinaten anzuwendenden Laplaceschen Operator 
und a? die Temperaturleitzahl des Korpers. 


Die Randbedingungen kénnen, entsprechend der Problemstellung, verschiedener Art sein. Wir 
wahlen sie in der bekannten, relativ allgemeinen Form 
A oT 
& snarls r) =L{T(P,i}p, = O,(P, 1), (2) 
Fy 


a, On 


mit n als aufserer Normalenrichtung, 2 als Warmeleitzahl und «, als Warmetibergangszahlen. 


1 J. M. C. Duhamel, J. Ec. polytech., Paris, 14 (1833), S. 20 
® A. Kneschke, Ing.-Arch. 24 (1956) S. 77. 
8 H. S. Carslaw und J. C. Jaeger, Conductions of heat in solids, S.19, Oxford 1947, 
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Sind die Teiloberflachen F;, mit der Temperatur @,(P, t) isotherm (x), > co), so geht (2) iiber in 


(T)z, = PD, (P, t). (3) 
Im Falle adiathermaner Abgrenzung (a, = 0) ergibt sich 
Ol 
(a), tn o 


Die Anfangstemperaturverteilung T(P, 0) ist durch f(P) gegeben. Sie wird die vorgeschriebenen 
Randbedingungen im allgemeinen nicht erfiillen. Deshalb wollen wir von der gesuchten Tempera- 
turverteilung T(P, t) nur verlangen, daB T(P, 0) im Inneren von K mit f(P) iibereinstimmt, wahrend 
T(P, t) fiir jede Zeit t > 0 die vorgeschriebenen Randbedingungen erfiillen soll. 


Nunmehr kann die gesuchte Temperaturfunktion T(P,t) aus folgendem Randwertproblem 
bestimmt werden 


D{T(P,)} =0, T(P,0)=f(P), ; 
Dee = OAP.1) Shir S00 (= 1;9,5087) | ©) 


b) Probleme mit innerer Warmeentwicklung. Wir haben bisher angenommen, daB im 
Inneren des Kérpers keine Warme erzeugt wird. Diese Voraussetzung lassen wir jetzt fallen und 
bezeichnen mit W(P, t) die im Kérper in der Volumen- und Zeiteinheit erzeugte Warmemenge, 
mit T(P, t) seine Temperatur, mit ®,(P, t) die Steuerfunktionen und mit f(P) die Anfangstempera- 
turverteilung. Die Aufgabe lautet dann bekanntlich 


D{T(P,1)} = ft CaN) al 


ze = F(P) 
TAT(P a a= OAPt) fix te 0 (k= 1,2... 


6 
eh (6) 


wobei ¢ die spezifische Warme und @ die Dichte des Kérpers ist. Das inhomogene Problem (6) 
kann in bekannter Weise durch den Ansatz 


T(P,t) = T(P,t) + U(P,t) (7) 


auf das halbhomogene (5) zuriickgefiihrt werden, wobei U(P,t) ein partikulares Integral der in- 
homogenen Differentialgleichung ist. Es wird dann 


F(R) =P) — U(P,0) 
und 
®,(P,t) = ®,(P,t) — L {U(P,1)} p, . 


Gilt insbesondere ®,(P,t) = 0 in (6), so wird diesem zunachst ungesteuerten Problem durch (7) 
ein gesteuertes, halbhomogenes Problem der Form (5) zugeordnet. Es geniigt deshalb, wenn wir 
uns im Folgenden mit der Lésung der Randwertaufgabe (5) beschaftigen. 


3. Lésung des Randwertproblems. Zunachst setzen wir voraus, da von samtlichen vorkommen- 
den Funktionen die ersten und zweiten Ableitungen nach den Ortskoordinaten sowie die erste Ab- 
leitung nach der Zeit und dem spater auftretenden Parameter 1 fiir alle in Betracht kommenden 
Werte der unabhangigen Veranderlichen existieren, endlich und stetig sind. Diese an die einzelnen 
Funktionen gestellten Bedingungen kénnen zum Teil noch abgeschwacht werden, worauf wir hier 
jedoch nicht weiter eingehen wollen. 

Wegen der Linearitat der Differentialgleichung und der Randbedingungen in (5) ist fiir T(P,t) 
der Ansatz méglich 

DE) Pt) an bee) 


mit 
D{T(P,H)} =9,  To(P,9) aie (8) 
PATO. akan to) 0 
und 
Deh) 0 P= 0, 0) 
L{T(P,)}p,=O(P,t) fir t>0 (k=1,2,...,7). 
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Die Lésung des vollhomogenen Problems (8) ist bekannt!. Sie lautet 


oo) are 
T(P,t) = YA,ul(Pye”, 
o( Pt) a (10) 


Dabei sind die u,(P) die normierten Eigenfunktionen der Randwertaufgabe 


Au(P) + 6?u(P) =0, (11) 
L {u(P)} » = 0 


und die 63, die entsprechenden EKigenwerte. Auf (11) wird man durch den Bernoullischen Produkt- 
ancatz~ 1 ((P.1) =u Pye" -setubrt, 

Um T,(P,t) aus (9) zu bestimmen, betrachten wir zunachst folgendes Teilproblem fir 
WS IBS iP 


D {T,(P,#)} = 05; T,,(P,0) = 0 (12) 
L {Tul P. 0} = 0,, O.CP.t)) funy t 0 ian 2?) 
Dabei ist 6,,= 0 fiir i--k und 0,,— 1 fir i= k, Offenbar gilt dann 
T(Pi) = 3 Ty(P.) (13) 


Deshalb geniigt es, (12) fiir beliebiges k allgemein zu lésen. Zu diesem Zweck wird vorerst an- 
genommen, es sei diejenige Funktion Gf bekannt, die fiir jeden beliebigen Punkt P des Kérpers K 
die Temperatur wiedergibt, die unter dem EinfluB der auf die Teiloberflachen F, wirkenden AuBen- 
temperaturen 0;,0,(P,t) (i = 1, 2, ...,r) im Laufe der Zeit t—1>0 entstanden ist. Gj ist eine 
EinfluBfunktion, die in spezieller Form erstmals von A. Kneschke zur Lésung bestimmter Warme- 
leitprobleme verwendet wurde. Sie hangt vom Ort, von dem friiheren Zeitpunkt 7 und der Zeit- 
differenz t—t ab, also 


* = GUP, r,1—1). 


Nunmehr ist es méglich, einen physikalisch sehr naheliegenden Ansatz fiir T,,(P,t) anzugeben 
und aus diesem in Verbindung mit (12) Gg zu bestimmen. T,,(P,t) setzt sich namlich zusammen 
aus der Temperatur Gf(P,0,t), die von ®,(P,0) herrithrt und der Summe aller Temperatur- 
ainderungen 

t 
bese iia 
a= Gi (P, t,t — 7) dt, 


ot T=T 


die durch ®,(P, t) wahrend der Zeit von t= 0 bis r= t hervorgerufen wurden. In den angefiihrten 
zwei Termen ist jedoch bereits die Temperatur mit enthalten, die zur Zeit t= t, verursacht durch 
P,(P,t), erst wirksam wird. Demzufolge muB noch Gf(P,t,0) subtrahiert werden. Also ergibt 
sich schlieBlich 

t 


T,(P,t) = GE(P, 0,1) + | E Gt(P, %,t—7) 


0 


dz —G#(P, t, 0). (14) 


T=T 


Da T,,(P, t) Lésung von (12) sein soll, bietet sich uns unmittelbar eine Moglichkeit, Gg zu be- 
stimmen. 


Wir bilden zuerst mit dem Ansatz (14) 


— 5 G(P, t,0) + PACK P, 1,0) =0. 


1 R. Courant u. D. Hilbert, Methods of Mathematical Physics. New York und London 1953, Bd. 1. 
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Beachtet man, daB 


a) 
= =, Gi(P, t, 0) 


t=t 


(lex 64. t ole 


ist und aus D {Gi (P, t. t)} = 0 fiir 0 <= t < tauch D{Gi(P, 0, t)! = 0 und D {GX(P, 7, t—t)| =0 
folgt, so ist die Differentialgleichung offensichtlich erfiillt, wenn gilt 


AG.(P,t,0)=0, D{R(P, 7,1} =0. 
Weiter ergibt sich aus den Randbedingungen 


t 


a) i 
L {Ty(P,)}n, = L {C8(P, 0, )} », + | [GLCP EO}, dz 


0 
— L{GEP, t, 0) p, =6,,9,(P,2) (¢=1,2,...,7). 


Diese werden befriedigt durch die Forderung 


L {Gi(P, t, O)}r, ee DWP) Mawes Lal as,.F) 
12G(P et.) ty == 0. furs 20. 
SchlieBlich ersieht man aus (14) sofort, daB die Anfangsbedingung T,,(P,0) = 0 in K erfiillt ist. 
Damit ist G(P, t,t) und folglich auch Gg(P, 1, t— 1) bestimmt. Es kann aus nachstehenden 
zwei Randwertaufgaben, deren Lésungsmethoden bekannt sind, berechnet werden. 
Zuerst ermitteln wir aus 
AH (Pt) =0.; 
L {H,(P,)}r, =—b,P(P,) (6 =1,2,..5 1 () 
eine Funktion H,(P,t). Danach setzen wir 
GLUE. 7, 0) ee, (Poa) 
und bestimmen aus 
DACA Part) 02 8 Gp(F, 7,0) — 1, (P. 7), (16) 
LAGE, 10), — 0 far t> 0 


die gesuchte Funktion Gx(P, T, t). 


Die Aufgabe (15) kann als ein von einem Parameter t abhingiges stationares Problem gedeutet 
werden, wahrend (16) ein von dem Parameter 1 abhangiges vollhomogenes Warmeleitproblem 
darstellt. Entsprechend (10) erhalten wir fiir Gj(P, 7, t) 

GE(P, 7,1) = Y BP(r) u,(P) e—* et , 


w=1 


B®(z) = LHP. 1) u(P) do. 


Aus (17) ist ersichtlich, daB G(P, 7,0) gleich der Entwicklung von H,(P,7) nach den Eigen- 
lésungen des Problems (16) wird. Diese Entwicklung stimmt aber nur im Inneren von K mit 
H,(P,1) itberein. Auf der Kérperoberflache geniigt sie den Randbedingungen von (16) und nicht 
den in (15) geforderten. Wir schreiben deshalb unseren Ansatz (12) in der Form 


t 


T (2.0) = CrP; 0,8) -- | leg CRP. #12) dz — H,(P, t) (18) 
P f ek 
und wollen kiinftig mit H,(P,t) immer die Lésung des Problems (15) und mit Gf(P, 7, 0) die 
Entwicklung von H,(P,7) nach den Eigenlésungen von (16) bezeichnen. Unter Beriicksichtigung 
dieser letzten Bemerkung iiberzeugt man sich leicht, daB die durch (18) gegebene Temperaturfunk- 
tion T,,(P, t) tatsaichlich die Lésung des Problems (9) ist, insbesondere erfiillt sie die geforderten 
Randbedingungen, gilt also auch auf der Kérperoberflache. 
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Damit haben wir das allgemeine, gesteuerte Warmeleitproblem (5) auf folgende Teilaufgaben, 
deren Lésungsmethoden bekannt sind, zuriickgefiihrt : 


das vollhomogene Problem (8), 
die r stationéren Probleme (15) 
die r vollhomogenen Probleme (16), 

) 


die einfachen Integrationen (18). 


und T(P,t) hat die Form 


PB) ee sce ae Oj! 
= 
‘2 oe) a 
+ >) | >) BRO) u,(P) en et + | 


eS i é Bel 


 u,(P) e— °6u'—*) dr—H,(P,t)|. (19) 


Die Eindeutigkeit der Lésung (19) des Rapdwertproblems (5) leuchtet physikalisch unmittelbar 
ein, sie bedarf jedoch des mathematischen Beweises. Wir verweisen diesbeziiglich auf die 
Literatur 1. 

Durch partielle Integration kann die Lésung (19) auf folgende, bisweilen vorteilhaftere Form 
gebracht werden: 


25 E 282 (¢—7 
,u,(P) et — @ f BM(r) u,(P) 53 6 8H dr 


ae Ce Va 


+ a ay BY(t) u,(P) — H,(P,1)} . (20) 


TP ae An 


Da der letzte Term fiir alle inneren Punkte von K Null wird, ist (20) besonders gut zur Bestimmung 
der Temperatur dieser Punkte geeignet. 


Wir bemerken noch, daB aus (20) fir f(P)=0 und L{T(P,t)\,, = O(P,t) das Duhamelsche 
Integraltheorem in der bei H. S. Carslaw angegebenen analytischen Form resultiert. Unter den 
eben angefiihrten Bedingungen folgt namlich entsprechend (20) fiir die Temperaturverteilung 
T(P,t) im Inneren des Kérpers K 


t 
foe) 


t 
t aa: A . 
T(P,) =a | SB,(2) 0 (P) eae = [ZS Bye) u,(P) MO ae 
6 eee 


6 a 
Das ist aber das Duhamelsche Integraltheorem. Es gilt im Gegensatz zu der allgemeinen Lé- 


sung (19) nur im Inneren des Kérpers K. 


4. Spezielle gesteuerte Probleme. Aus der allgemeinen Lésung (19) kann nunmehr in Verbindung 
mit (8), (15), (16) und (18) durch geeignete Wahl der ®,(P, t) eine ganze Reihe spezielle Ergebnisse 


abgeleitet werden, von denen wir jetzt einige nennen wollen. 


a) Hat die Steuerfunktion die Form 


P,P, t) = 74(P) Plt) 


fiir die Teilflachen F,,, dann setzen wir 
Gi(P, t,t — t) = 9,(t) G,(P,t— 7). 
Damit ergibt sich aus (18) 
t 
Ty (P, t) = 9,(0) G,(P, ¢) + J gult) GP, t— 1) de — (0) HP) , (21) 
wobei g(t) = = oy(7) ist. 


1 H.S. Carslaw u. J. C. Jaeger, a. a. O. S. 24. 
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Die H,(P) und G,(P, t) (k = 1, 2,...,1r) sind bestimmt durch 


AH,(P) = 0, 

EAC P) = — On u(r) (4,2, <.. 57), (22) 
DiG(P.)} =0, G,(P,0) = H,(P), 

TA CP = OF far” ot > 0° 123) 


Die Aufgaben (22), (23) unterscheiden sich von (15), (16) dadurch, daB jetzt die Randwertprobleme 
nicht Aa parameterabhangig sind. Fiir T(P, t) folgt somit entsprechend (19) 


T(P, t) = Y Ay 1,(P) ee! 


w=1 
(oa) 2 2 v —a at 3 
‘os PO) 3 BYP u,(P) ee + f(x a BYP u,(P) "#9 dr — yt) H(P)}. (24) 
Gres L= b= 


Der Fall, daB die Steuerfunktion fiir die gesamte Oberflache F gleich y(P) p(t) ist, kann aus (21), 
(22), (23), (24) unmittelbar abgelesen werden. 


b) Sind die Steuerfunktionen nicht ortsabhangig, also 


D,(P, t) = 7, (t) » 


dann brauchen wir in a) nur y;,(P) = 1 zu setzen. Die allgemeine Lésung dieses Problems ist daher 
ebenfalls von der Form (24). Sind in a) die y,(t) = 1, so ergibt sich die bekannte Lésung fiir nur 
ortsabhangige Temperaturfunktionen der umgebenden Medien. 


c) Ist schlieBlich die Steuerfunktion durch p(t) fiir die Gesamtoberflache F gegeben, so folgt 


T,(P, t) = p(0) G(P, N+ for ) G(P, t— 1) dr — 9(t) H(P) (25) 
mit 
AF) 0, (26) 
HCP) 
und 
J — = 
D{C(P,i)} =0, G(P,0)=H(P), (27) 
L{G(P,t)},;=90 fir t>0. 
Das stationare Problem (26) hat die einfache Lésung 
H(P) =— 1, 
daher wird 
T(P, t) = 2 A, uP) en?" + 9(0) 5 yal) os 
+ fot) 3 By uf P) HM de + 9) (28) 


wal 


Setzen wir noch T(P, 0) = f(P) = y(0) = T, = konst. (volliger Warmeausgleich zu Beginn des 
Prozesses), dann liefert ein Vergleich von 


(Pe oer (P) Ol il, 


L{T(P,t)}y = fire r= 0 
und 

D{G(P,)} =0, G(P,0)=—1 

LAC(P,)\,;=90 fir t>0 
sofort 


Dee Pi) = — TAP st): 
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Wir erhalten demzufolge nach (28) 
t ARO we 
T(P, 1) = (0) + f'(2) 3 By uP) de (29) 


Das ist aber das Ergebnis von A. Kneschke. 

Damit haben wir eine Fiille von technisch anwendbaren, gesteuerten Warmeleitproblemen 
erfaBt und einer analytischen Behandlung zuganglich gemacht. Die Lésungen gelten sowohl fiir 
Anwarm- als auch fiir Abkiihlvorginge, da die ®,(P, t) diesbeziiglich in keiner Weise eingeschrankt 
werden muBten. Daf die theoretischen Ergebnisse tatsachlich bei der Behandlung praktischer 
Aufgaben von Nutzen sind, werden wir sogleich an den Beispielen erkennen. 


5. Beispiele. a) Temperaturverteilung in einem Zylinder. Es ist die Temperaturvertei- 
lung in einem Vollzylinder mit der Anfangstemperatur f(r, z), dem Halbmesser R und der Lange | 
Av Bectimmens Die Grenzschicht des den Mantel umgebenden Mediums habe die Temperatur @(z, t) 
(Warmeiibergangszahl x). Die Grundflache des Zylinders werde adiatherman abgegrenzt, wahrend 
die Grenzschicht des die Deckflache umgebenden Mediums die Temperatur 4(7) p(t) haben soll 
(Warmeiibergangszahl ’). 

Verwenden wir Zylinderkoordinaten (r, p, z), so liegt entsprechend der Aufgabenstellung fol- 
gendes Randwertproblem vor: 


ot or? r or | @z2 


Us oa le oT 
ees ) = OH) (ze) _.=% 


(Patt) =m, (ZF) =0 


Um dieses zu lésen, berechnen wir zunachst nach (8) T,(r, z, ) aus 


oT ue fee TIES =a 
a — 


oT OT eT. weet, 

a rt A ae Fae ze) an 
eule 20153 () 0% 
a R or neal) 


1 oT, ie {pS () 
(73 aaa On & = 
I Oz 2=0 


Mittels des Bernoullischen Produktansatzes ergibt sich in bekannter Weise 


ea BRod pene 
T(r, od t) = Pia pS, An J(h,,, r) cos k,, Zz as (hin t+ kn)t 
m=l1n=1 


mit 


Ril 
Af ff(r, 2) 7 Jo(hm r) cos ky z dr dz 
00 


ie 2 A é A were 2 : , 
R (e hn) nm J ! + © sin® ky 1 T2(im R) 


Die h,, bzw. k,, sind die Wurzeln der Kigenwertgleichungen 


JER) eed a 


Jy(h R) Ritts ctgkl=- kl. 


Unter Jy) und J, sind die Besselfunktionen nullter bzw. erster Ordnung zu verstehen. Die fiir das 
Folgende erforderlichen Relationen zwischen diesen Funktionen fede man in der Literaturt. 


1G. N. Warns A Treatise on the Theory of Bessel Functions. Cambridge 1952. _ 
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Nun bestimmen wir die stationaren Probleme. Da die Steuerfunktion fiir die Grundflaiche 
identisch Null ist, sind nur zwei derartige Aufgaben zu lésen. Es ergibt sich fiir H,(r, z, t) nach (15) 


eH, ' 1 0H, , @#H, uty, 
Or? r Or | G22 2 
A OH, oH, 
ie @ th)  =— 99), (fe). = 
A OH, oH, 
(2 Bz +H) =0, (a). =° 
und daraus 
H(r, z, ie — Cn(t) Joli lbsie) cos lin Zz 


Te 1 - thn Jy(i kn R) — (Glen R) 


mt. 2 == y— 1 und 
l 
2 { De, t) cos kp z dz 
0 


C,,(t) = 2 r 
1+ — sin? k,, 1 
a 


Die zur Deckflache gehérige Steuerfunktion hat die spezielle Form x(r) p(t). Deshalb hangt H, 
nicht von dem Parameter t ab. Entsprechend (22) folgt fiir H,(r, z) 


OH, , 1 OH, , ®H, _ 0 

Oey Gr ee 
& a aie H,) =0, ie =0, 

or r=R or r=0 
A 0H, aol oH, et 
B de! Hy) hehe ( dz ‘he ae 
Wir erhalten 
22 iD, lee h I m 

RCD ee 7 al mT) cosh (hm 2) 


all = hm sinh (hy, 1) + cosh (hin l) 


mit 
R 
DA GAA edna) ele 
0 


DS te Wey 2 e 
R? (z=) A. | J? (lim R) 


SchlieBlich bestimmen wir Gi (r, z, T, t) nach (16) aus 


a(x 2% 1 aG* a2G* 
“i e eS r ar Z| =i, Gi(r, z, t, 0) = Hy(r, 2, 7) , 
ie OCT ct) 5) a =p 

a Or a 4 re ear ed 

1 eGr aG * fiir t= 0 

=7. 1 = GF =() ° 7 : =0. 

a’ Oz pa Gz ay 


und G,(r, z, t) nach (23) aus 
oe aa & | : oh | on == (() ? G(r, Bs 0) = H,(r, z) ? 


ot or? | ry ar | €2% 

(; 0G, == () laa = 0 

cA topee te sili Or} : 

a ae = fir t>0. 
A2eG; aG = 

faces se (7 ee =0 

(7 Oz +60) eles 
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Es wird 
Gir, 2, t,t) =— = SE Yee) ee rcoslin = (at (hin + hid! 
mat wel West ke (2 hn ++ 7 J(him R) 
und 
(rit) See S Ss Din fry sit fin U Jo(ltm 1) €08 Hen 4a" +) 


SS make ! m pointe 


Wegen (18) ergibt sich somit 


sites: Seis Cr(0) hm Jo(hm r) cos kn —a%(h?,-+-k2)t 
Ty(%)=—R > = oT ee Se = 
m=lin=I1 [h?, ote k2] (z hn) ap 1am R) 


t 
9 { 0 CAT) AmJo(hm 1) cos Kin 2 one (it Ka) poe 


~~ Te r 2 
R g main=l [h24 mall; hn) oe | Ji(hm R) 


oe Cr(t) Jo(t kn r) cos ky z i 
i # thy I(t kn R) — Jo(t kin R) 
a 


und wegen (21) 
eh ES Tp ia hig eal! Jo(him r) cos kin & one (hin 42) 


T,,(r, 2, t) = — 2 y(0) 4 
a a e+] f 2 sine 


t 


co ve) 2 A 2 ew: 
= | o'() Dig Pes it Hod Jinn 0) 008 Men ie le eae 


bishae dl re lt [h2, + k2 f -+- a sin? ken | 
0 04 


ee Din Jo(hm r) cosh (hm %) 


+9) > 
vt 2, us hm sinh (hm lL) -+- cosh (Am 1) 


a’ 


Summieren wir nunmehr die errechneten Temperaturverteilungen entsprechend (19), dann er- 
halten wir die gesuchte Lésung unserer Aufgabe. Es ist also 


Tre, t) ==01 (7, 2,0) i, et) 4 


Man itberzeugt sich ohne Schwierigkeit, daB das Ergebnis in der Tat alle geforderten Bedingungen 
erfillt. 


b) Temperaturverteilung in einem Betonquader wahrend des chemischen Ab- 
bindeprozesses. Es soll die Temperaturverteilung in einem Betonquader mit den Kantenlangen 
2s,, 2s,, 2s; wihrend des chemischen Abbindeprozesses bestimmt werden, wenn die Anfangs- 
temperatur T’, des Betons sowie die Temperatur des umgebenden Mediums konstant ist (Warme- 
ubergangszahl «). 

Ohne Beschrankung der Allgemeinheit kénnen wir annehmen, daB die Temperatur des um- 
gebenden Mediums 0° betragt. Fiir die je Volumen- und Zeiteinheit frei werdende Warmemenge 
des Betons W(x,, x5, x3, t) wird die von K. Hirschfeld’ angegebene Naherungsfunktion 


W= Wea e— tlhe 


verwendet, welche nur von der Zeit abhangig ist. W,,,, und t, sind hier nicht weiter interessierende 
Materialkonstanten. 


 K. Hirschfeld, Die Temperaturverteilung im Beton. Berlin/Gottingen/Heidelberg 1948. 
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Unter Zugrundelegung eines rechtwinkligen Koordinatensystems (x,, %g, Xz) ist somit nach (6) 
folgende Randwertaufgabe zu lésen: 


oT at ert OT. WW, a 
2 : ‘ Pia = max t/t, . : : = 
Ot ‘ & 1 Oxe | ) co : i D2, Xp, %, 0) = Ty, 
A oT = 
( C a 7) a a i : 
OX; x= 3; : 
a Tote 0) (i123) 
ae er Le = () 
a Ox; oi si = 
Wir erkennen unmittelbar, daB 
U(t) = — Ee 5, ete = — Wt, eth 


x 


ein partikulares Integral der inhomogenen Differentialgleichung ist. Entsprechend (7) liefert der 
Ansatz 


TGs Ans Xa. 1) = T(x, Xe, Xs, 0) — Wt, el 
das einfache gesteuerte Warmeleitproblem 


2 2 A2 = 
oT ie oto Fai) = ° Tr Xns ts, 0) Pa OW 4s 


at ax? ' 0x2 ' Ox 
A eT ard 
& Ox; a T). ae — W ty — t/t 


fie 0 est 
AT OL: \ ie 
e Same ate Tes, == W i, é-'h 


Nach (28) lautet die Lésung dieser Aufgabe 


——= ice) ice) ee) . 

Wt —a°(kT + Smt 3n)t 

Gs Xp, Xa, t) =| 1 = > > Aj mn ©O8 hy, %1 C08 ky, X2 cos kz, %3 e 
T, + Wty t=1m=l1n=1 


Dp wd 
ea tlto 9p —V (1p + ham + ks )t 


We é 
eta if W tas 3 ; S A, mn €08 Ky x1 C08 ky, X co8 kg, X3 Fa ehe ey 


| 


5 = es 
+ Wt, eto 
mit 
Sere W ty) sin ky] s, sin kym sz sin ks n $5 
Lee Ae ine lan (s 4 x sin? ky) 5} (+ + a sin?” hams (s + 2 sin? ie | 
Die k,, ke. kz, sind die Wurzeln der Eigenwertgleichungen 


Cs fs, a k, s, (ewe) 


Fir die gesuchte Temperaturverteilung Tey X», Xs, t) ergibt sich damit schlieBlich 


| a? (Ie) i komm ae kSn)¢ 
Cy, X55, b) => S S Aj myn ©O8 Hy) 1 C08 key», Xy C08 ky, X3 € 


SD ge eee 


ice) 


ies) 
° > Ai min ©O8 It, % 608 key, X_ C08 ks, X3 
Lo + W ty mn 


2(k2, 4 ke a 
ew tite — a(k? 2 + km + kz) to ES aaa 
x a? (2, + M2, + 2) t —1 


7 
a ie 


as 
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Wir bemerken noch, daB sich analoge Probleme mit zeitlich und értlich veranderlichen Temes 
peraturen der umgebenden Medien auf Grund der allgemeinen Theorie genau so einfach lésen lassen. 
Damit ist die Méglichkeit gegeben, die bei K. Hirschfeld eingefiihrte Beschrankung auf das Kin- 


dimensionale fallen zu lassen und dariiber hinaus auch verwickeltere Aufgaben zu lésen. 
My 


6. Zusammenfassung. Gesteuerte Warmeleitvorginge, d.h. solche, bei denen verschieden 
Teiloberflichen eines festen Kérpers von Medien umgeben sind, deren Temperaturen sowohl vom 
Ort als auch von der Zeit abhangen und auSerdem im Kérper selbst noch innere Warme entwickelt 
werden kann, sind fiir die Technik von besonderer Bedeutung. Es wird deshalb mittels gewisser | 
EinfluBfunktionen ein Theorem abgeleitet, welches die allgemeinste Aufgabe dieser Art lést. Sowohl 
die EinfluBfunktionen als auch der Lésungsansatz sind physikalisch sehr einfach zu deuten. Aus 
dem allgemeinen Resultat werden dann einige Sonderfalle hergeleitet, unter anderem das Ergebnis © 
von A. Kneschke, der erstmals ein spezielles, gesteuertes Warmeleitproblem auf diese Weise gelést 
hat. Am Beispiel eines endlichen Zylinders und eines Betonquaders ee abschlieBend die Brauch-— 
barkeit der vorgetragenen Theorie gezeigt. 


(Eingegangen am 1. August 1957.) 
Anschrift des Verfassers: Dr. H. Jackel, Colmnitz A47, Kreis Freital (Sachsen). 
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Die Bezieher des ,,Zentralblatt fiir Mathematik“ erhalten die ,,Ergebnisse der angewandten Mathe- 
matik“ gu einem gegeniiber dem Ladenpreis um 10% ermaBigten Vorzugspreis. 


Inhaltsiibersicht: Einleitung. — I. Physikalische und mathematische Grundlagen. Grundgleichungen. 
Linearisierung. Analysis unstetiger Funktionen. — I. Allgemeine Strémungsfelder. Lésungsprinzip und 
Darstellungsformeln. Dickenproblem bei gegebener Geometrie und Auftriebsproblem bei gegebener Druck- 
verteilung des Tragfliigels. Auftriebsproblem bei gegebener Geometrie und Dickenproblem bei gegebener 
Druckverteilung des Tragfliigels,, — III. Strémungsfelder mit konischer Symmetrie. Grundlegende Be- 
ziehungen. Dreieckstragfliigel mit konstantem Anstellwinkel. Dreieckstragfliigel mit nichtkonstantem 
Anstellwinkel. Uberlagerung konischer Strémungsfehler. — IV. Homogene Strémung. Verallgemeinerung 
der Formeln fiir konische Strémungsfelder. Die Behandlung homogener Strémungsfelder mit Unterschall- 
kanten nach dem Singularititenverfahren. — Literaturverzeichnis. — Namen- und Sachverzeichnis. 
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Inhaltsiibersicht: I. Allgemeine Grundlagen. — II. Vereinfachung der Strémungsdifferenzialgleichun- 
gen, Ahnlichkeitsgesetz fiir schallnahe Strémungen. — III. Linearisierte Behandlung schallnaher Stré- 
mungen. — IV. Exakte Lésungen der Potentialgleichung fiir schallnahe Strémungen. — V. Die Grund- 
lagen der Hodographenmethode. — VI. Diskussion schallnaher Strémungsfelder mit Hilfe des Hodo- 
graphen. — VII. Patikularlésungen der Tricomischen Gleichung. — VIII. Str6mungen mit der Machschen 
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Das Verhalten von Gasstrémungen zeigt weitgehende Verschiedenheiten, je nachdem die Strémungs- 
geschwindigkeit kleiner oder gréBer als die Schallgeschwindigkeit ist. Der Ubergang von einer Form zur 
anderen, das ,,schallnahe“ Gebiet, bereitet fiir das Verstiindnis besondere Schwierigkeiten, sowohl yom 
physikalischen als vom mathematischen Standpunkt aus. Jedoch ist in den letzten zwélf Jahren auch 
hier ein gesicherter Kern der Erkenntnis geschaffen worden. Das vorliegende Buch stellt erstmalig die 
Theorie schallnaher Strémungen im Zusammenhang dar. Die Forschung auf diesem Gebiete ist durch 
eine besondere enge Wechselwirkung zwischen physikalischem und mathematischem Denken gekennzeich- 
net. Das bestimmt auch den Charakter dieses Buches. Es zeigt die mathematische Theorie, wie sie aus 
physikalischen Vorstellungen herauswichst. Dabei fihrt es auch in die mathematischen Fragen ein, die 
hier entstehen und die gelegentlich iiber die klassischen Fragenstellungen hinausgehen. Doch wird auch 
dann der Versuch gemacht, durch physikalische Vergleiche die mathematischen Ergebnisse anschaulich 
zu machen. Die Darstellung ist fiir den Aerodynamiker, der einige mathematische Erfahrung hesitzt, 
bestimmt; aber auch dem Mathematiker kann das Buch als Einfiihrung in die Probleme des schallnahen 
Gebiets von Wert sein. 
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